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1.  INTRODUCTION 


Since  Che  60's  an  important  research  effort  has  been  devoted 
to  the  development  of  efficient  structural  optimization  techniques, 
which  could  amplify  the  power  brought  by  modern  numerical  analysis 
tools  such  as  the  finite  element  method. 


Structural  optimization  presents  a large  number  of  facets  for 
which  excellent  reviews  can  be  found  in  Che  references  ( G14,  S6, 
F8,  F13  ] . The  present  work  concentrates  on  one  of  these  aspects, 

that  is  Che  problem  of  minimum  weight  design  of  elastic  structures 
submitted  to  multiple  loading  cases  and  various  constraints. 

It  can  be  presented  as  a non  linear  mathematical  programming  pro- 


blem 


min  WCa^) 


(1.1) 


with  g . (a . ) < 0 

Bj  N j^/ 


where  W denotes  the  weight,  a^  the  design  variables  and  gj  an 


explicit  or  implicit  function  of  Che  design  variables,  eventually 
referred  to  a given  loading  case.  These  functions  describe  a 
number  of  characteristic  responses  of  the  structure  such  as  the 
stresses  in  various  components,  deflections,  dynamic  response, 
stability  limits,  etc  ... 


In  Che  last  two  decades  essentially  two  main  approaches  have 
been  explored  to  solve  this  problem.  The  first  one  is  based  on 
Che  rigorous  application  of  a certain  number  of  the  numerical 
methods  of  non  linear  mathematical  programming.  The  structural 
optimization  programs,  that  have  been  derived,  exhibit  excellent 
convergence  properties,  but  Che  cost  of  their  application  increases 
in  a very  discouraging  manner  with  Che  size  of  the  problems  I K7 , 
C13,  G16,  V3  ] . 


The  second  approach  has  been  developed  from  the  concept  of 
optimality  criteria.  For  a given  type  of  optimization  problem, 
explicit  expressions  of  the  constraints  in  terms  of  Che  design 
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variables  are  derived  which  are  exact  for  isostatic  structures. 
Their  choice  is  guided  by  engineering  intuition.  The  assumption 
is  introduced  that  these  explicit  forms  of  the  constraints  can 
be  used  as  approximations  in  the  case  of  hypeistatic  structures. 
Often  the  initial  problem  (1.1)  disappears  in  its  initial  form, 
and  is  replaced  by  a redesign  algorithm  derived  on  the  basis  of 
the  explicit  approximations  of  the  constraints. 


The  classical  Fully  Stressed  Design  (F.S.D)  procedure  is  the 
best  known  example  of  such  redesign  algorithms.  It  is  exact  for 
an  isostatic  structure^but  leads  eventually  to  non  optimum  designs 
in  hyperstatic  structures  ( S7,  R5,  G14,  K9  ] . However  such  opti- 

mality criteria  reveal,  in  many  practical  situations,  extremely 
good  behavior,  yielding  optimum,  or  near  optimum,  designs  in  a 
number  of  analysis  cycles  which  is  largely  independent  of  the 
number  of  design  variables.  This  attractive  characteristic  led  to 
the  development  of  a second,  distinct,  generation  of  computer 
programs  for  structural  synthesis  [ G15,  DS,  D6,  V4  ] . Their  ap- 

plication cost  is  reasonable,  but  their  convergence  properties 
are  disappointing  in  certain  cases  which  are  difficult  to  identify 
a priori. 


These  two  schools  of  thought  have  been  traditionaly  opposed 
in  the  past.  The  first  contribution  of  the  present  work  is  to 
establish  the  relations  that  exist  between  the  rigorous  mathema- 
tical programming  approaches  and  those  based  on  the  more  intuitive 
optimality  criteria.  The  discussion  concentrates  on  the  classical 
case  where  the  weight  is  linear  in  the  design  variables  and  where 
only  stress  and  displacement  limits  are  imposed,  plus  eventually 
minimum  and  maximum  sizes.  The  extension  of  the  conclusions  to  a 
broader  class  of  problems  is  straight  forward. 


An  essential  step  was  to  recognize  the  nature  of  the  approxi- 
mations that  leads  to  transform  the  original  mathematical  program- 
ming problem  into  an  approximate  one,  from  which  the  optimality 
criteria  derives.  The  next  step  was  to  use  the  understanding  of 
relations  between  the  two  approaches  to  improve  both  of  them. 


1 

M 
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A generalization  of  the  optimality  criteria  is  proposed  |j 

which  consists  in  deriving  them  from  the  KUHN-TUCKER  conditions 
written  for  an  approximate  problem,  in  which  all  the  constraints 
are  approximated  to  the  first  order.  It  is  sh>>wn  that  the  appli-  i| 

cation  of  such  a generalized  optimality  criterion  is  equivalent  jj 

to  solving  the  approximate  problem  exactly.  An  efficient  way  is  | 

provided  by  the  dual  algorithms  of  convex  mathematical  program-  I 

ming.  I 

Hybrid  optimality  criteria  are  defined  as  those  based  on  >1 

....  1 1 

explicit  approximations  of  the  constraints  that  are  of  different  | 

orders.  They  are  especially  attractive  in  the  presence  of  a large 

number  of  stress  constraints  where  a complete  linearization  is  j 

i 

costly  and  often  not  necessary. 

i 

Turning  to  the  mathematical  programming  approaches , it  is 
shown  that, in  the  space  of  the  inverse  of  the  design  variables, 
the  generalized  optimality  criteria  approach  can  be  viewed  as  j 

equivalent  to  a special  form  of  the  linearization  methods  applied 
to  the  original  problem.  This  allows  a better  understanding  of  | 

the  reasons  of  the  divergence  sometimes  met  with  optimality  crite-  : 

i : 

ria.  ^ 

On  the  other  hand,  a strict  projection  primal  algorithm  of 
mathematical  programming  is  derived  in  the  space  of  the  inverse 
of  the  design  variables.  This  space  is  advantageous  as  there 
the  constraints  are  very  shallow.  If  the  first  order  explicit 
approximations  of  the  constraints  are  used,  the  projection  algo- 
rithm can  take  different  forms.  When  the  explicit  approximations 
are  reevaluated  after  each  projection  step,  the  algorithm  is  a 
rigorous  application  of  the  primal  mathematical  programming  tech- 
niques. When  the  explicit  approximations  are  kept  unchanged  until 
the  exact  solution  of  the  approximate  problem  is  reached,  the 
method  is  strictly  equivalent  to  the  application  of  the  generalized 
optimality  criteria. 
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By  selecting  a limited  number  of  projection  steps  before 
recomputing  the  explicit  approximations  of  the  constraints, 
which  implies  to  analyse  the  structure,  the  projection  method 
becomes  a mixed  method  with  properties  lying  between  those  of 
the  optimality  criteria  approaches  and  those  of  pure  mathemati- 
cal programming.  Moreover  this  concept  of  mixed  method  is  easily 
extended  to  situations  where  the  order  of  approximation  of  the 
constraints  is  variable,  like  in  the  hybrid  optimality  criteria. 

The  efficiency  of  the  proposed  algorithms  is  proved  on  a 
rather  large  number  of  applications  taken  from  the  litterature 
on  structural  optimization. 

The  conclusions  of  the  present  invest igat i cn  are  the  follo- 
wing 

- efficient  algorithms  for  structural  optimization  seem  to  result 
only  from  the  rational  combination  of  mathematical  programming 
techniques  and  engineering  approximations; 

- the  approximations  of  the  constraints  introduced  in  the  deriva- 
tion of  optimality  criteria  are  quite  often  perfectly  satisfac- 
tory, even  those  based  on  F.S.D.,  provided  they  are  correctly 
used  in  the  weight  minimization  algorithms; 

- simple  tests  allow  to  detect  a priori  the  situations  where 
more  elaborate  approximations  are  necessary; 

- the  most  efficient  methods  seems  to  be  based  on  mixed  order 
approximations,  corresponding  to  what  is  defined  in  this  work 
as  hybrid  optimality  criteria.  The  solution  of  the  correspon- 
ding approximate  problems  are  best  obtained  by  dual  methods  of 
mathematical  programming  when  the  number  of  constraints  is 
small  and  when  their  non  linearity  is  not  too  pronounced. 
Otherwise  the  primal  methods  are  to  be  preferred; 

- if  the  mixed  approximations  improve  significantly  the  conver- 
gence properties  of  the  optimality  criteria,  it  remains  that 
their  convergence  can  never  be  taken  for  granted,  as  they  are 
equivalent  to  linearization  methods  of  mathematical  program- 
ming; 


1 

) 

when  the  hybrid  optimality  criteria  fail,  the  combination  of  ] 

the  new  mixed  method,  based  on  primal  projection  algorithms,  ! 

and  of  various  explicit  approximations  of  the  constraints, 
yields  a very  flexible  and  efficient  tool. 


2.  THE  TWO  CLASSICAL  APPROACHES 

As  mencionned  in  the  introduction  the  present  work  concentra- 
tes on  the  solutions  of  the  minimum  weight  design  problem  under 
three  types  of  constraints,  that  are , limitations  in  the  sizes 
of  the  design  variables  themselves,  limitations  in  the  stresses, 
and  limitations  in  deflections  or  flexibilities.  The  structures 
are  supposed  to  be  composed  only  of  bars  and  plates  in  a plane 
state  of  stress.  Their  cross  sections  or  thicknesses  are  the  de- 
sign variables  a^  and  the  objective  function,  that  is  the  weight 
W,  is  a linear  function  of  them.  The  problem  can  be  presented 
mathematically  in  the  form 


min  U (a.)  ■ I p.  1.  a. 

1 j i 1 1 


(2.1) 


with  the  constraints 


a . < a . < a . 
—Ill 


1,  n 


(2.2) 


iroblem 


u . . (a)  < u . 

2 


- 1,  n^ 

- 1. 


(2.3) 


< ‘’k 


1 , n 

1.  n. 


(2.4) 


where  n denotes  the  number  of  design  variables,  n^  the  number 
of  flexibility  constraints,  n^  the  number  of  loading  cases, 
the  mass  density,  a geometrical  quantity,  length  of  the  bars 
or  surface  of  the  plates,  such  that  a^  is  the  volume  of  the 
element,  u.,  denotes  a generalized  displacement  under  the 
loading  case  and  a stress  component, or  an  equivalent  stress, 

in  the  k^^  element  under  the  loading  case.  The  upper  (lower) 

bar  denotes  an  upper  (lower)  limit  prescribed. 


This  problem  is  labelled  problem  R (for  Reference  or  Real 
problem).  It  is  a non  linear  ^Mathematical  programming  problem. 
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Two  classical  approaches  have  been  used  to  solve  it. 

The  first  is  based  on  more  or  less  straight-forward  applications 
of  the  numerous  and  various  minimization  techniques  that  have 
been  developed  in  mathematical  programming.  The  second  class 
of  approaches  is  more  intuitive  and  replaces  the  problem  R by  a 
sequence  of  simpler  problems  obtained  by  introducing  various 
approximations  derived  by  taking  advantage  of  the  physical  or 
intuitive  knowledge  of  the  behaviour  of  the  structure. 

In  this  chapter  a brief  review  is  presented  of  some  of  the 
techniques  used  in  these  two  types  of  approaches.  It  is  not 
intented  to  be  exhaustive^nor  to  give  the  detailed  formulations 
that  can  be  found  in  the  references.  The  purpose  of  the  review 
is  simply  to  recall  some  aspects  of  these  methods  that  are  useful 
for  understanding  the  arguments  developed  in  the  next  chapters. 

2.1.  Direct  approaches  based  on  mathematical  programming 

2.1.1.  Classification  of  the  non  linear  mathematical  programming 
methods 

In  the  following  classification  of  the  methods  of  non 
linear  M.P.,  only  the  most  classical  methods  have  been  retained 
that  can  be  applied  in  the  context  of  the  present  work. 

The  statement  of  the  problem  R (2.1  to  2.4)  involves  a linear 
objective  function  and  non  linear  contraints.  This  problem 
however  can  be  transformed  in  different  manners, so  that  a much 
broader  class  of  methods  has  to  be  considered  chan  chose  that 
are  especially  adapted  to  the  type  of  functions  and  constraints 
initially  considered.  The  classification  adopted  here  is  based 
on  the  nature  of  Che  constraints. 

Unconstrained  minimization 

Almost  all  of  these  methods  follow  the  same  basic  scheme  which 
consists  in  defining  a search  direction  and  Chen  in  minimizing 


Ithe  function  along  Chat  direction, which  implies  the  computation 

of  a progression  step.  The  search  for  the  optimal  step  Chat 
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minimizes  the  objective  function  along  that  direction  is  called  | 

a linear  search.  When  this  step  has  been  achieved  a new  search  { 

direction  is  defined  and  the  process  is  repeated  until  conver-  { 

I 

gence  eventually  occurs.  An  important  concept  in  the  definition  | 

1 

of  the  search  directions, is  that  of  conjugate  directions 

( HI  , B1  ] . Although  some  methods  exist  that  do  not  require  ; 

the  computation  of  any  derivative  of  the  function,  like  the  I 

POWELL  method  [ P2  ] , the  only  choice  in  our  context  is  between 

methods  that  require  an  evaluation  of  the  first  or  second  deri-  j 

vatives  of  the  function.  Among  those  that  use  only  the  first 
derivatives, we  find  the  various  gradient  methods  like  that  of  | 

the  steepest  descent  [ Cl,  C2,  G1  ] , the  conjugate  gradient 

methods  [ HI,  H3,  Bl,  F4  ] and  the  variable  metric  methods  i 

j 

that  are  especially  useful  when  the  function  has  a pronounceed  | 

excentricity  [ F6  ] ; the  DAVIDON-POWELL  is  a typical  example  | 

[ PI,  FS,  B2,  P4,  P5  ] . When  the  second  derivatives  of  the  ] 

1 

objective  function  can  be  obtained,  powerful  methods  can  be  cons-  j 

tructed  by  using  quadratic  approximations  of  the  function.  ] 

ii 

Various  generalizations  of  the  NEWTON  method  have  been  proposed  j 

J 

[ C4,  G2,  L2,  M2  ] . | 

Minimization  with  linear  constraints 

The  minimization  of  a function  under  a set  of  linear  equality 

. 

or  inequality  constraints  plays  a foundamental  role  in  M.P. 

Even  if  one  does  not  consider  linear  programming  .which  is  an 
important  special  case  but  not  relevant  to  our  problem,  the 
fact  is  that  many  methods  of  non  linear  M.P.  lead  to  a sequence 
of  problems  with  linear  constraints.  The  classical  methods  of 
unconstrained  minimization  are  easily  transposed  to  this  type 
of  problem  from  the  point  of  view  of  the  choice  of  the  search 
directions  and  of  the  convergence  properties.  An  essential  modi- 
fication in  the  case  of  inequality  constraints  is  the  addition 
of  an  algorithm  which  decides  when  a contraint  becomes  active 
or,  inversely,  must  be  abandoned  . 

An  important  class  of  methods  is  constituted  by  the  projec- 
tion methods.  They  can  be  characterized  by  the  fact  that  they 

i 

— — — : d 


progress  along  descent  directions  that  remain  on  the  boundaries 
of  the  polyhedron  of  the  feasible  domain.  They  present  therefore 
the  advantage  of  producing  a sequence  of  feasible  points  that 
correspond  to  decreasing  values  of  the  objective  function. 

An  example  of  such  methods  is  the  gradient  projection  method  of 
ROSEN  [ R1  ] which  derives  from  the  steepest  descent  method. 

Better  convergence  is  obtained  by  using  projected  conjugate 
gradient  methods  [ M3,  F4  ] or  projected  variable  metric  methods 
[ G4 , G6  ] . These  methods  require  only  the  knowledge  of  the 

first  derivatives  of  the  functions.  Second  order  projection 
methods  can  also  be  applied  when  the  second  derivatives  of  the 
objective  function  are  available.  They  are  obtained  by  a genera- 
lization of  the  NEWTON  method  f G3,  G5  ] . 

Minimization  with  non  linear  constraints 

In  this  more  general  case,  two  essentially  different  classes  of 
methods  have  to  be  distinguished.  The  first  groups  the  primal 
methods  which  attempt  to  solve  directly  the  original  problem, 
i.e.  in  our  case  the  problem  R (2.1  to  2.4).  They  are  characte- 
rized by  an  iterative  process  that  produces  a sequence  of 
feasible  points  corresponding  to  decreasing  values  of  the  objec- 
tive function.  This  characteristic  is  important, as  it  allows 
to  stop  the  process  before  reaching  the  optimum  and  still  yield 
an  improved  design.  Among  the  primal  methods  one  can  distinguish 
between  two  subclasses.  The  first  one  uses  the  concept  of 
feasible  directions  introduced  by  ZOUTENDIJK  ( Z2  ] . The  feasible 

directions  are  selected  in  such  a way  that  none  of  the  constraints 
are  violated  during  the  progression.  The  second  one  consists  in 
a further  generalization  of  the  projection  methods  developed 
for  linear  constraints.  In  these  methods, each  minimization  leads 
to  a moderate  violation  of  the  active  constraints  and  has  there- 
fore to  be  followed  by  a restoration  phase  that  brings  back  the 
design  point  on  the  composite  constraint  surface.  First  order 
methods  are  available  like  the  ROSEN's  gradient  projection 
method  [ R2  ] , or  the  conjugated  gradient  methods  [ A2  ] , as  well 

as  second  order  methods  generalizing  the  NEWTON  method  ( S3,  S4  ] . 
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A difficulty  Chat  is  common  to  all  of  these  primal  methods 
comes  from  the  need  to  start  the  process  at  a feasible  point. 

When  the  constraints  are  non  linear, this  is  sometimes  far  from 
being  evident. 

In  the  second  group  of  methods,  the  original  problem  is  not 
directly  solved  but  transformed  into  another  problem  or  a 
sequence  of  auxiliary  or  approximate  problems.  Among  these 
methods  ve  find  the  dual  methods,  the  linearization  methods  and 
the  penalty  methods.  An  advantage  of  these  methods  is  that  the 
non  linearity  of  the  constraints  does  not  introduce  any  special 
difficulty.  Their  disadvantage  comes  eventually  from  the  fact 
that  the  design  points  that  are  produced  after  each  iteration 
are  not  necessarily  feasible  points, and  therefore  it  is  usually 
not  possible  to  stop  the  process  before  reaching  the  optimum. 

It  is  also  worth  mentionning  that  many  of  these  methods  require 
certain  convexity  properties  which  are  not  necessary  in  the 
primal  methods.  The  problem  R,as  it  stands, is  not  necessarily 
convex, and  therefore  optimization  methods  requiring  the  convexity 
of  the  design  space  can  only  be  applied  to  a modified  problem. 

The  simplest  method  to  solve  a non  linear  M.P.  problem  is 
perhaps  to  transform  it  in  a sequence  of  linear  programming 
problems.  Unfortunately  this  procedure  can  converge  only  if  a 
local  minimum  exists  at  a vertex  of  the  design  space  which  is, 
in  practice,  very  unusual.  A certain  number  of  more  elaborate 
linearization  methods  have  been  developed  to  avoid  this  diffi- 
culty. Among  them  we  shall  mention  the  KELLEY  cutting  plane 
method  [ K4 , W1 , LI  ] which  however  applies  only  to  convex  pro- 
blem. The  methods  of  approximation  programming  of  GRIFFITHS  and 
STEWART  [ G7  ],of  WILSON  [ W4  ] and  of  BEALE  ( B3  ],  introduce 
various  improvements  like  the  definition  of  move  limits,  or 
the  linearization  of  the  constraints  only, while  a quadratic  ap- 
proximation of  the  objective  function  is  used.  An  important  dif- 
ference with  the  primal  projection  me thods, which  at  the  current 
point  also  use  a linearization  of  the  constraints,  is  due  to 
the  fact  that  in  linearization  methods  the  design  point  obtained 
after  solving  an  approximate  linearized  problem  is  not  moved  back 


y 
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onto  Che  exact  composite  constraint  surface. 

The  penalty  methods  consist  in  transforming  Che  initial 
problem  in  a sequence  of  unconstrained  problems  by  adding  to 
Che  objective  function  an  auxiliary  function  chat  describes 
the  degree  of  non  satisfaction  of  the  constraints  [ FI  ] . 

The  choice  of  Che  auxiliary  function  can  be  made  such  chat 
Che  sequence  of  design  points  remain  inside  or  outside  Che 
feasible  domain  [ FI,  C5,  Fll,  A3,  R3,  LI  ] . 

In  Che  special  case  of  convex  programming,  Che  lagrange 
multipliers , that  are  associate  Co  the  constraints , have  Che 
meaning  of  dual  variables  in  terms  of  which  a new  formulation 
of  Che  problem  can  be  established.  This  dual  problem  consists 
in  maximizing  the  lagrangian  function  in  terms  of  the  dual 
variables  with  the  only  constraint  that  these  dual  variables 
have  to  remain  non  negative  [ FI,  LI,  W2,  F12  ] . These  methods 

play  an  important  role  in  the  context  of  optimality  criteria 
and  will  be  discussed  in  more  details  at  that  occasion. 


2,1,2.  Gradients  of  the  constraints 

In  most  of  the  efficient  minimization  techniques  it  is 
necessary  to  compute  the  value  of  the  various  functions  that 
define  the  problem  and  their  gradients.  As  Che  flexibility  and 
stress  constraints  are  implicit  functions  of  the  design  variables 
their  gradients  can  only  be  computed  numerically. 


In  Che  context  of  the  finite  element  method,  a flexibility 


u . 
J 


can  be  expressed  as  a linear  combination  of  the  nodal  D.O.F. 


corresponding  to  the  displacements  q 


Uj  - b.  q (2.5) 

where  b.  is  a line  matrix  of  constants.  The  gradient  of  Che 
flexibility  under  a loading  case  t is  therefore 

I 


(2.6) 
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For  Che  stress  constraints,  a stress  component  can  be  written 
in  the  same  manner  as  a linear  combination  of  the  D.O.F.  of 
Che  corresponding  finite  element,  using  in  practice  simply  a 
line  of  Che  stress  matrix  c 


'kf 


- t. 


and 
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kf 


da 


(2.7) 


If  Che  constraint  is  on  an  equivalent  stress,  Che  expression 
of  its  gradient  can  be  obtained  from  Che  gradients  of  the  3 
stress  components. 

The  computation  of  the  gradients  of  the  two  types  of  constraints 
reduces  to  the  evaluation  of  Che  gradients  of  the  generalized 
nodal  displacements  of  the  finite  element  model.  The  method  of 
"pseudo-loads"  [ F14,  GIO  ] is  rather  advantageous. 

Considering  the  equilibrium  equations 
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for  the  loads 


one  has 


0 


I ■ 1 , n 
' c 


i - 1 , n 

t - 1,  n^ 


from  which 


q 


i 


(2.8) 


(2.9) 


(2.10) 


In  the  elements  that  are  considered  in  this  work,  the  stiffness 
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matrices  are  linear  in  the  design  variables,  and  can  be  written 

K.  - a.  K.  (2.11) 

where  is  the  stiffness  matrix  for  a unit  value  of  the  design 
variable.  The  pseudo  loads  are  defined  by 

gi,  i - 1.  n (2.12) 

1 ■ 1 , n 
’ c 

so  that  the  computation  of  the  gradients  (2.10)  requires  only 

the  application  of  the  (n  x n^)  pseudo  loads  to  the  structure, 

in  addition  to  the  set  of  n actual  loads. 

c 

This  can  be  handled  rather  efficiently  by  the  finite  element 
programs, especial ly  if  the  total  number  of  loads  is  not  too 
large . 

2.1.3.  Characteristics  of  some  of  the  M.P.  approaches 

A large  variety  of  solutions  has  been  proposed  to  the 
structural  synthesis  problem  that  are  based  on  one  or  the  other 
of  the  M.P.  methods  mentionned  in  the  classification  given  above. 
The  earliest  algorithms  were  based  on  primal  methods, like  the 
alternate  step  methods  proposed  by  SCHMIT  [ S7  ] , GELLATLY  et  al. 

( G9,  GIO,  Gll,  G12  ] , the  ROSEN's  projected  gradient  method 

applied  by  BROWN  [ B4  ] , the  ZOUTENDIJK's  feasible  direction 

method  used  by  KARNES  [ K7,  T2  ] and  VANDERPLAATS  [ VI  ] etc. 

All  of  these  methods  solve  directly  the  problem  R in  its  form 
given  by  (2.1)  to  (2.4),  that  is  in  terms  of  the  design  variables. 
The  linearization  methods  have  been  applied  by  MOSES  [ M6  ] , 

REINSCUMIOT  [ R4  ] , SCHMIT  i S8  ] and  others.  It  is  worth 

mentionning  that  a certain  number  of  these  authors  have  noticed 
the  advantages  of  using  the  inverses  of  the  design  variables  as 
unknowns.  The  advantage  comes  from  the  fact  that  the  stress  and 
flexibility  constraints  are  linear  in  the  inverses  of  the  design 
variables  when  the  structure  is  isostatic , and  rather  shallow 


J 
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when  hyperscatic, so  that  a linear  approximation  of  the  cons- 
traints corresponds  to  a good  approximation  in  this  design 
space.  The  penalty  methods  have  been  applied  by  SCHMIT  [ SIO  ] , 

MARCAL  [ M9  ] , KAVLIE  [ K8  ] and  MOE  [ MIO  ] . 

All  of  these  algorithms  share  certain  common  characteristics. 

One  of  them  is  the  need  for  quite  frequent  evaluation  of  the 
derivatives  of  the  objective  function  and  of  the  constraints. 

These  can  only  be  obtained  after  a reanalysis  of  the  structure! 
which  is  a costly  operation.  Moreover  the  number  of  iteration 
steps  increases  with  the  number  of  design  variables!which  is  a 
well  known  behaviour  of  all  the  constrained  or  unconstrained 
minimization  algor ithms , when  only  the  first  derivatives  are  used. 

So  not  only  the  cost  of  the  analyses  increases  when  more  complex 
structure  are  optimized, but  also  the  number  of  analyses  and  of 
pseudo  load  cases.  In  practice, none  of  these  methods  revealed  | 

able  to  solve  large  problems  at  a cost  competitive  with  optimality  | 

.1 

criteria  based  methods.  However, the  main  advantage  of  the  H.P.  ] 

i 

methods  is  their  generality,  their  convergence  behaviour  which  i 

can  be  predicted  and  often  guarantied  due  to  the  sound  mathema-  | 

tical  bases  of  the  methods. 

An  important  exception  has  been  brought  recently  by  the  | 

method  of  SCHMIT  and  MIURA  [ S18,  1976  ],  which  is  characterized 

by  the  use  of  linear  explicit  approximations  of  the  constraints  i 

in  the  space  of  the  inverses  of  the  design  variables. 

Each  approximate  problem  is  solved  using  either  a feasible 
direction  method  or  a penalty  method.  This  method  reveals  very 
efficient, but; at  the  light  of  the  present  research, should  no 
longer  be  considered  strictly  as  a M.P.  method.  Due  to  the  type 
of  approximation  introduced.  It  becomes  in  fact  a mixed  method 
and  as  such  will  be  examined  later  on. 
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2.2.  Indirect  approaches  baaed  on  optimality  criteria 


2.2.1.  Problems  with  stress  constrsints 


When  only  stress  constrsints  ere  present  in  sddition  to 
the  side  constraints  (2.2)  the  most  popular  approach  consists  in 
using  the  F.S.D.  concept.  After  each  analysis  the  structure  is 
resized  according  to  the  stress  ratio  algorithm 


v+1 
1 . 
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a . max 
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v+1 

a . 
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> a . 
— 1 


(2.13) 


where  v denotes  an  iteration  number.  It  is  well  known 
[ Mil,  M12,  C7,  R5,  K9  ] that  the  approach  leads  to  an  optimum 
in  the  case  of  an  isostatic  structure.  In  such  a case, the  stress 
constraint  becomes  explicit  and  the  problem  R is  replaced  by  the 
problem 


min  W • £ p . I . a . 

1 1 


(2.14)  I 


Indeed  the  application  of  the  F.S.D.  does  not  require  to  explicit 
the  objective  function.  The  problem  (2.14)  being  a linear  pro- 
gramming problem,  its  solution  is  necessarily  at  a vertex  of  the 
design  space  which  is  unique  and  corresponds  to 


it  'x, 

a . ■ max  (a . , a . ) 

1 —1  ’ 1 ' 


(2.15) 
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In  Che  case  of  hyperstacic  structures,  the  solutions 
obtained  using  the  F.S.D.  are  no  longer  necessarily  minimum 
weight  designs  [ B7,  B8,  K9,  R5 , SI2,  S13  1 and  some  difficul- 
ties in  convergence  are  observed  in  certain  cases  [ B8,  M13,  CI3, 

R4t,  G14  ] . In  many  practical  situations  however  the  F.S.D. 

leads  very  rapidly  to  a design  that  is  very  close  to  the  minimum 
and  is  characterized  by  a number  of  iterations  (reanalyses) 
that  is  independent  of  the  number  of  Che  design  variables. 

2.2.2.  Problems  with  flexibility  constraints 


A large  number  of  publications  deals  with  the  optimality 
criteria  approaches  in  Che  case  where  flexibility  constraints 
are  imposed  ( B9,  B7,  PlO,  P13,  S14,  S15,  P14,  BIO  ] . 

Considering  a virtual  loading  case  g.  associated  to  each 
flexibility  constraint  u^ , its  virtual  work  is  written 


T T T ^ 

Uj  - qj  g - q gj  - q.  K q 


(2.16) 


where  q denotes  the  displacement  vector.  It  is  expressed  as  the 
sum  of  Che  contributions  of  each  element 


u . ■ I q . . K . q . 
J i ^ ^ 


1,  n 


(2.17) 
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1 1 


(2.18) 


with  c . . 

ij 


(q  . . K . q . ) a . 

1 1 


(2.19) 


The  flexibility  coefficients  are  constant  for  an  isostatic 
s c rue  Cure, which  is  the  situation  considered  to  derive  Che  opti- 
mality criteria.  Taking  the  case  where  only  one  flexibility 
constraint  is  imposed, which  is  consequently  an  i^quality  constraint, 
it  is  classical  Co  separate  the  n elements  into  a group  of  active 
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members  (i  ■ 1,  and  a group  of  passive  members  (i  - n) 

so  that  the  flexibility  can  be  expressed  as 


u,  ♦ u ■ u 
1 2 


"l  c. 

r -i 

i *i 
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n c . 

E -i 

Hj^  + l ®i 


(2.20) 


With  equality  constraints,  a lagrangian  method  can  be  used 
which  is  based  on  the  augmented  functional 


£ (a^,r)  • W(a^)  ♦ r I u(aj^)  - u ] 


(2.21) 


Its  stationarity  conditions  are  given  by  the  non  linear  system 
of  equations 

i - 1,  (2.22) 

1 c . _ 

E ^ - u - u (2.23) 

1 *i  2 

In  (2.22)  the  coefficients  c^  have  to  be  positivej which  yields 
a mean  to  distinguish  between  the  active  and  passive  members. 
Substitution  of  (2.22)  into  (2.23)  yields  the  value  of  the 
lagrange  multiplier 

“l  1/2  2 
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(2.24) 
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Reintroducing  (2.24)  in  (2.22)  yields  the  new  values  of  the 
active  design  variables  (i  >■  1,  n^) 


, 1 1/2  c.  1/2 

*i  ■ ^ ^ *^k^  . t . ^ 

u-u 2 1 11 


(2.25) 


This  redesign  formula  is  applied  after  each  analysis  cycle  in 
the  case  of  a hyperstatic  structure  and  the  choice  between 
active  and  passive  member  is  also  repeated  at  each  iteration. 
The  relation  (2.22)  can  be  written 


e . ■ — ■ z with  e.  ■ — 

1 r p . X . a . 1 a ■ 

111  1 


(2.26) 


The  optimality  criterion  states  that  the  minimum  weight  design 
is  obtained  under  one  flexibility  constraint  when  the  virtual 
strain  energy  density  e^^  per  unit  of  mass  is  the  same  for  all 
the  active  elements. 

The  criterion  is  exact  for  isostatic  structure  and  approximate 
when  hyperstatic  [ BIO,  G13  ] . The  inclusion  of  lower  limits 

to  the  design  variables  is  a straight^forword  addition. 

2.2.3.  Problems  with  stress  and  flexibility  constraints 

Such  problems  have  received  different  treatments.  In  the 
approaches  pioneered  by  BERKS  and  GELLATLY  [ G13,  Bll  ] the 
flexibility  constraints  are  considered  as  the  principal  ones, 
while  the  others  (stress  and  lower  limits)  are  treated  as  side 
constraints.  The  explicit  problem  solved  after  each  analysis  is 


problem 


with 


£ p . t . a . 

1 1 1 


c . . 

£ < u. 

^ J 


i ■ 1 , n 


j ■ 1.  n, 


(2.27) 


a . > a . 
1 1 
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where  a*  is  given  by  (2.15).  In  an  isostatic  structure  the 
problem  R (2.1)  is  identical  to  the  problem  A, while  in  the 
hyperstatic  case.it  is  expected  that  a sequence  of  approximate 
problems  A yields  to  the  solution  of  the  exact  problem  R. 

In  fact, most  of  the  classical  approaches  do  not  solve  directly 
the  problem  A, but  start  with  its  simplified  version 


min  W > £ p . I . a . 

1 1 1 

with  u . ■ u . 

J J 


(2.28) 


where  one  has  dropped  the  side  constraints  and  supposed  that 
the  flexibility  constraints  are  all  active  at  the  optimum. 
Taking  the  side  constraints  into  account  and  selecting  the 
active  flexibility  constraints  are  considered  as  auxiliary  pro- 
blems. The  BERKE  optimality  criterion  (2.22),  (2.23)  is  still 
applicable  and  can  be  generalized  for  multiple  flexibility 
constraints  by  introducing  a number  of  lagrange  multipliers 
equal  to  the  number  of  active  constraints  [ Bll  ].  It  yields 
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(2.29) 


(2.30) 


Unfortunately  it  is  no  longer  possible  to  solve  analytically 
this  non  linear  system  of  equations  in  the  case  of  multiple 
constraints.  The  envelope  method  of  GELLATLY  and  BERKE  ( G13  ] 

furnishes  an  approximate  solution  to  this  problem.  They  determine 
for  each  active  variable  the  flexibility  constraint  that  gives 
it  the  maximum  size  by  applying  the  relation  (2.25)  for  each 
flexibility.  The  resizing  relations  become 
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i - 1 , n, 


The  selection  between  active  and  passive  variables  is  different 
for  each  constraint  and  yields  to  different  values  of  U2  denoted 
U2  • The  minimum  size  of  variables  is  given  by  a*  defined  in 

(2.15).  In  the  case  of  multiple  constraints  this  procedure 
does  not  yield  the  exact  so lut ion , even  in  an  isostatic  problem. 
It  can  be  interpreted  as  a pseudo  lagrangian  method  where  the 
lagrange  multipliers  are  approximated  by 
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(2.32) 


and  where  the  design  variables  are  given  by  the  envelope 
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a.  • 1—  max(r.  c..) 

1 Pi^i  j ' J ij' 


i - 1 , n, 


(2.33) 


instead  of  the  exact  express  ion, solution  of  (2.29). 

Another  approximate  solution  has  been  proposed  by  TAIG 
and  KERR  (13).  They  have  used  an  exact  lagrangian  method 
but  in  terms  of  the  inverse  of  the  lagrange  multipliers 


(2.34) 


so  that  the  optimality  criterion  (2.26)  becomes 
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(2.35) 
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The  solution  of  the  simplified  minimization  problem  (2.28)  is 
given  by  the  system  of  (n  + n^)  non  linear  equations 
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I 
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(2 . 36) 
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(2.37) 


The  solution  of  (2.36)  is  obtained  by  an  iterative  NEWTON- 
RAPHSON  procedure.  A difficulty  arises  by  the  fact  that 
initial  values  for  the  X^  are  necessary  to  initiate  the  pro- 
cedure,which  are  sometimes  difficult  to  find.  In  the  case  of 
one  flexibility  constraint  the  TAIG  and  KERR  approach  yields 
results  identical  to  that  of  the  GELLATLY-BERKE  envelope  method. 

In  both  methods, most  of  the  difficulties  that  appear  in 
the  applications  can  be  related  to  the  fact  that^ instead  of 
solving  the  problem  A,  they  solve  the  problem  (2.28)  with 
equality  constraints. 

Ocher  solutions  have  been  proposed  [ V4,  B8,  B12  ] which 

proceed  from  similar  points  of  view.  They  are  not  detailed  here. 
The  two  methods  of  BERKE  [ BIO  ] and  TAIG  and  KERR  [ T3  1 are 
sufficiently  characteristic  of  the  optimality  criteria  approach 
to  allow  to  develop  the  comparisons  that  are  essential  for  the 
present  work. 


2.3.  Mixed  approaches 

Mixed  approaches  are  characterized  by  the  fact  that  they 
rely  simultaneously  on  the  minimization  algorithms  of  mathema- 
tical programming  and  on  Che  intuitive  approximations  of  the 
optimality  criteria. 


Such  approaches  have  been  proposed  by  VENKAYYA  [ V2,  V3  ) and 

others  [ 05,  S17  ] . We  note  especially  in  the  context  of  the 

present  work, the  method  of  "high  quality  explicit  approximations" 
of  the  constraints  proposed  by  SCHMIT  and  MIURA  ( S18,  S20,  S8  ) . 

Although  in  the  spirit  of  the  authors  the  method  pertains  to  the 
line  of  mathematical  programming  methods,  it  will  be  demonstrated 
later  that  it  is  indeed  a mixed  one.  The  approximations  of  the 
constraints  used  by  SCHMIT  are  obtained  by  a Taylor  serie  expansion 
limited  to  the  first  order  in  the  inverse  of  the  design  variables. 
The  approximate  form  of  the  problem  R that  they  obtain  is 


min  W(x.)  ■ Z p.  — with  x.  • — 

1 . 1 1 x . 1 a . 
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where  the  superscript  ( ) denotes  the  values  computed  at  a 
design  point  where  the  structure  is  analyzed. 

The  problem  R is  replaced  by  a sequence  of  such  approximate 
problems  which  are  solved  by  either  a feasible  directions 
method  or  a penalty  method. 
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3.  A GENERALIZATION  OF  THE  OPTIMALITY  CRITERIA 


An  essencial  ingredient  of  the  optimality  criteria  is  the 
use  of  explicit,  and  hence  most  often  approximate,  expressions 
of  the  constraints  in  terms  of  the  design  variables. 

These  expressions  are  exact  in  the  isostatic  case  only.  In 
hyperstatic  structures  the  degree  of  approximation  is  essentially 
different  for  stress  and  displacement  constraints  in  the  classi- 
cal approaches.  The  first  generalization  proposed  here  consists 
in  presenting  a way  to  derive  optimality  criteria  that  are  based 
on  approximations  that  are  consistent  for  all  the  constraints. 

It  is  to  be  noted  that  this  does  not  imply  necessarily  the  sugges* 
tion  to  use  these  generalized  optimality  criteria  in  that  form, 
but  they  represent  an  essential  step  for  the  understanding  of 
the  relations  between  the  various  structural  synthesis  techniques 
and  for  the  derivation  of  efficient  mixed  methods. 


3.1.  Scaling  of  the  design  variables 


Before  discussing  the  generalization  of  the  optimality 
criteria,  it  is  necessary  to  recall  a property  of  the  scaling 
operation  consisting  in  multipliing  each  design  variable  by 
the  same  factor.  Such  an  operation  evidently  does  not  introduce 
any  redistribution  of  the  internal  forces  if  the  stiffness  ma- 
trices are  linear  in  the  design  variables , what  we  have  assumed 


in  the  present  work.  If  a?  denotes  a point  in  the  design  space, 
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the  scaling  by  a factor  f yields  a new  point  a. 
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at  which  the  displacements  and  stresses  are 
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All  Che  points  obtained  by  (3.1)  are  on  Che  line  joining  a? 

Co  Che  origin  in  Che  design  space.  It  is  denoted  D(a?).  The 
intersections  of  that  line  with  Che  constraint  surfaces  yield 
Che  scaling  factors 
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u . 
J 


(3.3) 


that  are  to  be  applied  to  bring  a design  point  a?  on  those 
surfaces . 


3.2.  Order  of  approximation  of  the  constraints 

The  stress  ratio  redesign  algorithm  used  in  F.S.D. 
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which  is  applied  after  an  analysis  yielding  Che  stresses  o° , 
corresponds  to  an  explicit  approximate  expression  of  the  stresses 
in  terms  of  Che  design  variables 
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(3.5) 


where  the  denotes  an  approximate  expression.  The  exact  stress 
constraints 
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are  non  linear  implicit  functions  of  the  design  variables 


TUeir  approximation  (3.5)  is  however  equal  to  the  exact  one 
at  any  point  along  the  D(a?}  scaling  line,  that  is 

o.,  - a..  on  D(a?)  (3.7) 

il  li  1 

and, in  part icular, at  the  intersection  of  the  line  with  the 
exact  constraint  surface 

o..«o.,»o.  t-l,n  (3.8) 
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Geometrical ly, the  exact  curved  constraint  surface  (3.6)  is 
approximated  by  a plane  perpendicular  to  the  axis  a^  and 
passing  by  the  intersection  of  the  exact  surface  with  the 
scaling  line  D(a?).  As  the  approximation  (3.5)  preserves  only 
the  value  of  the  function, and  not  of  its  derivatives , it  is  a 
zero-order  explicit  approximation  of  Che  stress  constraints. 


Turning  now  to  the  displacement  or  flexibility  constraints^ 
as  introduced  for  instance  by  BERKE  [ BIO  ] , their  explicit 

approximation  is  given  by  (2.18,  2.19),  that  is 
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where  q?,  q^^  are  the  nodal  displacements  of  the  i^^  element 
under  a real  or  virtual  loading  case. at  an  analysis  point  a?. 
In  the  exact  constraint 


u . 
J 


E 

i 


a . 
I 


< u. 
J 


(3.10) 


28. 


the  flexibility  coefficients  are  implicit  functions  of  the 
design  variables.  By  the  same  arguments  as  for  the  stress 
constraints,  the  approximate  value  of  the  flexibility  (3.9) 
is  equal  to  the  exact  value  at  any  point  along  the  scaling 
line  U(a?).  Consider  now  the  first  derivative  of  a given 
flexibility  (j  fixed).  Recalling  (2.16,  2.17)  we  have 
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are  independent  of  the  design  variables 
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and  hence 
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Using  the  expansion  of  K in  terms  of  the  finite  element  stiffness 
matrices  which  are  supposed  to  be  linear  in  the  design  varia- 
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At  the  analysis  point  a?,  using  the  definition  of  the  flexibi- 
lity coefficients  (3.9),  we  have 
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which  is  identical  with  the  value  derived  from  the  approximate 
explicit  form  (3.9) 
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We  conclude  therefore  that,  at  any  point  along  U(a?),the  appro- 
ximate values  of  the  flexibility  and  its  first  derivative  are 
equal  to  the  exact  ones 
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Geome tr ical ly , the  explicit  approximation  (3.9)  is  a surface 
which  is  tangent  to  the  exact  one  at  the  intersection  point  of 
D(a?)  with  the  exact  constraint  surface.  It  is  therefore  a first 
order  approximation  of  the  constraint. 


3.3.  Interpretation  as  approximations  of  the  problem  R 

A graphical  interpretation  is  presented  in  figure  3.1  for 
a two  dimensional  problem.  The  exact  stress  and  flexibility 
constraints  are  indicated  which  define  the  exact  problem 
(problem  R)  and  lead  to  the  optimum  design  in  R.  From  a design 
point  A°  at  which  the  structure  is  analysed,  the  optimality 
criteria  approach,  using  the  F.S.D.  concept,  consists  in  using 
the  explicit  approximations  of  the  constraints  examined  above. 
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They  are  represented  on  the  figure  by  a straight  line  perpen- 
dicular to  the  axis  for  the  stress  constraint , and  a curved 
line  tangent  to  the  exact  flexibility  constraint.  The  problem 
R is  clearly  replaced  by  a problem  A as  defined  in  (2.27)  and 
yields  a solution  which  is  indicated  by  the  point  P. 

The  example  is  easily  generalized  and  leads  to  the  conclu- 
sion that  the  classical  methods  using  optimality  criteria 
replace  the  problem  R by  a sequence  of  approximate  problems 
which  are  characterized  by  the  use  of  explicit  expressions  of 
the  constraints.  In  the  case  where  the  stress  constraints  are 
created  by  the  F.S.O.  concept,  the  approximation  corresponds 
Co  the  problem  A.  It  is  expected  that  a sequence  of  problem  A 
leads  to  the  solution  of  the  problem  R.  This  is  not  necessarily 
the  case  and  it  is  clear  that  it  does  depend  of  the  quality  of 
the  approximations  introduced  for  the  constraints. 

It  is  important  to  note  that  most  of  the  optimality  criteria 
methods  proposed  in  the  past^do  not  always  solve  exactly  the 
problem  A,  especially  when  multiple  inequality  flexibility 
constraints  are  imposed. 


3.4.  First  order  approximation  of  the  stress  constraints 

The  zero  order  approximations  of  the  stress  constraints 
introduced  by  the  F.S.D.  can  easily  be  replaced  by  first  order 
approximations.  A given  stress  component  in  an  element  k can  be 
written  as  a linear  combination  of  the  nodal  degrees  of  freedom 
of  the  element 
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where  t,  is  a line  of  Che  stress  matrix  of  Che  element, 
k 

n is  Che  number  of  load  cases.  This  expre ss ion , al ready  intro- 
c 

duced  in  (2.7),  is  formally  identical  to  that  of  the  flexibility 
(2.5)  and  therefore  the  same  treatment  can  be  applied  to  compute 


first  order  explicit  approximations 
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with  the  stress  constraint  coefficients 
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The  first  derivative  of  the  explicit  approximation  of  the 
stresses  are 
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Such  approximation  exhibits  obviously  the  same  property  as 
mentionned  for  the  flexibility  constraints 
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which  characterizes  first  order  approximations  of  the  cons- 
traints . 

If  the  concept  is  easily  introduced,  a certain  number  of 
difficulties  arise  in  practice.  The  most  important  one  is  that 
the  computation  of  the  stress  constraint  coefficients 
requires  to  take  into  account  a number  of  virtual  loading  cases 
which  is  equal  to  the  number  of  first  order  approximated  stress 
constraints, and  can  be  quite  large.  The  solution  to  that  diffi- 
culty will  come  from  a selection  of  those  stress  constraints 
that  require  a first  order  approximation  and  will  be  examined 
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later  in  the  context  of  hybrid  methods. 

A second  difficulty  comes  when  the  need  exists  to  use 
equivalent  stresses  to  express  the  constraints.  They  can  not 
be  expressed  as  linear  combinations  of  the  nodal  displacements 
as  in  (3.19).  Two  methods  are  proposed  in  that  case. 
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Consider  the  quadratic  form  of  the  displacements 

- q’^  H q (3.24) 

where  H is  a (n  x n)  positive  definite  matrix  of  constants. 

A virtual  loading  case  is  defined  by 

g'"  - 7 H q (3.25) 

and  its  virtual  work  is  given  by 
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At  a given  analysis  point  a?  in  the  design  space, it  yields 
the  exact  value  of  the  form  (3.24) 
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Just  as  for  the  linear  form  in  the  displacements, 
write  a first  order  explicit  approximation 
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where  qY  are  the  displacements  of  the  element  for  the 

virtual  loading  case  (3.25)  evaluated  in  a?. 


The  exact  first  derivative  of  X are  given  by 


but 


V 

q 


+ 


T 3K 


V 

q 


+ 


3X 

da. 

1 


jT  iai 

^ da. 


and 


dX 

da. 

1 


V 

g 


so  that 


(3.30) 


dX  T dK  V 


and , 
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which  shows  that  the  exact  first  derivatives  are  equal  to  the 
ones  derived  from  the  approximate  expression  (3.28). 

We  conclude  that 
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along  D(a?)  (3.33) 


that  means  that  (3.28)  is  a first  order  approximation  as 
announced . 

This  procedure  can  be  applied  when  a constraint  exists  on 
the  module  of  a displacement  or  in  the  case  of  an  equivalent 
stress.  Taking  the  example  of  the  von  MISES  equivalent  stress 
in  two  dimensions 
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we  introduce  the  matrix  V 
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so  that 
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The  stresses  in  an  element  are  given  by  the  stress  matrix 
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and  the  square  of  the  equivalent  stress  is 
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which  allows  to  compute  the  stress  constraint  coefficients  dY. 
This  implies  to  take  into  account  one  virtual  loading  case 
for  each  equivalent  stress  constraint  and  each  loading  case. 

Di£e£t_c£m£u^a^io^n_of^  £h£  £r£d^ent£  stresses 

It  is  also  possible  to  compute  the  first  order  approximations 
for  each  stress  component 


'^ikl 


(3.29) 


35. 


where  are  defined  by  analogy  with  (3.29). 

In  the  case  of  the  von  MISES  stress  (3.24)| it  turns  out  that 
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This  procedure  requires  three  virtual  loading  cases  for  each 
equivalent  stress  constraint. 


3.5.  A rigorous  first  order  optimality  criterion 

When  all  the  constraints  defined  in  the  original  problem 
R (2.1  to  2.4)  are  replaced  by  their  first  order  explicit 
approximations, a first  order  approximate  problem  is  obtained 
which  is  denoted  problem  A1  and  is  written 
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A graphical  representation  of  the  probles'Al  is  given  in 
figure  3.2yWhich  is  to  be  compared  with  the  figure  3.1 
illustrating  the  problem  A. 

The  first  order  optimality  conditions  for  the  problem  R are 
given  by  the  classical  KUHN-TUCKER  conditions 
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The  variables  a^^  are  the  primal  variables,  while  r,  s,  t,  w 
are  the  dual  variables  which  have  the  meaning  of  lagrange 
multipliers.  They  are  positive  for  an  active  constraint  and 
zero  otherwise.  They  are  not  allowed  to  be  negative  like  in 
the  problems  with  equality  constraints.  Geometrically,  these 
conditions  means  that  the  gradient  of  the  objective  function 
is  orthogonal  to  the  plane  tangent  to  the  active  constraints 
at  the  optimum.  The  KUHN-TUCKER  conditions  are  only  necessary. 
They  become  also  sufficient  in  a convex  problem  which  is  not 
the  case  in  general  for  the  problem  R. 

The  KUHN-TUCKER  conditions  can  be  written  for  the  first 
order  approximate  problem  A1  (3.32)  by  using  the  explicit 
expressions  of  the  constraint  gradients  (3.17  and  3.22), 


It  yields  a rigorous  firsc  order  optimality  criterion. 
Its  expression  is  clarified  by  defining  three  classes  of 
variables 
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This  classification  is  analogous  to  the  classical  separation 
in  active  and  passive  members.  With  these  notations,  we  obtain 
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The  essential  differences  between  this  expression  of  the 
optimality  criteria  and  the  classical  ones  are  the  following. 

First, it  is  established  for  the  first  order  approximate 
problem  A1  formulated  in  terms  of  inequality  constraints. 

The  selection  between  active  and  inactive  constraints  is  therefore 
built  in  the  formulation  and  can  be  treated  rigorously, which  is 
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not  always  the  case  in  classical  formulations. 

Second, it  uses  a consistent  first  order  approximation  for 
all  the  constraints  while  the  classical  F.S.D.  concept  leads 
to  consider  the  stress  constraints  as  side  constraints. 

This  step  is  essential  for  convergence  considerations.  At  the 
optimum  the  necessary  conditions  (3.33)  are  in  general  not 
satisfied  using  the  F.S.D.  which  is  known  to  lead  sometimes 
to  non  optimum  solutions.  At  the  opposite, the  expression  (3.37) 
is  deduced  from  these  necessary  conditions  and  the  problem  is 
reduced  to  the  approximation  of  the  problem  R by  a sequence  of 
problems  Al,  which  does  not  necessarily  yield  a converging  process. 

It  is  to  be  mentionned  that  the  expression  (3.37)  corresponds 
exactly  to  previously  derived  optimality  criteria  [ BIO,  G13, 

Bll,  T3  ] when  only  one  flexibility  constraint  is  imposed  and 
no  stress  constraints.  In  the  case  of  multiple  constraints  the 
difficulties  of  solving  exactly  the  system  of  non  linear  inequa- 
lities (3.37)  joined  to  the  computational  handicap  of  using  first 
order  approximations  for  the  stress  constraints, led  many  authors 
to  using  various  additional  approximations.  These  are  responsible 
for  the  poor  convergence  properties,  or  even  divergence,  encoun- 
tered in  certain  cases. 

The  difficulty  of  solving  (3.37)  is  in  fact  essentially 
related  to  the  fact  that  the  problem  A1  is  presented  in  its 
primal  form.  Efficient,  rigorous  and  systematic  solutions  can 
however  be  derived  by  using  the  dual  formulation  of  the  same 
problem  as  shown  in  the  next  chapter.  Therefore  the  presentation 
in  the  form  (3.37)  is  only  intended  to  serve  as  a definition  of 
the  optimality  criterion  associated  to  the  problem  A1 . 


3.6.  Hybrid  optimality  criteria 

The  computational  handicap  of  using  first  order  approxima- 
tions for  the  stress  constraints  is  due  to  need  to  introduce  an 
eventually  large  number  of  virtual  loading  cases  in  the  analysis. 


which  is  avoided  when  Che  zero  order  F.S.U.  approximation  is 
used . 


In  practical  structures  it  is  observed  that  most  of  the 
stress  constraints  are  approximated  with  sufficient  accuracy 
by  Che  F.S.D.  while  some  other,  in  limited  number  most  of 
time,  reveal  critical  and  require  a first  order  approximation. 
It  is  not  possible  in  general  to  decide  a priori  which  cons- 
traint is  going  to  be  critical  and  Che  number  of  such  critical 
constraints  can  vary  at  each  iteration,  that  is  after  each 
reanalysis . 

An  automatic  selection  of  these  critical  constraints  can 
be  based  on  Che  preceding  interpretations  and  presented  as 
f ol lows . 


(3.38) 


Che  constraint  is  potentially  critical.  The  first  condition  is 
rather  evident  and  means  that  the  stress  in  the  element  is 
equal  or  close  Co  Che  prescribed  limit  which  implies  that  the 
constraint  is  active  or  potentially  active.  The  second  condition 
arises  from  Che  fact  that  a first  order  approximation  reduces 
to  a zero  order  approximation  when 
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Che  term  is  dominant, Che  situation  is  not  far  from  Chat 

existing  in  the  isosCatic  case  where 

In  ocher  words, if  in  Che  sum  (3.40)  Che  contribution  of  the 
other  elements  (i  ^ k)  is  negl igib le , the  F.S.D.  zero  order 
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approximation  is  a good  approximation. 

After  each  analysis  step^the  conditions  (3.38)  can  be  used 
to  define  which  stress  constraints  are  to  be  first  order  appro- 
ximated while  the  others  are  treated  by  F.S.D.,that  means  are 
treated  as  side  constraints  by  using  the  definition  (2.15)  for 
the  minimum  size.  This  procedure  leads  to  a sequence  of  problems 
A1  which  involve  each  a different  number  of  constraints.  The 
corresponding  optimality  conditions  can  be  derived  just  as  in 
(3.37)  . 
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. DUAL  FORMPLATIOM  AMD  OPTIMALITY  CRITSKIA 


4.1.  The  primal  problem 


The  classical  optimAliCy  criteria  have  been  derived  by 
using  a lagrangian  formulation  to  solve  a simplified  version 
of  the  problem  A (2.27)  where  a selection  between  active  and 
inactive  constraints  allows  to  define  the  problem  with  equality 
constraints.  The  generalized  versions  presented  in  the  prece- 
ding chapter  are  based  on  the  problem  Al  (3.32)  and  treat  the 
case  with  inequality  constraints.  Eventually  a hybrid  situation 
is  considered  where  only  a limited  number  of  stress  constraints 
are  first  order  approximated.  To  take  these  various  situations 
into  account  we  write  the  basic  problem  in  the  form 


problem 


i*l,n  (4.1) 

j - 1,  m (4.2) 

i*l,n  (4.3) 

i - 1,  n (4.4) 


where  we  consider  that  the  m constraints  (4.2)  group  all  the 
first  order  approximated  constraints,  that  is,  the  true  flexibility 
constraints  and  the  first  order  approximated  stress  constraints 
(3.20).  The  minimum  variable  sizes  are  defined  as  the  a* 

(2.15)^  that  is, result  from  the  choice  of  the  maximum  between  the 
true  minimum  sizes  or  the  values  computed  by  the  stress  ratio 
algorithm  of  F.S.D. 

In  such  a way, all  the  above  mentionned  situations  are  covered. 

The  problem  A'  is  called  the  primal  problem  because  it  is  ex- 
pressed in  the  primal  variables  a^.  It  is  not  necessarily  convex 
as  the  c^j  coefficients  in  (4.2)  are  not  always  all  positive. 
However  if  e change  of  variables 
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X . ■ — 

1 a . 
1 


(4.5) 


is  introduced  the  problem  becmmes  strictly  convex  and 
separable, which  allows  to  derive  easily  a dual  formulation. 


4,2.  The  dual  problem 

The  fact  that  the  classical  optimality  criteria  are  based 
on  a lagrangian  method  suggests  to  use  a dual  method  to  solve 
the  problem  A*.  The  change  of  variables  (4. 5), which  is  necessary 
to  have  a strictly  convex  problem, is  not  introduced  explicitely 
as  they  disappear , and  the  solution  will  be  expressed  in  terms  of 
the  dual  variables  only.  These  are  the  lagrange  multipliers  asso- 
ciated to  the  various  constraints. 


4,2.1.  The  complete  dual  formulation 

In  the  sense  of  WOLFE  [ W2  ] the  dual  problem  associate 
to  the  problem  A*  is 


max  £ (a,  r,  w,  t) 


- E a^  - I r.  (u  - u.)  - E w^  (a^  - a^^) 

i j ^ i 


- E t.  (a. 
i ^ ^ 


(4.6) 


with  the  constraints 


TIT  " ° 


(4.7) 


r.  > 0 
J 


w . > 0 
1 


tj  > 0 


(4.8) 


The  equality  constraints  (4.7)  insuring  the  stationarity  of 
the  lagrangian  function  £ yields,  as  the  problem  is  separable. 
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th«  explicit  relations 


£ r . c . . 

i ^ 


2 

i p . 1 . . -t . 


(4.9) 


I 1 


allowing  to  coapute  the  priaal  variables  a^  when  the  dual 
variables  r j , w^t  t^  are  known.  If  the  a^  given  by  (4.9)  are 
substituted  into  the  functional  (4.6), the  problem  reduces  to 
a maximization  problem  in  the  (m  * 2n)  dual  variables  r j , w^,  t^ 
subjected  to  the  (m  * 2n)  non  negativity  constraints. 

Formulated  in  this  way  the  dual  approach  is  not  advanta- 
geous as  it  involves  more  variables  than  in  the  primal  problem. 
However  the  simple  form  of  the  primal  constraints  (4.3)  and 
(4.4)  allows  to  use  a formulation  reduced  in  terms  of  the  only 

variables  r . . 

J 


4.2.2.  The  reduced  dual  formulation 

The  foundation  of  the  method  used  in  this  section  is  due 
to  FALK  [ F12  ] . Consider  the  primal  problem 


min  U * £ p.  t.  a.  i ■ 1,  n (4.10) 

i i 1 1 

_ c , . 

with  u.  ■ £ — ^ >0  j • 1,  m (4.11) 

J i *i 

•i  ^ 

where  C is  the  convex  and  compact  set  of  the  variables  dsfinmd 
between  the  minimum  and  maximum  values  by 

C - ( a ; sj  < a^  < I.  ) (4.12) 

The  definition  is  necessary  to  exclude  the  possibility  of 
infinite  solutions  in  the  minimization  process. 
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h 


The  dual  formulation  of  FALK  is  simply 
max  Y(r) 

j - 1,  m 

where  Y(r)  ie  en  auxiliary  function  given  by 
Y<r)  ~ 


(4.13) 


with  r . > 0 
J 


min  [ I p.  i.  a.  - £ r.  (u. 

• 1 1 1 1 J 

a . 1 J 

1 ■' 


I 1 

i ‘i 


(4.14) 


for 


a . < a . < a . 
—Ill 


or,  more  simply,  as  the  problem  is  separable, 

Y(r)  - £ p.  1.  a.(r)  - £ r.  (u,  - u.  [ a. (r)  ] ) 

jJj  J *■ 

(4.15) 

The  solution  of  (4.14)  leads  to  introduce  the  quantities 


P • (r)  ■ — r—  ^ c , , r . 
1 Pi^i  j 13  3 


j - 1,  m 
i ■ 1 , n 


(4.16) 


The  primal  variables  a.(r)  are  then  expressed  in  terms  of  the 


variables 

using  the 

quantities  6 

a . ■ 

/r. 

if 

a?  < B. 

1 

1 

—1  1 

a . ■ 

a . 

if 

B. 

1 

—1 

1 

a . ■ 

1 

a . 

1 

if 

•i  * «i 

-2 

a . 
1 

2 
a . 
—1 


(4.17) 


The  explicit  expression  of  the  dual  problem  still  requires 
the  definition  of  the  flexibilities 
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“j  ' *i‘'>  ' * ! ^ 


j - 1,  m 


(4.18) 


1 I 


The  gradient  of  lagrangian  function  can  be  obtained  as  follows 


3 n 3a.  m 3u 

tJt  - f ^ tf:  " ? " “j  - “j 

J 1 J k J 


(4.19) 


but,  using  (4.16)  and  (4.17) 


3a . , c . . 

i:  - i LI_ 

3r  . 2 p . i . a . 

J ^111 


1 ■ 1 , n 


(4.20) 


3a. 
1 

3r . 
J 


- 0 


i - n+1,  n (4.21) 


. n. 

if  n denotes  the  number  of  free  variables,  that  is  belonging  to 
the  interval  3 C • From  (4.20)  and  (4.21)  it  follows  that 


n 3a . - n c . . 

Z p.  ^ - 1 Z _iJ. 

1 I 3r . 2 . a . 

1 J 1 1 


(4.22) 


and  from  (4.18) 


3u, 


n c . , 3a  . 
k _ ik  1 

3r.  “ ■ . “2“  3r. 

J 1 *£  J 


(4.23) 


and 


m 
E r 
k 


k 3r . 
J 


n. 

, n c . 
- i I -i. 
2 . a . 

1 1 


( i.24) 


Finally  the  gradient  (4.19)  reduces,  using  (4.22)  and  (4.24)  to 


9r. 


U . “ u . 

J J 


(4.25) 
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All  the  ingredients  necessary  to  the  solution  of  the  dual 
problem  (4.10)  (4.11)  have  been  evaluated. 

The  reduced  dual  problem  involves  the  maximization  of  a 
lagrangian  function  (4.10)  under  the  constraints  (4.11)  in  a 
space  of  dimension  m equal  to  the  number  of  first  order  appro- 
ximated constraints.  A certain  number  of  methods  can  be  applied 
to  this  maximization  problem  and  some  of  them  are  described  in 
Che  next  section. 

It  must  be  emphasized  Chat  when  the  number  of  dual  variables 
is  much  smaller  than  the  number  of  design  variables,  which  is 
often  Che  case,  a significant  advantage  exists  for  the  dual  for- 
mulations. An  important  property  is  their  ability  Co  solve  exactly 
a problem  A*  associated  to  a given  rigorous  or  hybrid  optimality 
criterion, without  introducing  any  approximation  in  the  selection 
of  active  and  passive  constraints, nor  in  the  treatment  of  the 
side  constraints. 

Domains)  £f_d£f^n^t^on  £n£  £la.nes  of  d isc£nc inu^t^ 

The  relations  (4.17)  indicate  the  existence, in  the  dual  space, 
of  planes  of  discontinuity  for  the  second  derivatives  of  the 
dual  objective  function.  They  are  given  by 


Z c . . r.  - p,  I . a. 

j ij  j '^11-1 


(4.26) 


Z c..  r.  • p.  t . a. 

j ij  J ^111 


(4.27) 


When  crossing  these  planes, Che  definition  of  the  lagrangian 
function  (4.15)  is  modified  according  Co  (4.17).  In  each  sub- 
regions  of  the  dual  space,  called  domains  definition,  Che  dual 
objective  function  is  modified  due  to  different  distributions 
of  the  primal  variables  in  free  and  fixed  variables. 


1 
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4.3.  The  dual  algorithms 

Various  maximisation  methods  can  be  applied  to  solve  the 
reduced  dual  problem  (4.10)  (4.11)  corresponding  to  the  problem 
A*  (4.1  to  4.4).  Two  types  of  algorithms  have  been  developped 
and  used  in  the  numerical  applications.  The  first  one  is  a 
first  order  algorithm  which  requires  only  the  first  derivatives 
of  the  objective  function,  the  second  one  is  a second  order 
algorithm  requiring  the  second  derivatives  of  the  dual  objective 
function. 


4.3.1.  A first  order  dual  algorithm 

The  simple  algorithm  proposed  here  is  based  on  the  projec- 
ted gradient  method  described  by  ROSEN  in  [ R1  1 .It  takes  a 
very  simple  form  due  to  the  fact  that  the  constraints  reduce  to 
the  non  negativity  of  the  dual  variables. 

Initialization 

(i)  define  an  admissible  starting  point  such  that 


r?  > O 
J 


j - 1,  m 


^t£r£t^v^  £r£C£dur£ 


(ii)  Let  r^  the  point  observed  at  the  iteration  X 

Compute  - the  primal  variables  by  (4.16)  and  (4.17) 

- the  flexibilities  Uj  by  (4.18) 

The  search  direction  is  given  by  (4.25)  that  is 


z 


X 

j 


- 0 


if 


r^  ■ 0 and  u^  < u. 
J J J 


z 


X 

j 


if 


X __ 

r . > 0 or  u . > il. 
J J J 


i 
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(iii) 


if  I I z II  > e (given  tolerance)  go  Co  (iv) 


otherwise  r is  the  global  minimum  (of  Che  problem  A') 


determine  the  optimal  step  length  along  by  the 
linear  search 


Y(r^  ♦ z^)  - max  [Y(r^  + t z^)  ] 


where  3^  - ( t > 0 ; (r^  + t z^)  > 0 


During  this  linear  search,  a dual  variable  may 
eventually  reach  the  assigned  limit  r^  ■ 0 


Progress  to  the  step  X+1  by  computing 


,X+1  X ^ X X 
r • r + T z 


go  to  (ii) 

This  procedure  is  very  simple.  It  can  be  improved  significantly 
by  using  an  algorithm  of  conjugate  direction, which  is  not  pre- 
sented here  for  simplicity, as  it  is  rather  classical  [ FA  ] . 

In  the  application  described  later  in  this  work,  such  a conjugate 
direction  algorithm  is  used. 


A. 3. 2.  A second  order  dual  algorithm 

The  algorithm  developed  here  is  based  on  Che  generalized 
NEWTON  method  particularized  for  the  case  of  linear  constraints. 
A preliminary  problem  in  such  a method  is  to  determine  the  direc- 
tion of  Che  search.  Consider  the  case  where  the  m first  dual 
variables  are  positive  (active  primal  constraints)  and  the  m>m-m 
next  ones  equal  to  zero  (inactive  primal  constraints) 


j 


r . 
J 


0 


J 
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- m+l. 


m 


(4.28) 


The  gradient  g and  the  hessian  matrix  F are  partitionned 
accordingly 


g - Vy 


r * 1 

g 

2 

F 

F 

F - 7 y - 

_ 

g 

- 

L J 

(4.29) 


If  N denotes  the  matrix  defining  the  base  of  the  active  cons- 
traints so  that  (4.28)  is  written 


N*  r - 0 


(4.30) 


the  matrix  N has  the  simple  structure 


(4.31) 


where  I is  a (m  x m)  unit  matrix.  The  NEWTON  search  direction  is 
given  by  [ G3,  G5 , P9  ] 


with 


(4.32) 


P - I - N (N^  f'^  N)”^  N^ 
Introducing  the  notation 

- (N^  f"^  N)"^  N^  F~^ 


V 


g 


(4.33) 
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while  (4.33)  is  simply  written 


T 1 

— * *'<  — » 
-g-F  z-g-F 


i-1 

F g 


(4.38) 


As  expected  in  a second  order  method^ the  determination  of  the 
search  direction  z requires  the  inversion  of  the  hessian  matrix 
of  the  lagrangian  considered  as  a function  of  the  m free  dual 
variables  only. 

To  determine  the  active  constraints,  a very  simple  selection 
rule  can  be  established  if  the  sign  of  the  multipliers  (4.38)  is 
examined  only  at  a stationary  point  of  the  objective  function 
where,  by  definition,  the  gradient  g is  zero, which  implies 


g 


(4.39) 


which  does  not  require  the  evaluation  of  the  matrix  F*  present 
in  (4.38).  The  inversion  of  F is  therefore  only  achieved  when 
the  basis  of  active  constraints  is  known, which  avoids  the  numerical 
difficulties  that  could  appear  in  case  of  linearly  dependent 
flexibility  constraints. 

Algorithm 

Initialization 


(i) 


select  an  admissible  starting  point 


“VO  ^ 

r.  > 0 
J 

r . - 0 
J 


j - 1,  m 


j - m+ 1 , m 


such  that  the  hessian  matrix  f”  is  non  singular 
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(ii) 


(iii) 


(iv) 


£r£ce^dure^ 

Let  be  the  point  obtained  at  iteration  X 

. Compute  the  primal  variables  a^(r^)  by  (4.16)  and 
(4.17)  and  the  flexibilities  uj  by  (4.18) 

. Determine 

X — . , ^ 

J“l.® 

• If  1 I 8^  I 1 < t (given  tolerance)  go  to  (v) 
otherwise  compute 

Compute  the  optimal  step  length  along  by  the  linear 
search  in  the  space  of  the  m dual  variables 

Y(r^  + z^)  - max  [ y + t z^)  ] 

t€S^ 

where 

- { T > 0 : (r^  + T z^)  > 0 } 


I 

Compute 


~X  + 1 


•vX 

r 


^ X X 

♦ T Z 


1-1 


go  to  (ii)  with  eventually  m 
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js! 


:! 


f! 


[ 


(v) 


Compute 


-A  A - 
g . • u . - u . 
J J J 


j • m+1,  m 


If  max  t gj  ) < 0 


'v»A 

then  r is  the  global 
maximum  (of  the  problem  A') 


otherwise  go  to  (ii)  with  m ■ m+1  by  incorporating 
in  r the  dual  variable  which  has  the  largest 
positive  component  gj  . 

To  fully  determine  the  procedure^it  is  still  necessary  to  give 
the  explicit  form  of  the  Hessian  matrix  of  the  objective  function. 
From  (4.19)  to  (4.25)  we  obtain 


jIl 


3r . 3r, 
J k 


1 c . , c . , 

- i I 

2 • .3 

1 p . £ . a . 

Ill 


(4.40) 


where 


• -I 

1 ■ 1 , n 


J - 1.  m 


k - I, 


and  n is  the  number  of  free  primal  variables  as  defined  by  (4.17). 


4.3.3.  Linear  search  algorithm 

The  difficulty  of  the  linear  search  along  a given  direction 
arises  from  the  fact  that  the  dual  space  is  subdivided  into  various 
definition  domains  where  the  objective  function  takes  different 
forms.  Across  a plane  of  discontinuity  the  number  of  free  primal 
variables  n is  modified  and  the  second  derivatives  (4.40)  jump  to 
other  values. 

The  linear  search  can  be  represented  indifferently  by  the  maximi- 
zation of  the  lagrangian 

m _ 

max  y(t)  - E p.  £.  a.(T)  ♦ E r.(T)  [ u.(t)  - u.  J 

r i I 1 1 j J J 1 


1 


■j 


i 

y 

I 

i 

I 

i 

] 


1 


t 

♦ 

i 

j 
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or  by  the  problem  of  finding  a vanishing  first  derivative 


m 


y'(t)  - z V y(t)  - Z z.  [ u.(t)  - u.  1 - 0 (4.41) 

j J J J 


This  last  formulation  presents  some  advantages  and  can  be  used 
as  follows.  Introduce  the  constants 


m 


X - E z . u . 

j ^ ^ 


Y,  - 


E c.  . r. 


I Pi^i  j IJ  J 


(4.42) 


m 


Z.  - — E c,  , z, 

I p . 1.  . 11  1 

1 1 J •*  •* 


where  denotes  the  current  value  of  the  dual  variables. 
The  problem  consists  in  finding  the  zero  of  the  function 


n Z.p . i. 
y'(t)  - r ^ ^ ^ 


7 aTTO 
1 1 


- X 


(4.43) 


where  the  design  variables  aj^(T)  are  given  by 


a.(T)  ■ Y.  + T Z. 
1 11 


if 


2 —2 
af  < Y.  -f  T Z.  < aT 
—I  1 11 


(that  is  if  i - 1,  n)  (4.44) 


a . ■ a . 
1 —1 


or 


a . * a . 
1 1 


otherwise  . 


This  problem  can  be  solved  by  the  NEWTON-RAPHSON  or  by  the 
TCHEBYCHEV  method  [ B13  ] . The  higher  order  derivatives  required 

are  given  by 


y"(t)  - 


. I . i"l‘i 

2 , 3 

1 a . 


1 ■ 1 , n 


(4.45) 
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i 


I 

1 


Y (t) 


The  function  y'C'c)  is  continuous  across  the  planes  of  disconti- 
nuity^but  the  higher  derivatives  are  only  piecewise  continuous. 
They  exhibit  a discontinuity  each  time  a change  in  the  number  of 
free  variables  occurs.  This  is  the  case  for  the  following  values 
of  the  step  : 


(4.46) 


It  is  therefore  required  to  take  some  measure  to  avoid  a divert 
gence  of  the  NEWTON-RAPHSON  procedure.  A simple  technique  con- 
sists in  finding  by  (4.46)  the  values  of  the  step  where  a dis- 
continuity occurs, and  select  the  interval  in  which  the  maximum 
is  obtained. 

4.4.  Relation  with  the  classical  optimality  criteria 

The  optimality  criterion  for  a single  flexibility  constraint 
as  proposed  by  BERKS  [ BIO  ] appears  to  be  equivalent  to  a dual 
formulation  of  a problem  Al,  in  the  sense  of  FALK  [ F12  ] . 

In  this  simple  case  the  maximization  problem  (4.13)  is  solved 
analytically  as  the  auxiliary  function  y(<')  depends  only  of  one 
variable . 

Turning  to  the  optimality  criteria  of  TAIC  and  KERR  [ T3  ] , a 
change  of  variables  in  the  dual  problem  (4.13) 


X.  - 7-  (4.47) 
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and  the  application  of  the  second  order  NEWTON  method  to  the 
new  expression  of  the  problem. yields  the  iterative  scheme 


v+1 


X ♦t(X)  (u-u) 


(A. 48) 


where  v is  the  iteration  number  and  X the  Hessian  matrix  of 
the  new  dual  objective  function,  deduced  from  (4.40) 


Jk 


A. 


ax.  3x. 

j k 


_i 

2X 


c . . c . , 
1 j ik 


(4.49) 


• J. 

1 p . i . a . 
k ^111 


1.  n 


1 ''' 
1 • m 


k - 1. 


In  the  process  (4.48), the  step  length  has  to  be  computed  to 


V * X **  V 

maximize  the  lagrangian  along  the  direction  (X  ) (u  - u ) 


while  remaining  in  the  admissible  dual  domain  Xj  > 0. 

With  the  present  notation,  the  original  algorithm  of  TAIG  and 


KERR  replacea  (4.48)  and  (4.49)  by 


.v+1  .V  ^ ,„v.-l  V. 

X "X  ♦(X)  (u-u) 


and 


(4.50) 


X-i.  " 
Jk 


1 

2x: 


*^ij*^ik 
p . i . a . 

Ill 


i ■ 1 , n 


j ■ 1,  m 


k > 1 , m 


This  algorithm  consists  in  maximizing  the  lagrangian  using  the 
simple  NEWTON  method  without  computing  the  optimal  step  length 
which  is  equivalent  to  ■ 1 in  (4.48).  In  the  TAIG  and  KERR 
method  the  reference  strain  energy  densities  are  interpreted 
as  the  inverses  of  the  dual  variables  associated  to  the  flexi- 
bility constraints.  The  reasons  of  the  instability  of  this  more 
intuitive  algorithm  becomes  clear.  First  the  method  can  not  avoid 


the  Xj  Co  become  negetive,  chat  is,  ignores  the  non  negativity 
constrainCa  of  the  dual  variablea. 

Second, the  fact  Chat  the  optimal  atep  length  ia  not  computed 
can  lead  to  a divergence  when  the  atarting  point  ia  choaen  too 
far  from  Che  optimum.  It  ia  only  cloae  to  the  optimum  Chat  the 
atep  length  becomea  cloae  to  unity  in  the  generalized  NEWTON 
method . 

Third, they  ignore  the  fact  that  Che  dual  apace  ia  aubdivided 
into  varioua  definition  domain.  The  conaequence  ia  that  the 
choice  of  free  and  fixed  deaign  variablea  can  not  alvaya  be 
correct . 


The  dual  mechoda  propoaed  in  thia  chapter  appeara  aa  a 
aound  mathematical  way  to  aolve  the  problema  correaponding  to 
the  rigoroua  (or  hybrid)  optimality  criteria. 

Their  computational  performancea  are  excellent  and  will  be 
demonatraCed  in  the  chapter  dealing  with  applicationa . 

They  preaent  Che  fundamental  advanCagea  of  allowing  any  combi- 
nation of  zero  or  firat  order  approximated  conatrainta, and  of 
a correct  aelection  between  active  and  paaaive  conatrainta. 

The  determination  of  Che  valuea  of  the  lagrange  multipliera 
remaina  poaaible  even  when  the  conatrainta  are  linearly  dependent. 
The  convergence  propectiea  are  dependent  of  the  validity  of  the 
approximation  of  Che  'problem  R by  a aeqnence  of  problema  A'. 

Thia  impliea  Chat  if  the  convergence  to  the  solution  of  the 
problem  A*  ia  alwaya  guaranteed,  the  convergence  of  the  aequence 
of  problema  A*  to  the  aolution  of  Che  problem  R can  not. 

Thia  remaina  an  eaaencial  difficulty  uaing  optimality  criteria 
which  will  be  given  additional  attention  after  having  introduced 
the  mixed  methoda  in  the  next  chapter. 


4.3.  Two  aimple  illuatraciona 


A two  dimenaional  problem  ia  uaed  to  illuatrate  graphically 
Che  concepCa  introduced  in  the  dual  methoda. 

The  two  bar  truaa  achemaCized  on  figure  4.1  can  be  analytically 
formulated  aa  followa  : 
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min  U(a^) 


with  u 


V 


p 1 /2  (a^  ♦ aj) 

Pl/2  .1  + 1 N c 3 

2E  aj'  2 2E 

Pt/2  ,l  1.^1  Pt/T 

2E  'aj^  ■ a^^  2 2E 

1 < a^  < 2 
1 < a2  < 2 


(4.51) 

(4.52) 

(4.53) 

(4.54) 

(4.55) 


The  design  space  and  the  constraints  are  also  illustrated  in  the 

A 

figure  4.1.  The  optimum  corresponds  to  the  joint  * *3  • ^2  ** 
at  which  only  one  constraint  is  active  (horizontal,  u,  displace- 
ment) . The  formulation  of  the  dual  problem  involves  two  lagrange 
multipliers  r^  and  r^  associate  to  the  constraints  (4.52)  and 
(4.53).  The  planes  of  discontinuity  defined  in  (4.26)  (4.27)  are 
illustrated  in  figure  4.2  and  correspond  to 


ri  + r2  - 1 


- r. 


- r. 


- 4 


(4.56) 
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It  yields  6 subdomains  in  the  dual  space.  In  each  of  them,  the 
expression  of  the  objective  function  in  terms  of  the  dual  variable 
is  different  and  given  in  explicit  form  in  table  4.3. 

The  optimum  in  the  dual  space  occurs  in  the  subdomain  II  for 
r^^  • and  r2  * 0 and  a value  of  the  dual  lagrangian  objective 
function  yfr)  * j.  The  linear  search  along  a direction  r2  * •• 

constant  is  illustrated  on  figure  4.4  as  well  as  the  continuity 
properties  of  the  lagrangian  function  y(‘()  <>nd  its  derivatives 
defined  in  (4.41). 


Another  simple  example  is  provided  by  the  10  bar  truss  defi- 
ned in  figure  4.5.  In  this  problem  no  flexibility  constraints  are 
imposed  but  the  stress  constraints  are  < 25000  psi  in  all  the 
bars  but  in  bar  8 when  the  limit  Og  is  given  different  values. 

Using  the  classical  F.S.D.  approach,  as  reported  by  BERKE 
in  a similar  problem  [ B8  ] , the  solution  converges  rapidly  to 

the  optimum  for  Og  < 36800  psi  while  for  Og  > 45000  psi  a fast 
convergence  is  observed  but  to  a non  optimum  F.S.D.  design. 

For  36800  ^ ^ 45000  the  convergence  is  very  slow  and  does  not 

lead  to  the  optimum  solution  after  30  reanalyses  nor  even  to  a 
design  satisfying  the  F.S.D.  concept.  An  exception  exists  for 
Og  > 37500  psi  where  the  convergence  to  a non  optimal  F.S.D. 
design  is  obtained  in  11  reanalyses.  This  behavior  is  illustrated 
in  figure  4.6.  The  same  problem  has  been  solved  using  the  genera- 
lized optimality  criterion  based  on  first  order  approximation  of 
the  stress  constraints  and  using  the  dual  algorithm  described 
above.  The  results  are  given  in  figure  4.7.  Besides  the  fact  that 
the  optimum  is  obtained  for  > 36800  psi,  the  number  of  reanalyses 
remains  approximately  the  same  for  all  the  values  of  and  cor- 

O 

responds  to  a very  fast  convergence.  The  designs  corresponding  to 
Og  ■ 25000,  Og  ■ 30000,  Og  > 37500  are  presented  in  table  4.8. 


variable*  flexibilities  dual  objective  function 


TABLE 
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10  BAR  TRUSS  WITH  STRESS  CONSTRAINTS 

i 

t 

i OPTIMUM  SOLUTIONS 

.1 

-1 


‘’8  " 

25000 

Og  = 30000 

Og  > 37500 

bar 

section 

stress 

section 

stress 

section 

stress 

1 

8.0621 

25000 

8.0138 

25000 

7.9414 

25000 

2 

3.9379 

25000 

3.9462 

25000 

3.9586 

25000 

3 

0.1000 

- 15533 

0. 1000 

- 13462 

0. 1000 

- 10355 

A 

0.  1000 

- 15533 

0.1000 

- 13462 

0.1000 

- 10355 

5 

7.9379 

- 25000 

7.9862 

- 25000 

8.0586 

- 25000 

6 

0. 1000 

0 

0.1000 

- 10000 

0.1000 

- 25000 

7 

5.5690 

25000 

5.6373 

25000 

5.7397 

25000 

8 

5.  7AA7 

- 25000 

4. 7304 

- 30000 

3.7160 

- 37500 

9 

0.1000 

21976 

0.1000 

19038 

0.1000 

14645 

10 

5.5690 

- 25000 

5.5807 

- 25000 

5.5983 

- 25000 

Weight 

1593.18  (FSD) 

1545.91  (FSD) 

1500.82  (non  FSD) 
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5.  RELATIONS  BETWEEN  OPTIMALITY  CRITERIA  AND  MATHEMATICAL 
PROGRAMMING  APPROACHES 

5.1.  Interpretation  of  the  generalized  optiiaa  lity  criteria 

In  chapter  3,  it  has  been  established  that  consistent  first 
order  generalized  optimality  criteria  (G.O.C.)  could  be  defined 
by  the  application  of  the  KUHN-TUCKER  necessary  conditions  to 
an  approximate  problem,  denoted  problem  A1  in  (3.32),  obtained 
from  the  exact  problem  R,by  replacing  the  exact  implicit  cons- 
traints by  their  first  order  explicit  approximations. 

Mixing  aero  and  first  order  approximations  of  the  constraints 
leads  to  what  has  been  named  a hybrid  optimality  criterion  (H.O.C.) 
corresponding  to  a problem  A'  defined  in  (4.1).  The  classical 
optimality  criteria  (C.O.C.)  results  from  the  introduction  of 
additional  approximations  in  the  treatment  of  the  constraints, 
in  the  selection  of  active  and  passive  constraints  and  in  the 
minimization  algorithms  themselves.  These  additional  approxima- 
tions were  motivated  by  the  desire  to  obtain  simple  redesign  al- 
gorithms . 

In  chapter  4 it  has  been  shown  that  the  application  of  the 
dual  methods  furnishes  an  alternative  way  to  obtain  the  solution 
of  the  problems  A1  or  A'  which  is  rigorous  and  efficient. 


The  solution  of  the  problems  A1  or  A*  by  primal  methods 
of  mathematical  programming  does  appear  to  be  rather  inefficient 
due  to  the  non  linearity  of  the  constraints  even  in  their  first 
order  approximated  form.  However,  as  suggested  by  various  authors 
( R4,  Jl,  S8,  FIS  ] , the  change  of  varihbles 


makes  the  constraints  very  shallow.  Indeed  they  are  linear  in 
the  case  of  isostatic  structures.  The  fact  that  the  objective 
function  W(x^)  becomes  non  linear  is  not  a handicap  since  it  is 
known  in  an  explicit  form.  In  the  new  design  space  x^,  which  is 


called  the  inverae  deaign  space,  the  approximate  flexibility 
and  stress  constraints  (3.32)  are  written 


u.  • Z C. . X.  <u. 
J { 1 J 


j ■ 1. 

i ■ 1 , n 


(5.2) 


'kl 


I “ikt  ‘i 


k ■ 1 , n 
I - 1,  n^ 
i ■ 1 , n 


(5.3) 


The  problem  A1  is  rewritten  in  the  inverse  design  space  and 
denoted  problem  L 


min  W(x.)  ■ Z 


problem 


o 


with 


!iii 

*i 

(5.4) 

'ij  “l  * “j 

(5.5) 

“ikt  *1  * ”k 

(5.6) 

-i  » ii 

(5.7) 

M 

A 

Ml 

(5.8) 

When  the  first  order  approximations  (5.6)  of  the  stress  cons- 
traints are  totally  or  partially  replaced  by  their  zero  order 
% 

approximations  x^,  derived  from  F.S.D.,  we  denote 


XT  ■ max  { Xj^  , Xj^  } 


(5.9) 


The  equivalent  of  the  problem  A'  (4.1),  which  is  denoted  problem 
L'  in  the  inverse  design  space,  reads  as  follows 


f 


/ c . 


problem 


min  W(x.)  ■ £ '■ 

1 .X, 

I l 


with  £ C . . X.  < u . 
i iJ  1 J 


X . > X. 

1 —1 


X . < X . 

1 1 


(5.10) 

(5.11) 

(5.12) 

(5.13) 


The  filiation  of  the  various  problems  and  approximations  is 
illustrated  on  figure  5.1.  It  is  evident  that  an  exact  solution 
of  the  corresponding  approximated  problems  L or  A1  (or  L'  and  A') 
must  be  identical  if  one  recalls  that  these  approximate  problems 
are  convex  and  hence, that  there  is  no  danger  in  converging  to  a 
local  minimum. 


ThxA  the.  tqu4.v(U.tnce,  between  a Aotut^on  obtcuncd 

by  a G.O.C.  and  a iotatlon  obtoXnad  by  any  mathod  solving  exactly 
tka  CO Kae.6 ponding  paobtcmi  A1  ok  L.  The  4amc  conctiUKon  hotdi 
bctuieen  H,0.C.  and  paoblcmA  A*  ok  L*. 


i 


5.2.  Approximation  concepts  in  the  mathematical  programmin 


approach 


The  most  efficient  optimization  methods  based  on  mathematical 
programming  [ SCHMIT,  S18  ] use  high  quality  explicit  approxima- 
tions of  the  constraints  in  a design  space  where  they  are  as  shal- 
low as  possible.  In  the  present  context,  that  is  with  stiffness 
matrices  proportional  to  the  design  variables  a^,  using  the  space 
of  the  inverse  of  the  element  sizes  (5.1)  is  an  obvious  choice. 

In  the  spirit  of  the  classical  linearization  methods,  the  cons- 
traints are  developped  in  TAYLOR  serie  at  the  design  point  x? 
where  the  structure  has  been  analyzed 


V * I - <)  * 0 1 <«i  - 1 


f 

I 


rs. 


I 
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It  h«s  been  shown  in  (3.11)  to  (3. IS)  that  the  gradients  o£  the 
constraints  can  be  written  as 


3u. 


da . 

I 


C. . . (a.) 
ij  t 1 



do 


kl 


a . 

X 


da. 

X 


(5.15) 


with 


‘^iji 


■ (q . . K . q . ) a . 

’ij  1 1 


D..,. 

ikt 


(q..  K.  q.,)  a. 
^ik  1 ^xt  1 


or  in  the  inverse  design  space 


du . 

J—  - C . . (x . ) 
dx.  ij£  X 


do 


(5.16) 


Substitution  in  (5.14)  and  limiting  the  expressions  to  the 


terms  linear  in  x^  yields 


u.  - u?  - Z C? . X?  Z C? . X. 
jl  jl  ^ iji  1 ^ ijl  1 


(5.17) 


kt 


°kl  ■ ? '^ikt  *i  * ? *i 

1 1 


where  the  denotes  the  linearized  form  of  the  constraints. 


By  definition  of  the  coefficients  C...  and  D.,  . , one  has 
' ij£  ikl 


u?  - Z C?,  X? 
jt  ^ xji  1 


'kl 


Z D?.  X? 
ikl  1 


(5.18) 


and  hence  (5.17)  reduces  to 


(5.19) 


( 

In  classical  linearization  methods,  Che  objective  function  as 
well  as  the  constraints  are  linearized.  However,  as  mentionned, 

Che  simple  explicit  form  of  the  objective  function  makes  its 
linearization  useless.  The  linearized  form  of  Che  original  pro- 
blem R appears  to  be 


min 

w(Xi)  - / 

1 

(5.20) 

with 

u . . < u . 

J*  J 

(5.21) 

•v  - 

’kt  * “k 

(5.22)  ] 

X . < X . < X . 

—Ill 

(5.23)  i 

I 


i“l,n  j^lfUj  k*l,n  1*1, n^ 

which  is  identical  to  the  definition  of  Che  problem  L (5.4). 

It  i6  the.Kt£oKz  zvtdznt  that  the.  Aolutton  o£  the  pKoblem  Jl 
by  a Aequence  pKoblemA  L can  be  tnteKp/ieted  equivalently  aA 
a G.O.C.  approach  oK  aA  a mathematical  pacgaammlng  approach 
uAlng  a Apeclal  ^oA.m  the  llneaKlzatlon  method  In  the  InveKAe 
deAlgn  Apace. 

It  is  well  known  in  mathematical  programming  that  lineariza- 
tion methods  do  not  always  converge  and  if  they  do,  the  conver- 
gence is  not  often  monotonic.  This  is  due  Co  Che  fact  Chat,  at 
Che  optimum  of  the  linearized  problem,  Che  constraints  are  viola- 
ted Co  an  extend  that  depends  of  their  non  linearity.  This,  in 
general,  yields  a sequence  of  non  feasible  design  point.  Eventual- 
ly a restoration  of  a given  point  back  to  the  boundary  of  the 
feasible  domain  could  be  achieved,  by  a scaling  or  any  other  means. 
In  such  an  operation  however, the  weight  might  increase, which  cor- 
responds to  starting  an  eventually  diverging  process.  To  avoid 
this  difficulty  a classical  technique  is  to  avoid  to  be  driven 
too  far  away  from  the  exact  composite  constraint  surface.  This  is 
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achieved  easily  by  using  a nethod  Chat  remains, at  each  step, 
inside  or  at  the  boundary  of  the  linearized  feasible  domain. 

In  such  case. one  simply  stops  the  minimization  after  a limited 
number  of  steps, before  reaching  the  minimum  of  the  linearized 
problem.  What  is  obviously  a difficulty  is  to  evaluate  a priori 
Che  travel  that  can  be  achieved  without  violating  "coo  much"  the 
exact  constraints. 

Such  a concept  has  been  applied  by  SCHMIT  [ S18  ] with 

either  an  interior  point  penalty  method  or  a feasible  directions 
method.  It  will  be  applied  in  this  work  on  Che  basis  of  the 
projection  methods.  The  filiation  of  the  various  methods  is 
illustrated  in  figure  5.1. 
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6.  PRIMAL  PROJECTION  ALGORITHMS  FDR  SOLVING  THE  PROBLEMS 
L OR  L* 


6.1.  Fir*t  order  algorithm 

The  linear  constraints  (5.21)  (5.22)  are  written  in  matrix 
form, using  (5.19) 

X > - u (6.1) 


where  C is  an  (n  * m)  matrix  of  constants  (-  C?..  and  - D?,  .) 

ijt  ikf' 

evaluated  at  the  design  point  for  which  the  structure  has  been 
analyzed.  Besides  these  m linear  constraints,  we  have  also  2n 
side  constraints  corresponding  to  (5.22)  and  (5.23).  In  a first 
order  algorithm,  the  gradients  of  the  objective  function  g^  and 
of  the  constraints  are  required.  They  are  written 


p . I . 
1 1 


i • 1 , n 


(6.2) 


while  N denotes  the  matrix  of  the  gradients  of  the  active  cons- 
traints. The  method  considered  here  is  derived  from  the  projected 
gradient  method  of  ROSEN  [ R1  ] . Only  the  basic  results  contained 

in  ROSEN's  paper  are  recalled.  The  search  direction  is  given  by 


z--Pg--g+Nr  (6,3) 

where  P is  the  projection  matrix 

P - I - N (n'^^N)”^  n’^  (6.4) 


I and  where  r denotes  the  first  order  approximation  of  the  lagrange 

I multipliers 

si, 

i r - (N^  N)"^  n’^  g 


(6.5) 


r 


/o. 
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This  classical  algorithm  ( FI  ] can  not  be  applied  in  this  form 
since  the  size  of  the  matrix  N can  be  as  large  as  (n  x n)  which 
would  require  an  unacceptable  size  of  core  storage  in  a computer 
program.  It  can  be  modified  as  follows. 

Consider  first  the  elimination  of  the  side  constraints 
(5.22),  (5.23).  The  column  of  N corresponding  to  a given  constraint 
x^  ■ becoming  active,  contains  only  a factor  -*-1.  Let  A be  the 
matrix  deduced  from  N by  suppressing  the  line  and  column  correspon- 
ding to  the  constraint  and  the  variable  x^.  The  size  of  the  pro- 
jection matrix  P reduces  to  (n-1)  and  we  write 


- A (A^  A)“^  A^ 


n-1 


(6.6) 


To  obtain  more  informations  on  the  structure  of  the  matrix  P, 
let  us  rearrange  the  matrix  N such  that 


N - 


A 

0 

} 

% 

n 

B 

I* 

} 

n 

- 

- 

4 

m 

n 

(6.7) 


where  n denotes  the  number  of  free  variables  and  n the  number  of 

fixed  variables , that  is , constrained  to  be  constant  during  the 

• • • , \ 
Iteration.  The  q active  constraints  are  split  in  a group  of  m 

principal  (or  linear)  constraints  and  a group  of  n side  constraints. 
The  gradient  of  these  latters  are  simply  unit  vectors  ('*-1)  or  (-1) 
and  the  matrix  I*  contains  only  (-t-l)  or  (-1)  diagonal  terms,  depen- 
ding whether  the  constraint  is  a lower  or  upper  limit  (x.  or  x^) . 
Using  (6. 7), it  turns  out  that 


T - 1 
(N  N) 


(a'*^a)"^ 


-B*(A^A)"^ 


-(a’^a)“^b* 


I-+B*(A^A)"^B* 

n 


(6.8) 


with 


B*  • I*  B , and  finally 


T -IT 
(I-v-  - A(A  A)  A ) 
n 


(6.9) 


depends  only  of  the  matrix  A deduced  from  N by  suppressing  the 
lines  and  columns  corresponding  to  the  side  constraints. 

The  projection  relation  (6.4)  can  now  be  transformed  by 
using  the  same  partitionning  in  the  columns  g and  r 


r } m 
P 


Using  (6.8)  yields 


(6.10) 


/*T..-1  .T  -x. 
tp  - (A  A)  A g 

T*“  b*/»T..-1  .T'X'  t*/“  n \ 

r -I  g-B(AA)  Ag-I(g-B  r ) 

9 P 


(6.11) 


and  the  search  direction  becomes 


g + A r. 


(6.12) 


These  modifications  reduce  significantly  the  size  of  the  core 
storage  necessary  to  apply  efficiently  the  algorithm, since  the 
matrices  required  are  only 


(A*A) 


(n  X m ) 


(m  X m ) 


(6.13) 


where  m is  the  number  of  linear  (first  approximated)  constraints. 
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6.2.  Second  order  algorithm 


As  the  objective  function  is  separable,  a second  order  algo- 
rithm is  well  adapted.  Detailed  formulations  of  such  algorithms 
can  be  found  in  [ G3,  G5,  F9  ] . The  Hessian  matrix  F is  diagonal 

and  its  terms  are  given  by 


a^w  _ ^^i^i 
^ij  “ ax.  ax.  “ ^ij 


(6.14) 


i J 


X . 
1 


where  6^.^  is  the  Kronecker  symbol.  Its  inverse  is  therefore 
trivial  and  the  generalized  NEWTON  method  with  oblique  projections 
can  be  used  [ G3  ] . The  second  order  approximation  of  the  lagrange 

multipliers  turns  out  to  be 


T -1  -1  T -1 

r - (N^F  N)  N^  F ^ g 


(6.15) 


and  the  NEWTON  search  direction  is 


z - - P f"^  g - f"^  (-g  + Nr) 


(6.16) 


where  F is  the  projection  operator 


" -1  T -1  -1  T 

P - I - F N (N  F N)  N 


(6.17) 


The  computation  of  these  gener.l  expressions  is  significantly 
facilitated  if  the  same  partitioning  is  used, as  introduced  in 
(6. 7), not  only  for  the  matrix  N of  the  gradient  of  the  constraints, 
but  also  for  the  diagonal  Hessian  matrix  F so  that 


,-l 


»-l 


F-1 


The  projection  matrix  P turns  out  to  be 


p - 


(I-V. 

n 


'V'-l  T*''-!  -1  T 

F ^A(A  F ^A)  A ) 


(6.18) 


and  depends  only  of  the  n variables  that  have  not  reached  the 
upper  or  lower  limits,  or  x , The  lagrange  multipliers  take 
the  simplified  form 

. T *''-1  -1  T *''-1 

rp  - (A  F A)  ^ A^  F g 


(6.19) 


r ■ I*  (g  - B r ) 
s p' 


while  the  search  direction  becomes 


<\. 

z 


'V_1 

F ( 


g + A tp) 


(6.20) 


Z m 0 


This  form  requires  a reduced  amount  of  storage,  that  is  given  by 
(6.13)  plus  the  array  necessary  to  store  the  diagonal  Hessian 
matrix.  The  convergence  properties  of  a second  order  algorithm 
are  expected  to  be  significantly  better  than  the  first  order 
algorithm,  even  when  the  comparison  is  made  in  terms  of  computing 
time  rather  than  in  number  of  steps.  Examples  are  given  in  the 
following. 


6.3.  Auxiliary  algorithms 

A certain  number  of  auxiliary  algorithms  are  necessary  to  use 
either  one  of  the  basic  algorithms  describe'!  above. 


Conj^g£t£  £e£r£h_d^r£c£i£n£ 

The  basic  first  order  algorithm  of  ROSEN  is  significantly  improved 
if  the  search  directions  z are  selected  according  to  the  concept 
of  conjugate  directions  as  introduced  by  FLETCHER  [ F4  ) and 


a 


fek 


MIELE  [ M3  ] . The  algorithm  is  initiated  by  the  simple  projected 

gradient 


Z m - p o 

o o ®o 


The  subsequent  directions  are  extracted  from 


*k  " ” *k“l 


(6.21) 


with 


^k®k 


^k-l*k-l 


The  conjugation  algorithm  must  be  reinitiated  after  (n-q+l) 


Iterations . 


Restoration 


At  any  time  in  the  minimization^ the  current  design  point  has  to 
remain  on  the  composite  constraint  surface. 

A restoration  is  necessary  due  to  the  numerical  round-off  errors. 
In  practice, it  is  applied  each  time  the  set  of  active  constraints 
is  changed, and  each  time  the  conjugation  algorithm  is  reinitiated. 
The  restoration  consists  in  progressing  along  a direction  that  is 
normal  to  the  plane  tangent  to  the  constraints  at  the  current 
point.  If  denotes  a point  which  suffers  from  an  error,  the  cor* 
rected  position  x^  is  given  by 


-v  A / aT.  . -1 

X,  • X - A(A  A) 
i o 


(A^  i 7 . 1) 

o 


(6.22) 


for  the  first  order  algorithm,  and  by 


X '''-1  T '''-1  -1  T “V  T — — 

X - F A(A  F A)  (A  X ♦ B x ♦ u) 
o o 


(6.23) 
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for  the  newton's  algorithm  . It  is  evident  that  the  fixed 

variables  x do  not  have  to  be  restored, 
o 

I.ine£r_se^a£cli 

With  the  first  order  algorithm,  a cubic  interpolation  method 
has  been  used  [ Dl,  F21  ] while  in  the  second  order  algorithm  , 
advantage  has  been  taken  of  the  fact  that, in  such  methods, the 
optimal  step  is  close  to  unity  in  :he  neighborhood  of  a minimum. 
The  problem  can  be  presented  as  finding  the  minimum  of  the  func~ 
tion  of  the  step  length  t 


z 7 W(t)  - 


z . p . f . 

Z ■ 

i (x. -TZ . ) 
i 1 


(6.24) 


The  TCHEBYCHEV  method  is  used  with  the  initial  approximations 

“ 1 , which  would  be  exact  if  the  objective  function  was  quadra- 
tic. The  first  and  second  derivatives  of  (6.24)  are  easily  obtained. 
This  method  is  not  applied  in  the  projected  gradient  method  because 
the  initial  choice  of  the  step  length  reveals  rather  critical 
and  difficult  to  obtain. 

^tjra£e^y_f£r_s£l£cjt  i£g_the_ac  t£V£  £0£LSJt r£i n t£ 

Some  care  has  to  be  exercized  in  the  selection  of  the  active  cons- 
traints to  avoid  the  zigzaging  behavior  that  is  characteristic  of 
projection  methods.  The  approach  used  follows  the  suggestion  of 
ROSEN  [ R1  ] .It  consists  in  removing  a constraint  from  the  active 
set  only  if  it  leads  to  a sufficient  gain  in  the  decrease  of  the 
objective  function.  By  the  notion  of  marginal  cost, the  variation  of 
the  objective  function  can  be  deduced  from  the  change  in  the  value 
of  the  lagrange  multiplier  when  the  constraint  becomes  inactive 
[ FIO,  G5  ] . The  abandonment  of  the  j^^  constraint  leads  to  a 

gain  of  the  order  of 


where  Vj  is  the  diagonal  term  of  (N^N)  ^ or  (N^  F ^ N) 

Hence  the  strategy  consists  in  determining  among  the  constraints 
associated  with  negative  rj  components,  the  one  for  which  the 
gain  obtained  when  becoming  inactive  exceeds  y times  the  decrease 
of  the  objective  function  that  would  be  otherwise  obtained. 

For  a given  constraint , it  is  explicitly  written 


T 

z g 


(6.26) 


y being  a positive  constant.  For  first  order  algorithm  the  Hessian 
matrix  is  replaced  by  the  unit  matrix.  The  validity  of  this  test 
depends  of  the  quality  of  the  approximation  of  lagrange  multipliers. 
Therefore  it  is  applied  only  if 

g < E (6.27) 


where  e is  a given  tolerance.  This  insures  that  the  accuracy  on 
the  multipliers  is  sufficient  as  (6.27)  can  only  be  verified 
close  to  stationary  points. 


6.4.  Comparison  of  the  various  algorithms 

First  and  second  order  projection  algorithms  have  been  consi- 
dered for  the  solution  of  the  problem  L in  the  inverse  design 
space  x^.  In  the  original  design  space  a^,  dual  algorithms  of  the 
first  and  second  order  have  been  established  for  the  problem  A1 
in  chapter  4.  The  solutions  obtained  by  the  various  algorithms 
must  coincide  as  the  linearized  problem  L is  convex  and  therefore 
has  a unique  solution.  Two  special  examples  have  been  used  to  il- 
lustrate the  relative  performances  of  the  4 algorithms  in  competi- 
tion. The  comparison  is  made  on  the  CPU  time  for  programs  'written 
by  the  same  author  and  with  the  same  degree  of  generality.  This 
insures  that  all  the  aspects  of  the  computation  are  taken  into 
account  in  the  comparisons. 
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First  the  following  explicit  problem  has  been  solved 
1000 

min  W ■ £ a. 

1 

950  _ 1000 

with  u,  • £ — +10  £ 

^ 1 *i  951 

950  , , 1000 

u,  - £ lo"®  £ 

1 *i  951 

lo"^  < a.  < lo"^^ 

1 

The  analytical  solution  of  this  problem  is  easy  to  find  as 
a*  • 1 for  i - 1,950 

a*  - lo"^  for  i - 951,1000 

and  corresponds  to  a weight  W*  ■ 950.00005. 

The  convergence  of  the  4 algorithms  is  illustrated  on  figure  6.1. 

The  primal  algorithms  use  of  course  the  change  of  variables 

« 1/a^.  The  second  order  primal  algorithm  converges  much  more 
rapidly  than  the  first  order  one  which  was  stopped  after  394  ite- 
rations. In  general,  it  was  observed  that  the  fi.rst  order  algorithm 
requires  a number  of  steps  of  the  order  of  a multiple  of  the  number 
of  the  design  variables, while  the  second  order  algorithm  converges 
in  a number  of  steps  inferior  to  the  number  of  design  variables. 

Turning  to  the  dual  algorithms,  the  first  observation  is 
that  both  are  much  more  efficient  than  the  primal  ones.  This  con- 
clusion holds  usually  when 

/ 

n > > m (6.29) 

(here  n ■ 1000,  m "2).  The  advantage  of  using  the  second  order 


•=-  < 1000  (6.28) 
a . 

1 


— < 900 

a . 

1 
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algorithm  is  still  obvious  but  less  significant  (6  iterations 
instead  of  3) . 


The  test  seems  to  be  representative  of  the  problems  met  in 


structural  synthesis.  The  coefficients  c^j  of  the  constraints 
are  not  all  positive.  As  a consequence  some  of  the  design  varia- 


bles reach  their  lower  limit.  This  characteristic  reduces  the 
speed  of  convergence  of  the  primal  algorithms , but  not  of  the 
dual  ones. 


The  36  bar  isostatic  truss  of  figure  6.2  has  been  used  as 
a second  test.  A detailed  definition  of  the  problem  is  given  in 
table  6.3.  15  flexibility  constraints  are  imposed 


w. 


- w 


< e 


j - 1,  4 
k - 2,  5 


(6.30) 


and 


w . < w 
3 


1,  5 


among  which  only  5 are  independent. 

It  is  worth  pointing  out  that  the  classical  optimality  criteria 
( BIO,  G13,  T3  ] do  not  yield  the  optimum  of  this  problem.  As 
the  structure  is  isostatic, only  one  analysis  is  sufficient. 

The  comparison  of  the  performances  of  the  4 primal  and  dual  algo- 
rithms is  illustrated  on  figure  6.4.  The  test  is  different  from 
the  preceding  one, by  the  fact  that  the  number  of  constraints  is 
of  the  order  of  that  of  design  variables. 

Even  in  this  case  the  dual  methods  reveal  competitive  with  the 
primal  ones. 

The  conclusion  can  be  drawn  that  the  dual  algorithms  are  to 
be  recommended  in  general  when  the  number  of  variables  is  larger 
than  the  number  of  constraints.  With  primal  algorithms,  the  second 
order  brings  a significant  advantage , while  with  dual  algorithms 
the  benefit  is  counterbalanced  by  the  increased  complexity  of  the 
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solution  of  the  problem  L (or  L')  before  reaching  the  minimum.  j 

With  dual  methods  the  design  points  obtained  at  each  iteration 
are  not  feasible  points  and  Che  corresponding  weight  is  not,  in 
general,  decreasing,  as  it  is  the  dual  objective  function  which  i 

is  maximized.  Using  the  dual  algorithms  implies  therefore  to  I 

solve  completely  the  problem  A1  (or  L)  or  the  problem  A'  (or  L'). 


computation.  These  conclusions  have  been  verified  on  a lar 
number  of  applications. 


A significant  difference  still  exists  between  primal 
dual  algorithms,  which  will  be  used  at  the  advantage  of  th 
primal  ones  in  the  next  chapter  dealing  with  mixed  methods 
In  the  primal  algorithms  the  design  obtained  after  each  it 
is  a feasible  design, with  respect  to  the  approximate  const 
and  Che  weight  is  smaller  Chan  for  the  preceding  iteration 
It  is  therefore  possible, and  sometimes  advantageous , to  sto 


36  BAR  ISOSTATIC  TRUSS 
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aluminium 

E - 7200  hb 

p - 0.28  10  ^ k^/mm^ 

2 

£ ■ 0.5  mm 
a ■ 8 hb 

u - 10  mm  (at  nodes  1,  2,  3,  4,  5) 
z 

e ■ 0.5  mm 


Node  coordinates 


Node 

X 

Y 

Z 

Node 

X 

Y 

Z 

1 

25.00 

41.20 

0.00 

10 

1259.00 

266.4 1 

0.00 

2 

330. 75 

31.18 

0.00 

1 1 

1259.00 

264.43 

- 32.37 

3 

636.50 

2 1.18 

0.00 

12 

25.00 

548.00 

0.00 

4 

947.75 

11.00 

0.00 

1 3 

2 5.00 

516.01 

- 58.30 

5 

1259.00 

0.82 

0.00 

14 

636.50 

548.00 

0.00 

6 

25.00 

286.59 

0.00 

15 

636.50 

514.75 

- 60.60 

7 

25.00 

284.76 

- 29.91 

16 

1259.00 

548.00 

0.00 

8 

636.50 

276.59 

0.00 

17 

1259.00 

513.47 

- 63.95 

9 

330.75 

279.72 

- 30.52 

18 

947.75 

269.56 

- 31.75 

I 


Loading 


Node 

Direction 

Load  (daN) 

Node 

Direction 

Load  (daN) 

« 

Z 

10 

5 

Z 

10 

2 

Z 

- 20 

6 

z 

40 

3 

Z 

30 

8 

Z 

90 

4 

Z 

- 20 

10 

z 

40 

li 


TABLE  6.3 
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MIXED  METHODS 

^ . 1 . A projection  method  for  the  problem  R 

In  chapter  5,  it  has  been  shown  that  the  problem  R can  be 
replaced  by  a sequence  of  problems  L.  It  was  justified  as  a 
special  form  of  the  linearisation  methods.  Such  methods  are 
successful  in  that  case, due  to  the  fact  that  the  constraints  are 
near  to  linear  in  the  inverse  design  space.  On  the  other  hand 
the  sequence  of  problems  L corresponds  to  the  application  of  a 
generalized  optimality  criterion.  An  equivalent  solution  is  pro- 
vided by  the  dual  formulation  of  the  solution  of  a sequence  of 
problems  A1 . 

Considering  the  quasi  linear  characteristic  of  the  exact 
constraints  in  the  inverse  design  space,  one  could  contemplate 
a strict  application  of  the  projection  methods  directly  to  the 
problem  R.  Such  methods  have  the  fundamental  advantage  to  lead 

to  an,  at  least  local,  minimum  of  the  exact  problem  R without  any 
convergence  problem,  which  is  not  always  the  case  with  lineari- 
zation methods.  Another  advantage  of  the  projection  methods  is 
that  they  generate  a sequence  of  admissible  design  points  corres- 
ponding to  a sequence  of  decreasing  values  of  the  objective  func- 
tion. 


The  application  of  the  projection  methods,  like  that  of 
ROSEN  [ R2  ] , to  the  problem  R (instead  of  the  problem  L con- 

sidered in  chapter  6)  requires  the  following  steps. 

- At  a given  admissible  design  point,  analyse  the  structure 
and  compute  the  gradients  of  the  objective  function  and 
of  the  constraints; 


L- 

P 


Determine  a search  direction  by  a projection  of  the  gra- 
dient of  the  objective  function  on  the  plane  tangent  to 
the  active  constraints. 

Determine  a step  length  such  that  the  function  is  minimi- 
zed,but  without  violating  the  exact  constraints .which 
requires  a certain  number  of  structural  reanalyses; 


I 


i, 

I 

1, 

[■ 

!■ 


92. 


- As  the  constraints  are  non  linear,  a restoration  step  has 
to  be  applied  to  bring  the  design  point  back  on  the  com- 
posite constraint  surface,  which  is  a necessary  condition 
for  applying  the  projection  method.  This  step  requires 
also  to  reanalyse  the  structure  in  order  to  have  an  accu- 
rate estimation  of  the  constraints; 

- After  restoration  the  design  point  is  admissible  and  the 
process  is  repeated. 

It  is  well  known, and  it  has  been  illustrated  in  the  applica- 
tions of  the  projection  primal  methods  to  the  problem  L,  that 
the  number  of  steps  increases  in  proportion  to  the  number  of  de- 
sign variables.  As  each  iteration  requires  a certain  number  of 
structural  analysis,  the  application  of  such  methods  can  not  be 
considered  as  such  for  large  problems. 

Two  actions  can  be  taken  to  reduce  the  cost  of  computation 

- First,  reduce  the  number  of  minimization  steps; 

- Second,  reduce  the  number  of  structural  analyses  per  step 
or  for  a number  of  steps. 

The  concepts  used  in  establishing  the  projection  algorithms  for 
the  problem  L can  obviously  be  applied  here  to  the  problem  R. 
First, using  the  inverse  design  space  insures  to  have  nearly 
linear  constraints,  thus  allowing  larger  steps.  Second, the  expli^ 
cit  first  order  approximations  of  the  constraints  can  be  used  to 
compute  the  search  directions  and  the  step  lengths. 

Third, the  restoration  can  be  achieved  by  a simple  scaling  of  the 
design.  Such  a procedure  does  not  require  a reanalysis  and  it  has 
been  shown  in  section  3.2  that  it  preserves  the  value  of  the  cons- 
traints and  of  their  gradients 

^ 3u.  9u.  ^ 

“j  ■ “j  srl-  • 

^ • 1 e 


Consequently  the  analysis  after  scaling  is  the  only  one  which  is 
necessary  to  reinitiate  the  process. 
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The  algorithm  can  be  summarized  as  follows 


Analyse  Che  structure  at  an  admissible  point  x?; 

t-  - Compute  Che  explicit  first  approximations  of  the  constraints 

by  using  (3.9)  and  (3.20),  that  is,  in  the  x^  space,  linearize 
Che  constraints  in  the  form  (5.2),  (5.3) 


'V. 

u 


(7.2) 


Determine  the  search  direction  by  an  adequate  projection  of 
the  gradient  in  the  plane  tangent  Co  the  active  constraints 
and  compute  the  step  length.  These  operations  are  performed 
by  using  the  linearized  expressions  (7.2)  of  the  constraints 
and  thus  without  reanalyzing  Che  structure.  It  is  expected 
Chat  such  linearized  constraints  do  not  lead  to  a significant 
error  in  Che  step  length  so  Chat  Che  algorithm  still  works 
normally ; 


Restore  the  design  on  the  exact  constraint  surface  by  a simple 
scaling  by  a factor  f.  It  produces  a new  admissible  design 
point 


1 


X . 
1 


(7.3) 


at  which  an  analysis  is  achieved  and  Che  procedure  repeated. 
The  sealing  factor  is,  according  Co  (3.2) 

f ■ max  (f  , f ) 
u o 

with  ^ ^ (7. A) 
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The  algorithm  is  illustrated  in  figure  7.1.  If  it  is  compared 
with  the  projection  algorithms  proposed  to  solve  the  problem  L 
in  chapter  6,  it  appears  immediatly  that  the  only  difference 
is  that  here,  only  one  minimization  step  is  achieved  before  re- 
analyzing the  structure  and  thus  redefining  a new  set  of  approxi- 
mate constraints.  This  essential  difference  is  sufficient  to 
transform  the  linearization  method,  which  is  never  guarantied  to 
converge,  into  a primal  projection  method  which  is  known  to  mi- 
nimize the  weight  at  each  iteration  and  yield  to  convergence , if 
a sufficient  number  of  steps  is  achieved. 

7.2.  A mixed  method 

The  projection  method  outlined  in  the  preceding  section  has 
still  the  property  that  is  common  to  all  primal  methods,  that  is 
to  require  a number  of  iterations  proportional  to  the  number  of 
design  variables.  The  simplifications  introduced  have  reduced  the 
cost  per  iteration  but  have  not  changed  this  basic  property. 

It  is  illustrated  on  the  10  bar  truss,  in  the  configuration 
defined  geometrically  by  the  figure  4.5.  The  application  of  the 
projection  algorithms  described  above,  using  the  first  order 
(conjugate  gradient)  or  the  second  order  (NEWTOM's  method)  appro- 
ximations of  the  objective  function,  leads  to  the  results  illus- 
trated by  figure  7.2.  The  convergence  is  characteristic  of  the 
strict  applications  of  primal  mathematical  programming  methods, 
that  is  extremely  slow  with  the  first  order  algorithm. 

With  the  second  order  algorithm  the  number  of  steps,  and  thus  of 
structural  analyses,  is  significantly  reduced, but  still  a multiple 
of  the  number  of  design  variables.  In  figure  7.3,  the  performances 
of  two  classical  optimality  criteria  [ G13,  T3  ] and  of  the  con- 
sistent , first  order  generalized  optimality  criterion,  defined  in 
chapter  3,  are  compared.  The  most  sophisticated  optimality  crite- 
rion was  applied  via  a dual  formulation  as  described  in  chapter  4. 

The  comparison  of  the  two  applications  of  the  projection 
methods,  to  the  sequence  of  problems  L and  to  the  problem  R, 
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suggest!  a mixed  method  which  is  expected  to  retain  some  of  the 
advantages  of  both  approaches. 

The  simple, but  essential, idea  is  to  allow  for  a certain 
number,  denoted  k,  of  minimisation  stips  before  reanalyzing 
the  structure  and  updating  the  linearized  approximations  of  the 
constraints.  When  k is  limited  to  one  step,  the  method  is  a 
strict  application  of  the  projection  algorithms  to  the  problem  R 
and  presents  all  the  properties  of  strict  primal  mathematical 
programming  approaches,  that  is,  high  cost  but  guarantied  conver- 
gence. When  k is  not  limited,  the  problem  is  identical  to  the 
problem  L (S.A).  It  is  solved  completely  before  reanalyzing  the 
structure  and  defining  the  next  problem  L.  This  method  has  been 
interpreted  as  a linearization  method  in  the  inverse  design  space, 
or  as  a generalized  optimality  criterion  approach.  The  corres- 
ponding properties  of  fast,  but  uncertain  convergence  are  to  be 
expected.  When  k is  limited  to  a given  finite  number,  the  method 
is  mixed.  It  is  evident  that  the  larger  the  number  k is,  the 
cheaper  the  application  will  be,  as  the  number  of  reanalyses  is 
reduced.  The  reason  of  the  high  efficiency  of  the  optimality  cri- 
teria is  now  clearly  related  to  the  fact  that  the  explicit  appro- 
ximations made  for  the  constraints  are  not  updated  before  reaching 
the  solution  of  the  corresponding  approximate  problem.  This  consti- 
tutes the  power  and  the  weakness  of  the  optimality  criteria  approa- 
ches. 


The  figure  7.4  illustrates  the  mixed  method. 

At  the  design  point  x°,  the  structure  is  analyzed  and  the  set  of 
linearized  constraints  is  computed  (figure  7.4  (a)). 

After  scaling  (figure  7.4  (b))  a certain  number  of  minimization 
steps  are  achieved  along  the  surface  u * u,  leading  to  a point 
x^  or  x^.  The  point  x^  corresponds  to  the  minimum  weight  under 
the  approximate  constraint,  while  x*  is  the  minimum  weight  under 
the  exact  constraint.  A scaling  by  a factor  f brings  back  the 
point  on  the  exact  constraint  surface  in  X2,  where  the  weight 
is  smaller  than  in  x^,  but  larger  then  in  x*.  A smaller  step,  or 
in  a n-dimensional  space,  a smeller  number  IT  of  steps,  could  be 


selected  leading,  say  in  x|,  which  is  not  the  minimum  weight 
under  u ■ u.  The  scaling  by  the  factor  £'  < £ yields  the  admis- 
sible point  at  which  the  weight  is  smaller  than  in  X2« 

The  figure  7.5  (a)  and  (b)  illustrates  the  reason  why  the 

linearization  methods  (or  the  generalized  optimality  criteria) 

diverge  in  some  cases.  In  figure  7.5  (a)  the  constraint  u - u 

is  rather  shallow.  It  is  characteristic  of  weak  hyperstaticity 

as,  in  an  isostatic  structure,  it  would  be  exactly  linear.  The 

solution  of  the  problem  L (that  is  k unlimited  in  the  mixed 

method)  yields  the  minimum  at  x^^  which  after  scaling  comes  in  X2 

where  the  weight  W(x2)  is  smaller  than  W(x^) . In  figure  7.5  (b) 

the  non  linearity  of  the  exact  constraint  is  more  pronounced  and 

the  process  diverges  as  W(x»)  > U(x  ).  A limitation  of  k would  be 

z o 

necessary  to  keep  the  process  converging. 

Thli  InttKpKttcutlon  attorn  to  comtdeA.  the.  numben.  oi  itep^ 

T.  ah  a convergence  controt  parameter  that  hhoutd  be  ahhtgned 
htgh  vatueh  £or  economy  and  reduced  tahen  divergence  occurh. 

The  effect  of  the  Ic  convergence  control  parameter  is  illus- 
trated in  Fig.  7.6  a & b for  the  10  bar  truss  defined  in  Fig.  4.5. 
It  is  recalled  that  all  the  constraints, including  therefore  those 
on  the  stresses, are  exactly  linearized  in  this  problem.  The  in- 
crease in  convergence  speed  for  an  increase  of  the  number  of 
steps  k is  very  significant,  especially  using  the  conjugate  gra- 
dient algorithm.  For  k > 50  (or  k > 4 in  the  NEWTON's  method) 
the  problems  L are  exactly  solved  and  the  performances  become 
similar,  as  it  corresponds , in  fact, to  that  of  the  generalized 
optimality  criterion,  which  has  been  presented  in  the  chapter 
dealing  with  dual  formulations.  The  table  7.7  gives  the  values 
of  the  weight  and  of  Che  scaling  factor  using  the  conjugate  gra- 
dient algorithm  for  k ■ 1,5  and  10.  It  shows  that  the  scaling 
factors  remain  closer  to  unity  when  k is  small,  and  also  that 
Che  scaling  remains  moderate  for  k <■  10,  which  indicates  that 
the  exact  constraints  are  not  seriously  violated  after  10  steps. 


7.3.  Extension  of  the  mixed  method 


In  chapter  5 we  have  considered  a sequence  of  problems  L' 
to  approximate  the  problem  R.  The  problem  L'  defined  in  (3.10)  in 
the  inverse  design  space, corresponds  to  the  problem  A'  defined 
in  (4.1)  in  the  direct  design  space.  These  problems  have  been 
shown  to  correspond  to  the  hybrid  optimality  criteria,  introduced 
in  section  3.6,  and  characterized  by  the  use  of  zero  and  first 
order  approximations  of  the  constraints.  It  is  recalled  chat  a 
first  order  approximation  of  all  the  stress  constraints  can 
reveal  very  costly  for  the  analysis, due  to  Che  need  to  incorporate 
a large  number  of  virtual  loading  cases.  The  zero  order  approxima- 
tion of  Che  stress  constraints  by  F.S.D.  does  not  introduce  any 
increase  in  the  analysis  cost. 

Indeed  Che  zero  order  approximation  for  the  stress  constraints 
can  reveal  satisfactory  if  the  conditions  (3.38)  are  met.  These 
conditions  were  proposed  as  a basis  for  selecting  Che  potentially 
critical  stress  constraints  requiring  an  exact  linearization  in 
Che  context  of  the  generalized  optimality  criteria  approach. 

In  the  preceding  section, it  has  been  shown  that  Che  possible 
divergence  of  Che  various  algorithms  considered  is  related  to 
essentially  two  factors 

- Che  quality  of  the  explicit  approximations  of  Che  constraints; 

- the  number  k of  minimization  steps  achieved  with  a given  set 
of  approximate  constraints. 

This  implies  Chat, in  Che  search  for  the  most  efficient  minimiza- 
tion algorithms,  these  two  factors  should  be  taken  into  account. 

It  suggests  Co  consider  the  partial  solution  of  Che  problem  L' 
using  Che  same  concepts  as  developed  in  the  preceding  section. 

Doing  so  the  method  becomes  mixed  in  two  senses 

- zero  order  approximations  are  used  instead  of  the  first  order 
required  by  a strict  application  of  the  projection  methods  to 
Che  problem  R; 
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- more  than  one  step  of  minimization  is  achieved  with  a given 
set  of  approximate  constraints. 

The  benefit  expected  from  such  a technique  is  to  reduce  the 
analysis  cost  by  avoiding  many,  if  not  all,  the  virtual  loading 
cases  associated  with  higher  quality  approximations.  The  draw- 
back  is  a possible  reduction  of  the  number  k of  steps  that  can 
be  allowed  and  still  avoids  the  divergence.  The  situation  is 
evidently  less  critical  than  in  using  the  various  optimality 
criteria  as  the  problem  L'  is  not  necessarily  solved  in  the  mixed 
method.  The  relations  between  the  various  methods  (linearization, 
mixed  or  by  optimality  criteria)  are  illustrated  in  figure  5.1. 

An  example  is  given  in  figure  7.8  where  the  10  bar  truss 
defined  in  figure  4.5  is  used  again.  The  stress  constraints  are 
all  treated  by  F.S.D.  while  the  displacement  constraints  are 
still  necessarily  first  order  approximated. 

The  convergence  has  been  investigated  for  various  values  of  the 
k convergence  control  parameter,  using  the  first  and  second  order 
algorithms  (conjugate  gradient  and  generalized  NEWTON).  The  effect 
of  k is  striking.  For  k unlimited  (k^>  100  or  k2  ^ 5)  the  instabi- 
lity observed  when  applying  the  classical  optimality  criteria  re- 
appears (compare  with  figure  7.3).  Limiting  k to  k^^  > 40  or  k2  * 2 
brings  back  the  monotonous  convergence.  It  is  however  important 
to  note  that  the  generated  design  corresponds  to  a weight  of 
5076.7  lbs  which  is  slightly  heavier  than  the  best  known  optimum 
of  5060.85  lbs  obtained  in  figure  7.2  and  in  figure  7.6. 

It  indicates  the  presence  of  a local  minimum  close  to  the  global 
one.  As  is  well  known,  there  is  no  way  of  insuring  to  reach  the 
global  minimum  and  the  convergence  of  one  algorithm  or  the  other 
to  a local  minimum  can  be  due  to  a specific  peculiarity  of  the 
problem. 

7.4.  Controlling  the  convergence 

If  the  k convergence  control  parameter  allows  effectively  to 
insure  the  convergence  when  a sufficiently  small  value  is  taken. 
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the  determination  of  its  highest  admissible  value,  that  still 
yields  to  convergence,  is  not  possible  a priori.  When  divergence 
occurs,  the  normal  procedure  is  to  restart  the  last  minimization, 
with  a reduced  number  of  allowed  steps.  This  is  an  a posteriori 
measure, which  is  closely  related  to  the  concept  of  "intermediate 
design  vectors"  introduced  by  KHOT  [ KIO  ] . 

The  limitation  of  Ic  would  be  a better  procedure  if  applied  as 
an  a priori  measure. 

The  intermediate  design  vector  technique  proposed  in  [ KIO  ] 
considers  two  design  points  a^  and  a^  corresponding  to  the  results 
at  two  successive  steps.  If  the  process  diverges 

W(a.)  > W(a)  (7.5) 

i o 

The  illustration  of  figure  7.8  (a)  suggests  that  between  the 
scaling  lines  and  D(a^)  there  exists  an  optimal  scaling 

line  D(a*)  such  that 

W(a*)  < W(a^)  < W(a.)  (7.6) 

o 1 

One  such  point , denoted  a, is  obtained  by  determining  the  interme- 
diate scaling  line  D(a)  as  a normalized  average  between  D(a^) 
and  D(a^),  and  next  by  computing  the  intersection  of  D(a)  with 
the  exact  composite  constraint  surface.  This  procedure  is  applied 
in  the  direct  space  by  KHOT  [ KIO  ] each  time  the  weight  increases 
in  two  successive  steps.  It  was  proposed  in  the  context  of  the 
classical  optimality  criteria.  It  relies  on  the  implicit  assump- 
tions that  the  composite  constraint  surface  is  convex. 

The  same  technique  can  be  applied  with  advantage  in  the 


inverse  design  space,  in  conjunction  with  the  mixed  method  and  is 
illustrated  in  figure  7.8  (b) . Starting  at  an  admissible  point  x^, 
the  solution  of  the  linearized  problem  yields  such  that 
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and  it  could  happen  that  WCx^)  > W(x^)  which  leads  to  divergence. 
A family  of  intermediate  design  vectors  can  be  generated  in  the 
space  containing  D(x^)  and  D(x^).  It  is  written 

X - x^  + T (Xj  - x^)  0 < T < 1 (7.7) 

and  produces  the  restored  design  points 

•v. 


where  f is  the  scaling  factor. 

A simple  procedure  would  be  to  transpose  in  the  inverse  design 
space  the  procedure  of  KHOT  [ KIO  ] and  write 


X 


X + X, 

o 1 


(7.8) 


A better  procedure  is  to  achieve  a linear  search  along  the  direc- 
tion (xj^  - x^)  for  determining  the  optimum  intermediate  design 
vector.  It  is  written 


min  W(t)  ■ f(T)  W(t) 


(7.9) 


Suppose  first  that  the  same  constraint  u u remains  active 
for  all  the  x possible  intermediate  designs.  Use  can  be  made  of 
the  knowledge  of  Che  gradients  of  the  constraints  at  both  x^ 
and  x^  points.  The  sc'hling  factor  and  its  derivative  are  given  by 


f(o  - uUiiil 

u 

(7.10) 
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We  have  however 


u°  - L C?  X? 

i ^ " 


- E 

i " ^ 


r°  m 

1 9x . 
1 


(7.11) 


and  therefore 


fo  - f(o)  - ^ - 1 
u 


- f(l)  - — 


from  which  we  derive 


f;  - i (E  C°  %[  - u°)  - 0 
u , i 

o 1 . 

as  X and  x satisfy  the  linearized  constraint 


E C . X . * u 
i " " 


f\  - - (u^  - E X?) 

V u i ^ ^ 

V • 

The  function  fCt)  can  be  approximated  by  a cubic  polynomial 


f(T)  » (t[  - 2 f^  ♦ 2)  T 


(7.12) 


. ♦ (3  f^  - f • - 3)  ♦ 1 


which  is,  in  practice,  a very  accurate  representation  for  such  a 

function.  For  practical  applications,  one  has  to  take  care  that 

the  same  constraints  do  not  necessarily  remain  active  between  x° 

1 ‘ *■ 

and  X and  that,  consequently,  the  scaling  factor  has  to  be  written 
as  the  envelope 
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f(T)  - max  I f.(T)  1 j - 1,  n,  (7.13) 

j ^ 

taken  over  the  imposed  constraints.  In  addition  when  a combi- 
nation of  zero  and  first  order  approximations  of  the  constraints 
are  used,  polynomials  of  lower  degree  have  to  be  used  for  f(T) 
in  (7.12)  when  the  gradients  (7,11)  are  not  known  at  one  or  at 
both  ends  of  the  segment. 

This  more  sophisticated  procedure  for  selecting  intermediate 
design  vectors  in  the  inverse  space  has  been  applied  to  the  10 
bar  truss  defined  in  4.5.  Cubic  approximations  of  f ( t)  have  been 
used  for  the  flexibility  constraints,  while  linear  ones  have  been 
taken  for  the  stress  constraints  which  are  all  treated  by  F.S.D. 

The  results  are  compared  to  those  obtained  by  KHOT  ( KIO  ] in 
figure  7.9. 

The  proposed  procedure  had  to  be  applied  only  once,  between  itera- 
tions 7 and  8,  and  reveals  very  efficient  as  the  minimum  of 
5076.77  lbs  was  obtained  rapidly  without  any  additional  search 
for  intermediate  designs.  The  advantage  is  important  over  the 
more  intuitive  procedure  in  the  direct  space. 

It  remains  however  that  the  use  of  intermediate  design  vectors 
is  only  a palliative  measure.  A better  answer  would  be  given  by  an 
a priori  estimation  of  the  maximum  number  of  steps  k that  can  be 
achieved  in  the  coming  iteration.  It  seems  that  such  a forecast 
could  be  achieved  on  the  basis  of  an  extrapolation  of  the  scaling 
factors  affecting  the  various  constraints.  It  certainly  deserves 
further  study. 


0 


50 


100 


150 


FIG.  7.2  ; 10  BAR  TRUSS-PROJECTION  METHODS 


I 

J 


i 


CONVERGENCE  CONTROL  IN  THE  MIXED 

- - - - --  1, 

'J 


METHOD 


I 


FIGURE  7.4 


WEIGHT  (lb) 


FIGURE  7.6.B.  10  BAR  TRUSS  - THE  MIXED  METHOD  CONCEPT 


110 


k 

- 1 

k 

- 5 

k - 

10 

iter at  ion 

WEIGHT 

SCALE 

WEIGHT 

SCALE 

WEIGHT 

SeVLL 

(lbs) 

FACTOR 

(lbs) 

FACTOR 

(lbs) 

FACTOR 

1 

8266 

1.0 

8266 

1.0 

8266 

1.0 

2 

6321 

0.9996 

6337 

0.9924 

6203 

0.9002 

3 

6626 

0.9993 

6153 

0.9907 

597i) 

0.9364 

4 

6541 

0.9993 

6060 

0.9971 

5306 

0.9923 

5 

6497 

0.9995 

5939 

0.9974 

5766 

0.9925 

6 

6456 

0.9994 

5939 

0.9935 

5651 

0.9925 

7 

6425 

0.9996 

5394 

0.9995 

5492 

0.9364 

3 

6306 

0.9996 

5853 

0.9983 

5396 

0.9920 

n 

6373 

0.9997 

5810 

0.9972 

5306 

0.9975 

10 

6349 

0.9997 

5791 

0.9993 

5259 

1.0291 

11 

6329 

0.9997 

5772 

0.9999 

5133 

()."995 

12 

6307 
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5756 
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5113 

1.O035 

13 
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5736 
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14 
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16 
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0.9937 
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13 

6086 

0.9932 

5625 
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5077 

1.0)00 

19 

6073 
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0.9993 

20 

6053 

0.9996 

5573 

0.9993 

21 

6041 

0.9999 

5567 

0.9992 

22 

6036 

0.9998 

5530 

0.9994 

23 

6031 

1.0000 

5510 

0.9996 

24 
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25 
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29 
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O.P999 
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J1 

5903 
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32 
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0.9P99 

5443 

1 .0000 

33 

5937 

1 .0000 

5436 

1.0000 

34 

5932 
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5427 

0.9993 

35 
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0.9999 

5421 

1.0000 

10  BAR  THUSS 

EFFECT  OF  THE  K CONVERGENCE  CONTROL  PARAMETER 
(First  order  algorithm) 


Table  7.7 
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8.  NUMERICAL  APPLICATIONS 
i ■ 

I 

I 

The  applicaClons  have  been  divided  in  3 groups.  The  first  one 
concerns  trusses^ which  have  the  advantage  of  yielding  rather  pure 
problems  for  which  the  optimization  results  are  not  influenced  by 
the  modeling  technique.  The  second  group  of  applications  contains 
box  beam  problems  and  the  third, more  complex,  near  to  industrial 
problems . 

Two  different  programs  have  been  used.  One  is  an  experimental 
program  using  only  bar  elements  and  working  in  core.  It  is  denoted 
OPTBAR.  It  allows  to  use  all  the  algorithms  described  in  this 
work,  that  are  based  on  primal  or  dual  methods,  using  first  or 
second  order  minimization  techniques.  It  also  allows  to  select 
automatically  the  most  critical  stress  constraints,  according  to 
test  (3.38),  by  specifying  two  tolerances  and  £2  corresponding 
to  the  two  parts  of  the  test. 

The  second  program  is  the  optimization  module,  named  OPTIM, 
that  has  been  written  for  the  SAMCEF  general  purpose  finite  ele- 
ment program  [ S19  ] . It  allows  also  for  primal  and  dual  approa- 

ches, using  first  and  second  algorithms,  but  does  not  include  the 
automatic  selection  of  the  critical  stress  constraints.  For  compa- 
rison purposes,  the  envelope  method  of  GELLATLY  [ G13  ] is  also 

incorporated . 

The  following  abbreviations  are  used 

- C.O.C  Classical  optimality  criteria,  as  proposed  in  the 

litterature ; 

- G.O.C  Generalized  optimality  criteria.  It  implies  that 

all  Che  active  constraints  are  linearized  and  chat 
the  corresponding  problem  L (or  Al)  is  solved 
exactly  by  a primal  (or  dual)  method; 
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- H.O.C 


- S.O.C 


- MMl(k-xx) 


} 

- MM2(k>xx) 


Hybrid  optimality  criteria.  It  implies  that  a 
limited  number  of  stress  constraints  are  linearized 
in  addition  to  the  flexibility  constraints. 

The  approximate  problem  L*  (or  A')  is  solved 
exactly  by  a primal  (or  dual)  method; 

Simple  optimality  criteria.  All  the  stress  cons- 
traints are  zero  order  approximated  by  F.S.D.  while 
the  flexibility  constraints  are  linearized. 

This  corresponds  often  to  what  the  C.O.C  claim  to 
achieve.  The  corresponding  problem  A'  (or  L')  is 
however  solved  exactly  by  dual  (or  primal)  methods; 

The  mixed  method  is  used  with  first  (MMl)  or 
second  (MM2)  order  projection  algorithms.  These 
abbreviations  can  follow  G.O.C,  H.O.C  or  S.O.C. 

It  means  that  the  corresponding  problems  (L  or  L') 
are  not  solved  completely  and  therefore  Che  method 
is  no  longer  equivalent  to  an  optimality  criteria 
method.  The  convergence  control  parameter  has  been 
fixed  Co  k ■ XX  and  the  approximations  of  Che  cons- 
traints are  chose  defined  by  the  abbreviations  G.O.C, 
H.O.C,  S.O.C. 


Unless  otherwise  specified, 
uniform  distribution  of  the 


Che  starting  point 
design  variables. 


corresponds 


to  a 
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8.1.  Truss  problems 


8.1.1.  Simple  cantilever  10  bar  truss 


This  simple  and  classical  example  has  aircfdy  been  mencion- 
ned  to  illustrate  various  aspects  of  the  methods<'proposed  in  the 
preceding  chapter.  The  geometry  is  defined  in  figure  8.1. 


^t£e£s_c£n£t£a£n^s_onl^  (Problem  l.A)  ’■ 

'"'I 

This  problem  has  been  examined  in  chapter  4 and  ,t(ie  results  illustra- 
ted by  figure  4.6.  It  is  recalled  that  vhen  the  stress  limit  is 
25000  psi  in  all  the  bars,  the  F.S.D.  yields  the  optimum. 

When  the  limit  in  bar  nbr  8 is  increased,  the  F.S.D.  reveals  unable 
to  yield  the  optimum  from  Og  > 36^800.  Between  36,800  and  45,000  psi 
the  F.S.D.  does  not  converge  in  a reasonable  number  of  steps. 

When  the  linearization  of  stress  constraints  is  introduced,  the 
optimum  is  easily  obtained  for  any  value  of  the  Og  limit. 


jPr£b2.em 

The  loads  and  the  constraints  are 


F at  nodes  2 > - 100  K 

y 


3 - - 100  K 


|u^|  at  all  nodes  < 2. in 


The  results  obtained  for  this  problem  are  summarised  in  tables 

* \ 

8.2  and  8.3.  Convergence  curves  have  been  presented 


- for  the  6.0.C  in  figure  7.3; 


- for  the  MM2(k>2)  in  figure  7.6  where  certain  rj^sults  extracted 
from  the  litterature  on  C.O.C  have  also  been  plotted; 


- for  the  S.O.C  using  the  intermediate  design  vectdr  technique, 
in  figure  7.9. 
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The  essential  difficulty  in  this  problem  comes  from  Che  element 
nbr  6.  The  stress  constraint  in  this  element  is  the  only  one 
that  reveals  critical  and  its  linearization  is  necessary  to  obtain 
Che  lowest  weight  (5061  lbs). 

This  solution  has  been  obtained  by  a H.O.C,  that  is  by  using  F.S.D. 
for  all  the  elements,  but  element  6.  It  seems  that  this  design 
can  not  be  reached  by  the  other  methods  which  converge  to  the  local 
optimum  existing  in  the  neighborhood  (5076  lbs) . 

■ 

P^r£bJ_em  l,.^ 

The  constraints  are  identical , but  the  loads  are 

- 150  K 

- 150  K 

+ 50  K 

+ 50  K 


The  stress  constraint  in  element  6 is  more  critical  than  in 
the  preceding  problem, which  leads  to  a more  severe  divergence 
of  the  C.O.C  and  of  the  S.O.C  and  in  fact  for  any  method  relying 
on  F.S.D  in  element  6.  Convergence  curves  are  given 

- for  the  S.O.C,  using  MM2(k  ■ 1,2  and  ••  ),  in  figure  8.4; 

- for  the  H.O.C,  using  stress  constraint  linearization  in  element 
6 only,  for  the  envelope  method  of  GELLATLY  ( G13  ] which  is 
also  available  in  the  SAMCEF  program,  and  for  various  reference 
works,  in  figure  8.5. 

The  results  are  presented  in  tabular  form  in  tables  8.6  and  8.7. 
The  best  results  are  those  of  the  H.O.C  closely  followed  by  those 
of  SCHMIT  ( S18  ] . 

Problem  l.D 


The  loads  are  applied  in  nodes  2 and  3,  as  in  problem  l.B,  but 
equality  constraints  are  imposed  at  nodes  4 and  5 


r 
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The  problem  has  been  solved  by  Che  H.O.C  using  a linearization 
of  Che  stress  constraint , only  in  Che  element  3. 

Note  Chat  the  procedure  of  automatic  solution  of  these 
stress  constraints  that  require  1 inearization, has  been  used 
and  Che  two  tolerances  for  the  test  (3.38)  are  given 

on  the  figure. 

The  results  are  given  in  digital  form  in  table  8.8. 

It  is  observed  that  the  convergence  is  obtained  in  8 iterations 
and  Chat  the  equality  constraints  are  satisfied  at  the  last  ite- 
ration only, while  the  stress  in  element  3 has  reached  its  upper 
limit.  This  problem  was  suggested  by  VENKAYYA  ( V7  ] but  no 
reference  solution  was  found  in  the  litteraCure. 
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8.1.2.  Double  cantilever  22  bar  truss 

As  mentionned  in  chapCer  3,  Che  FSD  leads  Co  poor  resulcs 
when  all  Che  sCrucCural  componencs  do  noC  have  Che  same  scress 
limic,  which  is  Che  case  when  differenC  maCerials  are  considered. 
BERKE  ( BIS  ] proposed  a special  problem  Co  illuscrace  Che  diffi- 
culCy.  Ic  is  made  of  Cwice  Che  10  bar  cruss  of  seccion  8.1.1, 
connecCed  by  Cwo  bars  as  indicaCed  on  figure  8.9.  The  daCa  and 
Che  conscraincs  are  given  in  Cable  8.10.  Three  cases  have  been 
examined . 

£rob2_em 

3 

All  Che  elemenCs  have  Che  same  mass  densicy  (0.1  Ibs/in  ). 

The  sCress  limic  is  Che  same  for  all  Che  elemenCs  (o  - 23000  psi) 
buc  for  Che  2 bars  connecCing  Che  10  bar  Crusses  Co  the  load,  for 
which  A scress  limit  of  SOOOOOpsi  is  imposed. 

In  this  case  Che  FSD  yields  Che  opCimum  in  46  iceracions. 

Using  Che  G.O.C.  does  not  improve  Che  convergence.  The  Cable  8.11 
presence  Che  design  obCained^ while  the  table  8.12  contains  Che 
evolution  of  the  weight  for  various  configurations  and  methods. 

^r£b^em 

The  stress  limit  in  the  left  part  of  Che  truss  is  increased  up  Co 
50000  psi.  The  optimum  design  in  this  case  is  not  FSD  and  Che 
stress  ratio  algorithm  diverges.  Using  the  H.O.C.  with  automatic 
selection  bf  the  elements  requiring  stress  constraint  linearization 
according  cb  (3.38)  leads  Co  Che  linearization  of  10  stress  cons- 
traints. This  number  could  probably  be  reduced  by  using  more  severe 
tolerances.  The  optimum  design  obtained  by  the  H.O.C.  is  given  in 
cable  8.13  and  corresponds  to  a weight  of  654  lbs.  It  differs  dras- 
tically from  Che  heavier  (1044  lbs)  FSD  in  which  the  larger  sec- 
tions are  given  to  the  wrong  side  of  truss.  The  evolution  of  the 
weight  is  recorded  in  Cable  8.12  and  illustrated  in  figure  8.14  (a). 
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P^ro^b^em 

This  case  differs  from  case  B by  a higher  mass  density  in  the 
left  part  (p  > 0.4  Ibs/in  ) which  also  has  the  higher  stress 
limit.  The  data  are  given  in  table  8.10.  In  thu  F.S.U.  approach, 
the  sections  in  the  left  part  are  assigned  smaller  values  due  to 
their  higher  allowable  stress.  So  the  F.S.D.  brings  the  material 
back  to  the  left  side,  but  it  does  this  for  a wrong  reason  ! 

The  results  obtained  by  S.O.C.  and  by  H.O.C.  are  presented 
in  table  8.12  and  in  figure  8.14  (b) . The  optimum  is  a F.S.D. 
but  it  is  obtained  much  more  rapidly  by  the  H.O.C.  than  by  the 
simple  F.S.D.  used  in  the  S.O.C.  Note  that  all  the  elements  reach 
either  the  maximum  stress  or  the  minimum  size  and  that  the  design 
is  completely  reversed  from  that  obtained  in  case  B.  The  final 
design  is  detailed  in  table  8.15. 
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This  classical  problem  [ V2  ] is  defined  in  figure  8.16 
and  in  cable  8.17.  The  symmetry  of  Che  problem  is  used  Co  reduce 
Che  number  of  design  variables  Co  105.  Each  design  variable  cor- 
responds Co  a group  of  two  bars. 


Pr£b^em  _3.A 

No  stress  constraints  are  imposed.  The  S.O.C,  H.O.C  and  G.O.C. 
are  identical  in  this  case.  They  differ  from  C.O.C.  by  their 
rigorous  treatment,  by  dual  formulation,  of  the  explicit  problems 
jenerated  after  each  analysis.  The  figure  8.18  presents  Che  con- 
vergence of  the  weight, given  in  digital  form  in  table  8.19.  This 
Cable  also  contains  the  computing  time  for  Che  envelope  method 
( G13  1 and  Che  second  order  dual  algorithm.  They  have  been  ob- 

tained for  a programming  that  presents  the  same  level  of  sophis- 
tication in  both  cases.  It  shows  that  Che  more  general  and  rigo- 
rous dual  methods  are  not  significantly  more  costly  than  the  simple 
and  intuitive  envelope  method.  The  optimum  design  obtained  by  Che 
H.O.C.  is  presented  in  table  8.20.  Its  weight  (29037  lbs)  is 
smaller  Chan  that  reported  by  VENKAYYA  [ V2  ] (31020  lbs). 


Stress  constraints  (o  < 10000  psi)  are  added  to  the  preceding 

problem.  The  S.O.C.  is  now  less  efficient  than  Che  H.O.C.  j 

il 

The  results  are  illustrated  in  figure  8.21  and  in  tables  8.19  j 

and  8.20.  Note  Chat  Che  final  design  weights  29163  lbs  and  is  j 

not  very  different  from  Che  preceding  one.  j 
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8.1.4.  Wing  carry  through  box  (63  bare] 


This  problem  is  also  classical  [ B8  ] . It  is  defined  in 

figure  8.22  and  cable  8.23.  Two  cases  are  examined. 

P^r£b^em  ^.A 

Stress  constraints  only  are  imposed.  The  S.O.C., which  reduces 
here  simply  to  F.S.D., leads  to  instability  in  convergence.  The 
H.O.C.  leads  Co  the  linearization  of  58  out  of  the  126  constraints 
and  is  therefore  close  to  the  G.O.C.  It  converges  quite  rapidly 
to  a design  that  has  also  been  found  by  SHMIT  [ S18  ] . Figure 
8.24  illustrates  the  convergence , which  is  detailed  in  table  8.23. 
The  final  design  (4975.1  lbs)  is  presented  in  table  8.26. 

It  is  evident  that,  for  this  problem,  the  linearization  of  the 
stress  constraints  is  an  essential  ingredient  for  fast  convergence. 
The  method  of  SCHMIT  [ S18  ] is,  in  this  case,  very  close  Co  the 

G.O.C.  as  it  solves  almost  exactly  the  problem  L generated  after 
each  analysis. 

Prob2_em 

The  relative  displacement  at  nodes  1 and  2 is  limited  to  1.0  in. 
Again  the  S.O.C,  using  F.S.D.  for  the  stress  constraints,  exhibits 
poor  and  unstable  convergence  (6159  lbs  in  50  iterations)  while 
Che  G.O.C.  converges  very  rapidly  (5120  lbs  in  12  iterations) 
just  as  Che  SCHMIT  [ S18  ] algorithm  (5121  lbs  in  13  iterations) 

which  is  closely  related.  The  results  are  displayed  in  figure  8.27 
and  in  tables  8.25  and  8.26.  The  performances  of  Che  G.O.C.  is 
obtained  at  the  expense  of  a large  number  (126)  of  virtual  loading 
cases  corresponding  to  each  linearized  stress  constraints.  An 
investigation  has  been  done  using  the  automatic  selection  option 
of  the  program  OFTBAR  for  the  stress  constraints  chat  are  most 
critical,  as  based  on  Che  criterion  (3.38).  The  tolerance  for  the 
selection  has  been  made  more  and  more  severe,  corresponding  to 
0,  6,  7,  10  and  27  linearized  stress  constraints  so  Chat  a progres- 
sive transition  from  S.O.C.  to  G.O.C.  is  achieved. 
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The  results  are  illustrated  in  figure  8.28  as  a function  of  the 
number  of  reanalyses  and  in  figure  8.29  as  a function  of  the  total 
CPU  time.  It  appears  that  the  economic  optimum  is  obtained  when 
Che  10  most  critical  stress  constraints  are  linearized.  This  inves- 
tigation indicates  that  the  generalized  use  of  the  selection  crite- 
rion (3.38)  for  Che  stress  constraints  is  strongly  recommended. 
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8.1.5.  Four  level  72  bar  space  truBs 

This  problem  has  been  investigated  by  many  authors 
I T3,  G13,  V3,  B8,  S8,  S18  ] . It  is  defined  in  figure  8.30 

and  in  table  8.31.  Two  loading  cases  and  four  flexibility  cons- 
traints are  specified.  The  stress  constraints  arc  uniform  at 
o - 25000  psi  but  reveal  not  very  critical.  The  optimum  is  obtained 
in  5 iterations  by  the  S.O.C.  and  is  identical  to  the  solution 
obtained  by  TAIG  ( T3  ],  which  is  quite  understandable  as  in  both 
cases  the  same  explicit  problem  A is  solved  exactly  at  each  itera- 
tion. The  convergence  is  illustrated  in  figure  8.32  and  in  table 
8.33.  The  final  design  is  given  in  table  8.34.  The  results  obtained 
Lj  by  various  authors  are  also  presented.  It  appears  clearly  that  this 

example  is  well  adapted  to  the  F.S.D.  approximation  of  the  stress 
constraints.  Any  additional  sophistication  increases  the  cost  of 
computation  without  improving  the  global  efficiency. 
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8.1.6.  Transmission  Cowar 


This  example  is  very  popular  in  Che  scruccural  opcimizacion 
licCeraCure  ( F27,  G12,  M9 , V2,  D5  ] . IC  is  defined  by  Che  figure 

8.33  and  Che  Cable  8.36.  Using  Che  symmeCry  of  Che  problem  25 
design  variables  are  necessary.  Two  loading  cases  and  Cwo  flexibi- 
licy  consCrainCs  are  specified.  The  sCress  limiCaCions  are  given 
differenc  values  in  craccion  and  compression.  The  scress  limics 
used  by  GELLATLY  ( G13  ) and  SCHMIT  ( S8,  S18  ] have  been  used. 

They  differ  slighcly  from  Chose  mencionned  by  VENKAYYA  [ V3  ] 
buC  Chis  difference  does  noC  seem  Co  influence  Che  resulcs. 

In  figure  8.37  and  Cable  8.38  are  compared  Che  performances 
obCained  by  using  Che  S.O.C.  and  Che  H.O.C.  In  Che  H.O.C.  Che  scress 
conscrainc  in  elemenC  20  for  Che  second  loading  case,  is  Che  only 
one  linearized.  Ic  leads  Co  a moderaCe  improvemenC.  All  Che  algo- 
richms  proposed  in  Che  liCCeraCure  yield  almosC  idenCical  designs, 
which  are  given  in  Cable  8.39. 


125. 


8.1.7.  Simplified  airbrake  structure 


This  problem  reproduces  on  a very  simple  equivalent  truss, 
the  essential  difficulties  met  on  a more  complex  airbrake  struc- 
ture which  is  examined  in  section  8.3.  The  problem  is  reduced 
to  the  hyperstatic  57  bar  truss  illustrated  in  figure  8.40  and 
precised  in  table  8.41.  Only  one  loading  case  is  applied. 

The  difficulty  comes  from  the  flexibility  constraints  which  impose 
that  the  difference  of  the  w displacements,  along  the  z axis, 
between  three  couples  of  nodes  1,  2,  3 should  be  less  than  0.5  mm. 
This  is  equivalent  to  specifying  that  the  trailing  edge  of  the 
airbrake  has  to  remain  straight.  A maximum  stress  is  imposed  as 
well  as  a minimum  cross  section. 


All  the  optimality  criteria  failed  on  this  problem. 

The  solution  was  obtained  only  by  using  stress  constraint  lineari- 
zation and  the  mixed  method  with  limitation  of  the  number  k of 
steps.  The  results  are  illustrated  by  the  figures  8.43  and  8.44 
(note  the  difference  in  the  scale  for  the  weight)  and  by  the  tables 
8.45  and  8.46. 


Consider  first,  in  figure  8.43,  the  results  obtained  by 
using  the  F.S.D.  zero  order  approximation  for  the  stress  constraints. 
As  an  example  of  C.O.C.,  the  envelope  method  [ G13  ] is  used  and 

diverges  immediatdy,  yielding  designs  that  remain  more  than  10 
times  heavier  than  the  best  found  subsequently. 

The  S.O.C.  is  used  via  the  primal  (k  •■  •>  ) and  the  dual  methods. 

It  reveals  unstable, but  converges  after  20  analyses  to  a design 
weighting  20.5  kg,  which  is  still  far  from  the  best  design  (15.6 
kg).  Using  the  same  F.S.D.  approximation  for  the  stress  constraints 
but  limiting  the  number  of  steps  k improves  the  speed  and  stability 
of  the  convergence,  but  still  does  converge  to  the  same  non  optimum 
design.  The  number  k was  limited  to  57  using  the  first  order  primal 
algorithm.  The  results  are  denoted  SOC-MMl  (k  ■ 57).  Using  the 
first  order  algorithm,  it  seems  that  a good  estimation  of  the 
number  of  steps  k is  given  by  the  number  of  elements,  which  is  57 


here . 
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Consider  now,  in  figure  8.4A,  the  results  obtained  with 
linearization  of  the  stress  constraints.  If  all  of  them  are 
linearized,  the  exact  solution  of  the  corresponding  explicit 
approximate  problem  corresponds  to  a G.O.C.  The  convergence 
of  a sequence  of  such  problems  is  still  non  monotonic , but  yields 
in  20  analyses  a design  weighting  15.6  kg  which  is  the  minimum 
found  for  this  problem.  During  the  next  10  to  13  iterations  the 
results  oscillate  slightly  around  this  design.  Linearization 
of  all  the  constraints  being  costly,  the  H.O.C.  has  been  inves- 
tigated using  the  autoratic  selection  algorithm  of  the  program 
OPTBAR  (c^  ■ 0.8,  C2  * 0.02).  Indeed  the  algorithm  leads  to 
linearize  only  one  of  the  five  stress  constraints  which  are 
active  at  the  optimum. 

The  corresponding  H.O.C.  is  therefore  cheap  to  apply. 

If  k is  unlimited  however,  the  behavior  is  similar  to  using 
F.S.D.  as  in  figure  8.43  (S.O.C.).  If  the  mixed  method  is  applied 
(H.O.C.-MMl,  k ■ 57)  the  algorithm  converges  very  rapidly  to  a 
near  optimum  design  (17.32  kg)  in  7 analyses.  During  the  next 
10  to  15  analyses  the  design  is  only  slightly  udapted  and  the 
weight  reduces  to  15.57  kg.  After  the  7^^  reanalyses,  Che  limi- 
tation of  k has  been  suppressed  (H.O.C.-MMl,  k - after  7^^ 
step).  This  change  reintroduces  two  large  oscillations  of  Che 
weight  but  leads  to  an  identical  final  design  weighting  15.57  kg. 

In  this  example,  Che  15.6  kg  design  was  obtained  only  with 
stress  linearization.  It  illustrates  the  need  for  an  automatic 
algorithm  for  detecting  Che  critical  stress  constraints.  If  all 
the  constraints  are  linearized  (G.O.C.)  Che  optimality  criterion 
converges  to  Che  lighter  design  but  is  very  expensive.  The  combi- 
nation of  one  stress  constraint  linearization  and  of  Che  mixed 
method  yields  a faster  convergence  at  much  cheaper  computing  cost, 
as  shown  in  the  cable  8.47  which  summarized  the  computing  times. 

Note  also  that  the  2 designs  obtained  by  S.O.C.  and  by  H.O.C. 
(or  G.O.C.)  are  compared  in  figure  8.46.  Although  the  weights  are 
quite  different  (20.5  kg  and  15.6  kg)  as  well  as  the  design  varia- 
bles themselves,  Che  same  constraints  are  active.  They  are  the 
two  flexibility  constraints  and  5 stress  constraints. 
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8.1.8.  Computing  times 

The  execution  times  for  most  of  the  truss  problems  inves- 
tigated in  this  chapter,  are  recorded  in  table  8.47. 

The  methods  and  the  program  used  are  indicated.  It  is  recalled 
that  OPTBAR  is  a simple,  experimental,  all  in  core  program  used 
for  research  essentially,  while  SAMCEF  is  a general  purpose  pro- 
gram with  much  larger  analysis  capability,  but  wliich  is  evidently 
less  efficient  in  small  scale,  academic  applications. 


The  execution  times  are  split  in  preprocessing,  which  exists 
only  for  the  SAMCEF  program,  in  analy s is , includ i ng  eventually  the 
computation  of  the  gradients,  and  in  minimization  time.  The  time  is 


^■total 


t 

preproc . 


+ n , 


Iter 


(t 


analysis 


+ 


t . . 

minim . 


) . 


This  splitting  allows  to  appreciate  the  increase  of  computing 
time  due  to  the  use  of  more  sophisticated  minimization  methods. 

In  some  cases,  especially  with  the  G.O.C.  where  all  the  constraints 
are  linearized,  this  increase  in  minimization  time  is  clearly  pro- 
hibitive, as  it  is  not  counterbalanced  by  a sufficient  reduction 
of  the  number  of  iterations. 


The  comparison  of  the  computing  times  obtained  for  the 
various  methods  reinforces  the  conclusion  drawn  from  the  examina- 
tion of  the  convergence  properties,  which  is,  that  the  H.O.C. 
appears  as  the  best  compromise, especially  if  the  selection  of  the 
critical  stress  constraints  is  achieved  automatically  and  if  the 
mixed  method  limiting  k is  used. 
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8.2.  Box  beam  structures 


8.2.1.  Simple  wing  box 

The  simple  box  beam  illustrated  in  figure  8.48  is  a clas- 
sical test  case  proposed  in  the  litterature  ( G12,  G13,  S18  1 . 

A detailed  statement  of  the  problem  is  given  in  table  8.49. 

Two  load  cases  are  considered.  Maximum  allowable  stresses  are 
given  for  the  bars  and  membranes  elements  which  are  grouped  as 
indicated  in  table  8.49^ so  that  only  16  design  variables  are  invol- 
ved. Minimum  values  are  also  specified. 

By  symmetry,  only  one  half  of  the  structure  is  modeled  in 
finite  elements.  Simple,  constant  strain,  conforming  displacement 
elements  have  been  used,  yielding  to  45  D.O.F.  In  the  membrane 
elements,  the  von  MISES  equivalent  stress  is  limited.  The  lineari- 
zation of  these  constraints  has  been  achieved  according  to  the 
method  presented  in  (3.35). 

The  S.O.C  has  been  applied.  It  reduces  here  to  the 
straightforward  F.S.D.  The  convergence  is  rather  fast  and  monotonic, 
as  illustrated  in  figure  8.50  and  in  tables  8.51  and  8.52. 


When  the  G.O.C  is  applied,  one  realizes  that  the  optimum 
is  not  a F.S.D  and  that,  with  the  linearized  stress  constraints, 
the  convergence  is  much  faster  and  leads  to  a significantly  lighter 
design.  Comparing  the  results  in  table  8.52  indicates  that  the  two 
designs  differ  also  significantly. 


The  computing  time  for  this  problem,  using  the  SAMCEF  pro- 
gram, are  as  follows 


S.O.C  20  iterations  preprocessing 

analysis 

minimization 

total 


12.0 

8.8 

0.5 

198  sec. 


G.O.C  10  iterations  preprocessing 

analysis 


14.1 

36.9 
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minimization  0.8 

total  390  sec. 

It  indicates  again  that  the  large  number  of  virtual  load  cases 
and  the  additional  computing  cost  induced  by  the  G.O.C  lead  to 
almost  double  the  computing  time  although  the  number  of  iterations 
is  just  half. 

The  H.O.C  should  bring  a more  efficient  way  of  reaching  the 
200  lbs  design,  but  it  was  not  used  here.  The  reason  is  that  the 
comparison  of  the  results  with  those  presented  in  the  litterature 
is  rather  difficalt  as  the  finite  element  models  used  by  the 
various  authors  are  not  the  same.  In  such  a simple  model  the  diffe- 
rences due  to  using  one  element  or  another  are  quite  important  and 
make  the  comparisons  non  significant. 


8.2.2.  A composite  box  beam 


The  simple  box  beam  illusCraCed  in  figure  8.53  and  detailed 
in  table  8.54  has  been  proposed  by  KHOT,  VENKAYYA  and  BERKE  | K15  ] 

as  an  example  of  composite  material  design.  The  upper  and  lower 
skin  elements  are  supposed  to  be  made  of  composite  material  (boron- 
epoxy  as  defined  in  table  8.54)  while  the  transverse  members  and 
the  reinforcing  bars  are  in  aluminium.  The  transverse  members  are 
idealized  by  shear  panels,  while  the  skin  elements  are  obtained  by 
superposition  of  4 constant  strain  displacement  finite  elements. 
Each  of  them  corresponds  to  a different  orientation  of  the  fibers, 
that  is  0*,  90*,  ♦ 45*  and  - 45*.  The  superposition  of  these  4 
layers  yields  a unique  definition  of  the  design.  The  material  pro- 
perties in  each  layer  are  orthotropic  and  correspond  to  a unidirec- 
tional reinforcement  of  the  resin.  The  tickness  of  each  layer  is 
a design  variable. 

In  the  design  of  composite  structure,  a frequent  problem 
is  to  limit  the  displacement  or  to  tailor  the  flexibilities  accor- 
ding to  given  laws.  Consequently,  no  stress  limitations  have  been 
introduced  in  this  problem  which  is  governed  entirely  by  two  diffe- 
rent flexibility  constraints. 

I.imi£a£i£n_in  ^end^n^  (Problem  2. 2. A) 

In  this  case  4 tip  loads  are  applied  symmetrically  at  nodes  1,  2, 

3,  4 and  the  displacements  at  each  of  the  36  nodes  are  limited 
(table  8.54).  The  double  symmetry  of  the  problem  has  been  used. 

It  allows  to  reduce  the  number  of  design  variables  to  45  and  the 
number  of  constraints  to  9. 

The  S.O.C  has  been  used  and,  as  the  number  of  flexibility  cons- 
traints is  rather  large,  the  primal  method  has  been  selected  with 
the  second  order  (NEWTON)  minimization  algorithm.  The  convergence 
is  illustrated  in  figure  8.55  and  tabulated  in  table  8.56  where  it 
is  compared  with  the  results  obtained  by  KHOT  [ K15  ] . The  slight 

differences  in  the  final  designs,  presented  in  table  8.57,  are  due 
to  the  differences  in  the  finite  element  models. 
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Limi£a£i£n_in  £0£8ion  (Problem  2.2.B) 

The  loading  is  now  asymmetric  and  induces  primarily  bending 
deformations  but  also  a limited  amount  of  torsion.  The  limita- 
tions in  deflections  are  14  inches  at  nodes  1 and  2 where  the 
loads  are  1000  lbs  and  of  IS  inches  at  nodes  3 and  4 where  the 
loads  are  975  lbs.  These  limitations  require  evidently  a non 
symmetric  design.  As  the  tip  deflections  arc  the  only  one  limited, 
the  shear  webs  keep  the  same  thickness  along  the  span.  This  allows 
to  reduce  the  number  of  design  variables  to  30. 

The  S.O.C  has  been  used  and  the  results  are  presented  in 
figure  8.55  and  in  tables  8.58  and  8.59  together  with  those 
obtained  by  KHOT  [ K15  ] . 

In  both  problems  (2. A and  2.B)  the  convergence  of  the 
process  is  monotonic  and  significantly  faster  than  reported  by 
KHOT  ( K15  ] . The  computing  times  are  the  following,  using  the 

SAMCEF  program, still  on  IBM  370/158. 

problem  2. A 10  iterations  preprocessing  34.3 

analysis  57.3 

minimization  1.1 

total  617  sec. 

problem  2.B  10  iterations  preprocessing  34.9 

analysis  27.4 

minimization  0.9 

total  317  sec . 
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8.3.  Airbrake  structure 

The  airbrake  structure  illustrated  in  figure  8.60  has  been 
used  in  an  exercise  to  determine  the  capability  of  the  methods 
proposed  in  the  SAMCEF  program,  when  applied  to  a problem  with  a 
number  of  design  variables  and  of  degrees  of  freedom  larger  than 
those  met  in  the  classical  tests  extracted  from  the  litterature 
on  optimization.  Moreover  the  flexibility  constraints  revealed 
rather  difficult  to  take  into  account  in  the  classical  algorithms, 
although  they  could  be  frequently  encountered  in  practical  problems. 
Only  the  essential  results  obtained  for  this  problem  are  presented 
here.  A more  detailed  description  has  been  given  by  FLEURY  in 
( F18,  F29  ] . 

The  structure  illustrated  in  figure  8.60  is  classically 
designed  in  light  alloy  aluminium  sheet.  The  main  (front)  spar 
and  the  secondary  spar  are  joined  by  twelve  ribs  and  covered  by 
two  skins  reinforced  by  stringers.  The  airbrake  is  hinged  at  three 
points  and  actuated  at  one,  in  the  midspan. 

The  loads  consist  in  pressure  distributions  on  both  faces, 
corresponding  to  two  flight  configurations.  In  one  of  them,  a 
flexibility  constraint  is  imposed  which  stipulates  that  the  trai- 
ling edge  has  to  remain  straight  within  a tolerance  e ■■  0.5  mm 
to  insure  aerodynamic  efficiency. 

Stress  constraints  are  imposed  and  take  different  values  depending 
of  the  material  used  and  of  other  technological  considerations. 

For  membrane  elements,  the  von  Mises  equivalent  stress  is  used. 
Minimum  thicknesses  are  assigned  which  also  differ  from  places  to 
places  depending  of  manufacturing  considerations. 

Re^u£t^orj  £^_t]ie_n£mb£r__o^  ^ ^ 

In  a finite  element  model,  as  illustrated  in  figure  8.61,  the 
flexibility  constraint,  as  initially  stated,  has  to  be  written 


1 w.  - Wj^  1 < e j - 1,  N-1  (8.  XI) 
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were  w.,  w.  are  the  transverse  deflections  at  any  node  j or  k 
J 

along  the  trailing  edge.  If  N is  the  number  of  nodes  along  the 
trailing  edge  in  the  finite  element  model,  the  total  number  of 
constraints  is 

N (N  - 1)  / 2 

which  is  prohibitive,  when  it  is  considered  that  10  to  12  nodes 
are  to  be  used  along  the  trailing  edge.  After  having  solved  the 
problem  using  a simplified  model  where  N - 5,  it  was  recognized 
and  confirmed  by  physical  intuition  that  the  constraints  (8.3.1) 
could  be  replaced  by  the  reduced  set 


where  M and  H*  are  two  nodes  selected  so  that  the  distances  (1-M) 
and  (M*-N)  are  approximately  the  same.  Eventually  a certain  number 
of  couples  of  nodes  (M-M*)  could  be  selected.  It  was  verified  a 
posteriori  that,  using  the  reduced  set  of  constraints,  all  the 
initial  constraints  (8.3.1)  are  satisfied. 


£ijrs£  simplified  model  (2  7 elements) 

In  order  to  determine  the  behavior  of  the  algorithms  in  the 
presence  of  the  flexibility  constraints  (8.3.1),  a very  simple 
model  illustrated  in  figure  8.62  was  built  up.  It  consists  in  26 
membrane  elements,  plus  one  bar  used  to  model  the  control,  which 
leads  to  43  D.O.F.  Five  nodes  are  taken  along  the  trailing  edge 
and  the  10  exact  flexibility  constraints  (8.3.1)  are  imposed, 
but  no  stress  constraint. 

This  problem  reveals  the  obvious  need  to  limit  the  global 
deflection  of  the  trailing  edge. 

w . < w 
J 


T 


1 j'* . 


Indeed,  if  one  consider  a design  yielding  a perfectly  straight 
deflected  trailing  edge,  the  application  of  scaling  factor  does 
not  affect  the  satisfaction  of  the  constraints  (8.3.1). 

Thus  the  process  converges  to  a zero  weight  design  with  a very 
larger  deflection.  This  happens  because  the  stress  constraint 
are  not  taken  into  account  in  this  simplified  problem.  However 
it  was  felt  inconsistent  to  leave  the  limitation  of  the  global 
deflection  to  the  stress  constraints,  as  it  is  not  known  a 
priori  whether  they  become  active  before  reaching  an  unacceptable 
large  deflection.  Therefore  a set  of  five  constraints  (8.3.3) 
was  added  to  the  problem,  with  w - 10  mm. 


The  envelope  method  [ G13  ] has  been  applied  as  an  example 

of  C.O.C.  Next  the  S.O.C  has  been  used,  which  is  here  equivalent 
to  a C.O.C  as  all  the  constraints  are  flexibility  constraints, 
and  thus  are  exactly  linearized.  The  results  are  presented  in 
figure  8.63  and  revealed  extremely  poor  (if  any)  convergence. 


This  implies  that  the  constraints  (8.3.1)  are  very  non  linear,  | 

1 

even  in  the  inverse  design  space.  Using  the  mixed  method  (MMl), 
even  with  k being  small,  did  not  improve  significantly  the  conver- 
gence. Looking  more  carefully  at  the  optimization  process,  in  the 
mixed  method,  it  appears  that  the  difficulty  arises  from  the  fact 
that,  at  each  iteration,  the  algorithm  selects  a design  where  the 
maximum  deflection  constraint  (8.3.3)  is  active.  After  reanalysis,  j 

the  differential  constraints  (8.3.1)  are  seriously  violated  and 
the  subsequent  seal ing  , which  has  to  satisfy  them, leads  to  increase  j 

the  weight  in  such  way  that  the  divergence  starts  (see  the  comments  I 

of  figure  7.5).  i 


It  suggests  the  introduction  of  a kind  of  move  limit  in  the 
form  of  a limitation  to  the  increase  of  global  deflection  between 
two  iterations. 


with 


Wj  < min  (w*  , w) 


w*  ■ 6 max  (w.) 

j-l,N  J 
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where  6 is  a constant  close  to  unity.  In  this  way  the  deflection 
cannot  exceed  6 times  the  value  obtained  after  the  current  ana- 
lysis and  scaling  cycle.  The  constant  d defines  the  amplitude 
of  the  "move"  that  is  allowed  in  the  deflection. 


The  move  limit  parameter  6 and  the  convergence  control  para- 
meter k can  now  be  use  to  improve  the  convergence  behavior  of 
the  problem.  Taking  k > 20  in  the  MMl  for  20  iterations  and  then 
k ■ 10,  and  allowing  6 - 1.1  brings  a nice  convergence  as  indicated 
in  figure  8.64,  which  also  contains  the  history  of  the  global  de- 
flection,and  that  of  the  trailing  edge  deflection.  Note  that  40 
iterations  are  necessary  to  stabilize  the  weight.  The  change  of  k 
at  iteration  21  was  introduced  after  a beginning  of  divergence 
was  detected.  This  weight  should  not  be  compared  to  that  of  the 
final  model  since  no  attention  was  payed  to  assigning  minimum 
values  for  the  design  variables  which  are  consistent  with  the 
full  scale  problem. 

The  investigation  on  this  first  simplified  model  revealed 
the  need  of  a careful  control  of  both  the  number  k of  steps  in 
the  mixed  method  and  of  the  move  limit  5.  The  optimality  criteria, 
even  in  the  sophisticated  form  of  the  G.O.C,  failed  completely. 

It  also  revealed  that  the  set  of  constraints  (8.3.1)  can  be  repla- 
ced by  the  reduced  set  (8.3.2)  with  only  one  couple  of  points 
(M.M*) . 

^e£ond_s_im£l^f_ie^  (i25  elements) 

A second,  still  simple  model  was  used,  which  is  illustrated 
in  figure  8.65.  42  membrane  elements  and  83  bars  are  used  leading 
to  103  D.O.F. 


The  flexibility  constraints  are  written 
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I I < 0.5  mm  I | < 0.5  mm 

1 Wj  ~ Wg  I < 0 . 5 mm  | | < 0 . 5 mm 

Wj  < 10  mr  Wg  < 10  mm 

The  stress  constraints  and  the  minimum  sizes  of  the  full  scale 
problem  are  taken  into  account.  The  stress  constraints  are 
treated  by  F.S.D. 

The  mixed  method  (MMl)  has  been  applied  with  k * AO  up  to 
the  sixth  iteration  and  unlimited  after.  The  move  limit  on  the 
global  deflection  was  still  defined  by  6 > 1.1.  The  initial 
design  was  obtained  from  the  first  simplified  model.  The  conver- 
gence curves  are  illustrated  in  figure  8.66. 

Note  that,  due  to  the  introduction  of  the  stress  constraints 
the  deflection  does  not  reach  the  maximum  value  of  10  mm  but  flat- 
tens at  3.6  mm.  However  the  move  limit  is  still  necessary  to  avoid 
the  divergence.  Although  the  convergence  is  rather  fast  (9  itera- 
tions) it  is  felt  that  even  faster  convergence  could  be  obtained 
by  a larger  value  of  k,  around  100.  Note  that  the  second  order 
algorithm  (MM2)  was  not  available  at  the  time  of  the  investigation. 

It  should  also  be  pointed  out  that  the  weight  cannot  be 
compared  with  those  obtained  by  the  other  models  as  no  attention 
was  payed  to  the  exact  equivalence  in  the  statement  of  the  problem 
and,  more  precisely,  of  the  stress  and  minimum  size  constraints. 

T^he^  ^ina^  (627  elements) 

The  final  model  is  illustrated  in  figure  8.67. 

It  is  made  of  378  membrane  elements  and  249  bars.  Second  degree 
displacement  fields  are  used,  leading  to  2300  D.O.F.  The  sizes 
of  all  the  elements  have  been  left  free,  yielding  627  design 
variables.  In  practice  a certain  number  of  elements  should  have 
been  given  identical  sizes.  This  is  achieved  in  the  SAMCEF 


program  by  defining  groups  of  elements, which  are  assigned  Che 
same  design  variables.  It  has  Che  advantage  of  reducing  the 
number  of  design  variables  and  of  yielding  more  realistic  designs. 
However,  the  present  problem  being  essentially  an  exercise,  it 
was  decided  to  solve  it  with  the  maximum  number  of  design  variables 
The  results  obtained  by  grouping  the  elements  according  to  various 
design  options  are  reported  in  [ F29  ] . 

The  627  stress  constraints  are  treated  by  F.S.D  while  the 
flexibility  constraints  are  written 


»!  - We  1 < 0.5 


Wj^2  " Wg  1 < 0.5  mm 


< 10  mm 


Wi-w?  I <0.5  mm 


Wj^2  " w^  I < 0.5  mm 


Wj^2  ^ iO  mm 


As  expected  from  Che  experience  gained  on  the  simplified  pro- 
blems, a good  convergence  was  obtained  with  the  mixed  method  MMl , 
using  the  conjugate  gradient  algorithm,  by  setting  6 ■ 1.1  and 
k.  « 500,  that  is  slightly  below  the  number  of  design  variables. 
The  choice  of  k was  guided  by  the  fact  Chat  the  convergence  of 
first  order  minimization  algorithms  is  achieved  in  a multiple 
of  Che  number  of  unknowns  and  hence  the  same  ratio  between  k 
and  Che  number  of  design  variables  was  kept  as  in  the  previous 
problems . 

The  evolution  of  the  weight  is  illustrated  by  Che  figure 
8.68.  Again  Che  convergence  is  very  satisfactory  in  10  to  12 
iterations.  The  maximum  deflection  of  Che  trailing  edge  is  not 
reached.  The  satisfaction  of  the  exact  constraints  (8.3.1)  can 
be  verified  a posteriori.  The  initial  scale  up  of  the  weight  is 
due  Co  the  fact  chat  the  original  design  of  the  airbrake  did 
not  satisfy  the  flexibility  constraints. 


1 


The  computing  times  for  this  final  application  are, 
IBM  370/158 


13  iterations 


preprocessing 

282 

analysis 

813 

minimization 

196 

total 

13  203 

p ■ O.I  Ib/in 

0.  ■ 25000  p«i 

^ 2 
a . ■ 0.1  in 
-1 

initial  valuaa 
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10  BAR  TRUSS  - PROBLEM  l.D 


FINAL  DESIGN 


element  section  element  section 


21.579 

3.434 

0.100 

1.401 

22.663 


WEIGHT  CONVERGENCE 


section 

0 

.100 

14 

.545 

12 

.695 

1 

.932 

11 

.928 

iteration 


Weight 

(lbs) 


7963.06 

5141.56 
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4043.96 
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FIGURE 


DOUBLE  CANTILEVER  22  BAR  TRUSS 


GENERAL  DATA 


E - lO'  psl 
a - 0.001  In^ 

a°*  1.0  in  (InltlAl  design) 


Loading  : node  13 « component  Y,  - 100  000  lbs 


ALLOWABLE  STRESS  LIMITS  AND  MASS  DENSITY 
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0 (psl) 

1 

2A 

1-20 

25000 

21,22 

500000 

2B 

1-10 

5000') 

11-20 

25000 

21,22 

5000)') 

2C 

1-10 

50000 

11-20 

25')0') 

21,22 

5000)0 

p (lb /111  ) 
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DOUBLE  CANTILEVER  22  BAR  TRUSS  - PROBLEM  2. A 
FINAL  DESIGN 


FSD 

GOC 

elements 

secti^on 
( In") 

Stress 

(psl) 

section 

(In^) 

Stress 

(psl) 

1,11 

2.00069 

- 25000 

2.01023 

- 25000 

2,12 

1.09937 

- 25000 

1.000.38 

- 24900 

3,13 

0.00100 

153 13 

). 01101 

20733 

4,14 

0.00100 

1581R 

0.01101 

20783 

5,15 

3.09032 

25000 

3.90004 

24993 

6,16 

0.00100 

1374 

0.0»)100 

3131 

7,17 

2.32745 

- 25000 

2.81506 

- 24990 

3,13 

0,00100 

24313 

0.01731 

20.106 

9,10 

0.0010(3 

- 22370 

0.01634 

- 207 3M 

10,20 

2.32753 

25000 

2.31553 

24000 

21,22 

0.10000 

500000 

0.10002 

40^003 

Weight  (Ibs^ 

1224 

.135 

1226 

.031 

number  of 
analyses 

50 

30 
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DOUBLE  CANTILEVER  22  BAR  TRUSS  - PROBLEMS  2A.  2B.  2C 
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DOUBLE  CANTILEVER  22  BAR  TRUSS  - PROBLEM  2.B 
FINAL  DESIGN 


elements 

FSD 

HOC 

section 

(in^) 

stress 

(psl) 

section 

(in^) 

. _ 

stress 

(psl) 

1 

0.90030 

- 46524 

1. 99009 

- 5o;)oo 

2 

0.72350 

- 48127 

1.03037 

- 5(V''’0 

J 

0.00207 

40234 

0.t)0li)0 

31661 

4 

0.00207 

40234 

0.0010) 

31661 

5 

1.53960 

47654 

3.07031 

500;)0 

6 

0.00103 

- 16756 

0.00100 

2375 

7 

1.02267 

- ^3127 

2.31331 

- 50<)00 

3 

0.13100 

40314 

0.00100 

43563 

9 

0.13148 

- 40284 

0.0010) 

- 44776 

10 

1.02313 

43127 

2.31339 

50000 

11 

2.73063 

- 24167 

0.02060 

- 250)0 

12 

2.21350 

- 25000 

0.01937 

- 250  )0 

13 

0.23533 

20930 

0.00100 

15313 

14 

0.23533 

20930 

0.00 100 

15313 

15 

4.72201 

24754 

0.03931 

25000 

16 

0.00103 

3720 

0.00100 

1374 

17 

3.13090 

- 25000 

0.02731 

- 25000 

13 

0.40411 

20934 

0.00100 

24313 

10 

0,40353 

- 20930 

0.00100 

- 22370 

20 

3.13744 

25000 

0.02739 

25-*O0 

21 

0.03165 

472340 

0.2)729 

43  )0()  ) 

22 

0.12519 

4907  P 

O.OOlOO 

50  )000 

Weight  (lbs) 

1044 

.42 

654.049  (nnn  FSD) 

number  of 

t 

.7  i 

analyses. 

i 

153 


TAliLli  3 . 15 

DOUBLE  CANTILEVER  22  BAR  TRUSS  - PROBLEM  2.C  - FINAL  DESIGN 


FSU 

HOC 

olements 

section 

(In) 

stress 

(psi) 

section 

(in^) 

stress 

(psl) 

1 

0.27941 

- 48151 

0.01029 

- 50f)0() 

2 

0.27813 

- 48144 

0.00397 

- 50iK)0 

3 

o.ooiol 

30461 

0.00100 

31035 

4 

0.00101 

30461 

0.00100 

31635 

5 

0.55632 

43146 

0.01851 

50000 

6 

o.ooiol 

2645 

0.00100 

2743 

7 

0.39326 

- 43144 

0.01260 

- 50)0') 

3 

o.ooiol 

46819 

o.'vnoo 

43626 

0 

0.00101 

- 43078 

0.00100 

- 44739 

10 

0.39333 

43144 

0.01263 

50)00 

11 

3.46404 

- 24999 

3.93149 

- 25'))0 

12 

3.46255 

- 25000 

3.93017 

- 25‘)Oi) 

13 

o.ooiol 

16175 

0.00100 

16477 

14 

0.00101 

16175 

0.00100 

16477 

15 

6.92560 

25000 

7.96091 

25000 

16 

0.00101 

3 1'W 

0.00100 

4553 

17 

4.89656 

- 25000 

5.6  2373 

- 25')0O 

18 

0.00117 

23451 

0.00100 

24313 

19 

0.00101 

- 22875 

0.00 100 

- 22370 

20 

4.39679 

25000 

5.62331 

25  000 

21 

0.02787 

481476 

O.fVUOO 

480  ),)!) 

22 

0.17316 

499091 

0.199f)4 

50)0)0 

Weight  (lbs) 

1391 



.37 

1230.0) 

(I'SD) 

nu'ober  of 
analyses 

50 

27 

TABLE  3.17 


200  BAR  PLANAR  TRUSS 


General  data  : 


E ■ 3 10  psi 

p - 0.233  Ib/ln^ 

o - lOOOn  psl  (in  problem  3.B) 

^ « 0,1  in^ 

a = 10.0  in  (initial  desif^n) 


Multiple  load  ins  caae 


Load 

cases 


nodes 


load  components  (lbs) 


1,6,15,20,29,34,43,43,57,62,71 


1 , 2 , 3 , 4 , 5 , 

6 , 3 , 10  , 12  , 14  , 


15  , 

16  , 

17  . 

13  , 

19 

20  , 

22  , 

24  , 

26  , 

28 

• • • 

71  . 

72  , 

73  , 

74  , 

75 

load  cases  1 plus 


Flexibility  constraints  (max.  values) 


components  (in) 


- 10010 


N*”"  ANALYSES 

-j I j — j I — ^ 

4 6 8 10  12 


BAR  PLANAR  TRUSS 


PROBLEM  3. A 


13B 


1.4 

2.3 

5,17 

6,16 

7.15 

d,14 

'7.13 

10,12 

11 

13,25 

10,24 

20,23 

21,22 

26.33 
n,'il 
23,36 
29,35 

30.34 
31,33 
32 

39,42 

40,41 

43,55 

44,54 

45,53 

46,52 

47,51 

43,50 

49 

56,63 

57,62 

53,61 

59,60 

64,76 

65,75 


.147 
.100 
.333 
.230 
. 100 
.241 
.100 
.103 
.777 
O.lfX) 
0.100 
0.100 
0.100 
6.363 
0.100 
0.313 
3.137 
0.133 
0.144 
3.931 
0.332 
0.210 
7.960 
0.640 
0.100 
4.323 
0.291 
0.234 
4.754 
0.100 
0.100 
0.100 
O.lfK) 
9.239 
0.100 


0.100 

2.021 

0.100 

0.115 

2.795 

0.100 

0.100 

0.100 

0.100 

6.965 

0.100 

0.437 

2.371 

0.133 

0.125 

4.009 

0.573 

0.256 

7.806 

1.015 

0.100 

4.244 

0.308 

0.210 

4.749 

0.100 

0.111 

0.111 

0.147 

9.124 

oa22 


scctloi 

( in") 

case  A 

case  B 

66,74 

0.741 

1.139 

133,133 

67,73 

4.975 

4.370 

134,137 

63,72 

0.251 

0.216 

115,136 

69,71 

0.316 

0.327 

140,  152 

70 

5.525 

5.520 

141,151 

77,30 

0.341 

1.019 

142, 150 

78,79 

0.373 

0.242 

141,149 

31,93 

9.273 

8.94) 

144,141 

32,92 

1.441 

1.735 

145,147 

83,91 

0.100 

0.100 

146 

84,90 

6.527 

6.327 

153,150 

85,39 

0.443 

0.441 

154,155 

36,83 

0.403 

0.367 

157,169 

37 

6.009 

5.945 

153,163 

94,101 

0.100 

0.100 

159,167 

95,100 

0.150 

0.196 

160,166 

96,99 

0.150 

0.196 

161,165 

97,98 

0.205 

0.227 

162,164 

102,114 

10.266 

9.954 

163 

103,113 

0.169 

0.203 

170,177 

104,112 

1.586 

1.913 

171,176 

105,111 

7.097 

7.402 

172,175 

106,110 

0.415 

0.339 

173,174 

107,109 

0.470 

0.470 

173,190 

103 

6.626 

6 . 560 

179,189 

115,118 

1.511 

1.571 

180,138 

116,117 

0.523 

0.394 

181,137 

119,131 

9.230 

3.733 

132,136 

120,130 

2.548 

2.675 

133,185 

121,129 

0.100 

0.100 

184 

122,123 

9.326 

9.962 

191,194 

123,127 

0.473 

0.543 

192,193 

124,126 

0.654 

0.576 

195,20!) 

125 

6.931 

6.846 

190,199 

132,139 

0.100 

0.100 

197,193 

case  I case  b 


stress 
limit  (psl) 

nuithod 

Final  v’elght 
(16s) 

/ 

10000 

soc 

HOC 

29037 

29163 

0.177 
0.177 
0.243 
10.046 
0.10} 
2.705 
9.365 
0.671 
0.5  >0 
7.470 
2.460 

0.766 
7.639 
4.oi<) 
0.10) 
13.133 
1.647 
0.591 
7.753 
0.100 
0.130 
0.130 
0.169 
3.329 
0.133 
4.109 
13 .646 
0.629 
1.702 
3.240 
5.345 
3.549 
10.657 
17.723 
7.745 


nuidber  of 


0.274 
0.274 
0.315 
9.613 
).27  0 
2.909 
10.514 
o.()15 
0.593 
7.334 
2.495 
0.042 
7.045 
4.145 
0.100 
13.937 

1.460 

0.722 

7.554 

0.100 

0.256 

0.256 

0.274 

7.743 

0.231 

4.342 
14.464 

0.301 

1.557 

3.540 

5.343 
3.015 
9.367 

13.313 

7.330 


f 


TABLE  3.23 


WING  CARRY  THROUGH  BOX 


General  data  : 


! 


Nodal  coordinates 


R = 1.6  10  psi 
3 


p ■ O.ib  lb /In 
0 


lOOOOO  psl 
2 


a = 0.01  In 


a°  =«  50.0  In^ 


161 . 


L 


(Continue) 


TABLE  3.25 


WING  CARUV  through  box 


WEIGHT  CONVERGE"CE 


iteration 


Weight  tibs) 


no  flex. -constr.  case  B;  witn  flex,  const. 


Ucrke 
I B8I 


30214.0 

6360.0 

5386.0 

5615.0 

5335.0 

5262.0 

5255.0 

5234.0 

5272.0 

5239.3 

5201.1 

5164.1 
5131.8 

5104.6 

5032.3 

5064.1 

5049.3 

5049.7 
5051.5 

5053.2 

5054.3 

5056.2 

5057.4 

5053.4 

5059.2 

5059.3 

5060.4 
5060.3 
5060.0 

5059.7 


SAIICEF 


30214.3 
7630. 1 
6590.6 

6393.0 

6269.9 

6246.4 

6199.1 

6158.5 

6126.2 

6123.0 

6121.0 

6119.3 
6119.0 

6113.4 

6117.9 

6117.6 

6117.5 


Schmlt 
[ S13) 


1 1022.. 3 

9550.6 

7544.2 

6306.3 

6456.7 

6234.9 

6201. s 

6160.4 

6159.9 
6140.6 
6123.3 

6120.9 


ISerkc 
[ 831 


10214.0 

7577.0 

6134.0 

6923.0 

6301.0 

6609.0 

6473.0 

6333.0 

6333.0 

6292.5 

6262.6 
6240.5 

6230.7 

6215.7 

6220.1 

6253.7 

6236.3 

6300.4 

6301.7 

6296.0 

6236.1 
6273.9 

6261.1 

6243.4 

6236.5 

6225.7 

6216.2 

6207.7 
6200.4 
6194.1 


5042.8 

5041.1 

5039.5 

5037.9 

5036.2 

5034.5 


6160.9 

6 1 60 . 5 
6160.2 
6159.3 

6159.6 
6159.3 


TAiJLE  8.26 


WING  CARRY  THROUGH  BOX  - ‘'INAL  DESIGN 


cross  sections  (in^) 
element  flex,  constr.  ' 


case  B : with  flex,  constr. 


SAMCEF 


38.282 
35.018 
51.680 
54.4S7 
24.904 
23.443 
17.624 
20.547 
25.173 
26.344 
7.526 
8.772 
24 . 203 
20.630 
4.097 
2.454 
37.051 
37.133 
0.010 
0.010 
0.139 
0.103 
0.132 
1 . 100 
0.052 
2.592 
0.797 
4.607 
0.641 
2.630 
2.563 

5.345 

5.345 
6.130 

5. 345 


ochmlt  Berke 
I S18]  I 83) 


33.23 

35.03 

51.69 

54.49 

24.93 

23.46 

17.64 

20.52 

25.21 
26.82 

7.535 

8.301 

24.21 
20.63 

4,123 
2.495 
37.07 
37.14 
o.OlO 
0.010 
- 0.151 
0.067 
0.137 
1.085 
0.065 
2.574 
0,304 
4.532 
0.670 
2.651 
2.530 
5.829 
5.330 
6.122 
5.339 


25.44 

28.69 

17.73 

20.75 

25.32 

27.49 

7.62 

8.32 
24.62 
20.03 

4.16 

2.30 

37.53 

37.65 

0.01 

0.01 

0.07 

0.01 

0.03 

1.22 

0.07 

2.33 
0.31 
4.87 
0.51 

2.69 

2.70 
5.39 
5.32 
6.10 
5.32 


SA’ICLF 

ichm  It 

1 .;ni 

IJf-rkf 

1 B.iJ 

37.629 

37.55 

36 .8o 

36.386 

36.40 

36. 9v) 

52.535 

52.00 

53.33 

53.302 

53.76 

53.31 

23.672 

23.79 

24.13 

28.965 

28.9 

5 

27.82 

17,255 

17.26 

17.35 

21.231 

21.40 

22.00 

26.1  )0 

26.06 

23.42 

25.126 

25.15 

25.95 

3.704 

3.734 

9.44 

9.007 

3.''66 

9.32 

23.422 

23.43 

22.37 

1''.504 

n.57 

13.59 

5.172 

5.165 

5.79 

2.94S 

2.956 

4.47 

37.121 

37.0 

7 

36.39 

37.345 

37.3 

) 

37.52 

0.010 

0.0 

10 

0.01 

0.010 

0.0 

lo 

o.Ol 

0.099 

0.2 

13 

0.15 

0.061 

0.1 

70 

0.01 

0.010 

0.0 

10 

1.01 

0.010 

1.0 

1 ) 

0.18 

0.010 

0,0 

10 

0.01 

4.  142 

4.1 

91 

6.11 

0.911 

0.9 

15 

0.01 

3.251 

3.2 

85 

4.43 

0,010 

<).o 

lo 

1.15 

0.263 

7.8 

61 

6.94 

3.964 

7.7 

99 

9.76 

3.394 

9.3 

)0 

11.03 

3.932 

9.2 

29 

3.09 

9.342 

9.7 

69 

11.59 

8.983 

9.2 

30 

.3.09 

(continue) 


cross  sections 

1 

c lemcnt 

case  A:  no  flex,  constr. 

SAMCEF 

SclimiC 

I S13J 

Uerke 

1 33] 

2.762 

2.32 

2.568 

2.70 

2.704 

2.71 

2.604 

2.603 

2.70 

40 

5.731 

5.736 

5.33 

41 

5.314 

5.321 

5.30 

42 

16.437 

16.45 

16.60 

43 

13.734 

18.30 

13.99 

44 

10.936 

11.01 

11.25 

45 

13.331 

13.40 

13.66 

46 

11.434 

11.42 

11.60 

47 

5.948 

5.961 

6.03 

48 

12.143 

12.16 

12.24 

40 

14.266 

14.25 

14.40 

50 

7.223 

7 . 240 

7.26 

51 

7.451 

7.416 

7.35 

52 

5.510 

5.501 

5.62 

53 

0.407 

0.566 

0.01 

54 

3.630 

3.639 

3.67 

55 

9.636 

9.631 

‘J.97 

56 

4.392 

4.375 

4.54 

57 

0.245 

0.310 

0.03 

53 

0.127 

0.051 

0.01 

59 

0.129 

0.125 

0.01 

60 

0.010 

0.010 

0.01 

61 

0.010 

0.010 

0.01 

62 

0.010 

0.010 

o.ol 

63 

0.010 

0.015 

0.01 

We  iglit 
(ll)S) 

4075.1 

4076,0 

5034.5 

nbr.  of 
.ina  lyses 


15 


14 


50 


(In^) 


case  B;  with  flex. 

constr. 

S dim  It 

1 3131 

IJprko 

1 1)31 

3.679 

3.257 

7.30 

8.066 

7.301 

9.77 

3.133 

7 ,33  5 

6.03 

3.034 

7.352 

0.77 

3.034 

0.134 

10.02 

a.9b3 

n.  131 

i.OO 

25.475 

25.23 

24.61 

27.310 

27.07 

24.54 

19.047 

10.35 

19.78 

20.038 

21.13 

21.03 

16.503 

16.81 

17.19 

12.462 

12.65 

12.93 

18.732 

13.55 

13.00 

20.039 

10.91 

10.67 

6 ,666 

D .65) 

0.73 

5.636 

5.753 

7.40 

3.102 

3.12:5 

0.50 

3.35b 

5.64  2 

0.01 

5.946 

5.086 

6.79 

11. '’66 

11.03 

8.90 

5.377 

5.343 

3.30 

1.503 

1.601 

3 . 38 

0.010 

0.01) 

0.01 

0.010 

0.010 

O.ol 

0.010 

0.01  ) 

).oi 

0.010 

0.010 

0.  )1 

0,010 

0.010 

o.ol 

0.010 

0.010 

0.01 

6117.5 

6120.0 

6150.3 

17 


13 
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r ^ 


FOUR  LEVEL  SPACE  TRUSS 

General  data  ; 

E = 10^  psl 
0 “0.1  Ib/in^ 

0 “ 2500)  psl 

a - 0.1  lii^ 
a°  ■ 1.0  in^ 


TABLE  8.33 

FOUR  LEVEL  SPACE  TRUSS  > WEIGHT  CONVERGENCE 


Iteration 


I Weight  (lbs) 

SAMCEF, 
Taig  & 
Kerri  T31 

Gellatly 
& Berke 
IG13] 

Venkayya 

IV31 

Berke 
& Khot 
IB31 

■ ■ - ■ - > ■ ' 

Scimilt 
Far  sill 

1 SI) 

Jcli:-Iit 
<5  Mlur.i 

1 SI  !1 

— 

656.77 

337.25 

379.7vS 

379.63 

379.66 


656.77 

416.07 

406.21 

399.06 

396.32 

396.25 

396.02 

395.97 


656 .3 

478.6 
455.0 
446.9 
445.5 

445.4 

401.7 

391.5 

333.6 


331.2 


656.77 

337.01 

379.67 

379.37 


309 . 1 2 
333.09 
796.16 
763.61 

736.69 

716.63 

703.77 

645.07 
616.97 
525.29 
49  1 .96 

463.69 
450.22 

433.77 
423.94 
413.65 

404.08 
397.43 
393.33 
333.14 

333.63 


731.15 
477.95 
397.43 
3 ’.3.27 
3 ; '.47 
379.36 

379.63 

379.64 
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TABLE  8.34 

FOUR  LEVEL  SPACE  TRUSS  - FINAL  DESIGN 


Cross  sections  (In  ) 


ilements 

SAMCEF, 
Taig  & 
Kerri  T3] 

Gellatly 
& Berke 
IP13] 

Venkayya 

[V31 

Berke 
& Kliot 

1 B8] 

Sclunlt 
& Pars  lit 

1 So] 

[ Sclimlt 
u Miura 

1 318] 

1-4 

1.8973 

1.4636 

1.313 

1.8931 

2. 07')  4 

1.335 

5-12 

0.5158 

0.5207 

0.523 

0.5171 

0.5034 

0.5125 

13-16 

0.1000 

0 . 1000 

0 . 100 

0 . 1000 

0 . 1000 

0.1000 

17,18 

0.1000 

0.1000 

0.100 

0 . 1000 

0 . 1000 

0.1000 

19-22 

1.2801 

1.0235 

1.246 

1.2793 

1.1067 

1.267 

23-30 

0.5143 

0.5421 

0.524 

0.5149 

0.5792 

0.5113 

31-34 

0.1000 

0 . 1000 

O.lOO 

0.1000 

0.1000 

o.iooo 

35,36 

0.1000 

0.1000 

0.100 

0 . 1000 

0 . 1000 

0.10)0 

37-40 

0.5067 

0.5521 

0.611 

0.5032 

0.2643 

0.5233 

41-43 

0.5200 

0.6034 

0.532 

0.5196 

0.5430 

0.5173 

49-52 

0.1000 

0 . 1000 

0 . 100 

0 . 1000 

0.10)0 

0.1000 

53,54 

0 . 1000 

0 . 1000 

0 . 100 

0.10)0 

0.1500 

0.1000 

55-58 

0.1571 

0.1492 

0.161 

0.1571 

0.1535 

0.1565 

59-66 

0.5356 

0.7733 

0.557 

0.5335 

0.5936 

0.5453 

67-70 

0.4099 

0.4534 

0.377 

0.4156 

0.3414 

0.4105 

71,72 

0.5690 

0.3417 

0.506 

0.5510 

0.6076 

J.569Q 

ight  (lbs) 

379.66 

395.97 

331.2 

379.67 

333.63 

379.64 

number  of 
analyses 


5 


8 


11 


22 


TABLE  3,36 


TRANSMISSION  TOWER 


General  data  : 


E =*10  psi 
p “0,1  lb/ In 


a^  “ 0,01  in 
a°  - 2.0  tn^ 


Allowable  stresses 


e lements 


stress  limits  (psi) 

traction 

compress  ion 

40000 

- 35092 

40000 

- 11590 

4f)000 

- 17305 

40000 

- 35092 

stress  limits  (psi) 


elements 


12,13 

14-17 

13-21 

22-25 


traction 


40000 

40000 

40000 

40000 


compress  Ion 


load 

node 

case 

1 

1 

2 

3 

6 

2 

5 

6 

load  component  (lbs) 


10000 

10000 

0 

0 


20000 

20000 


Flexibility  constraints 


max,  displacement  (in) 


1 

2 


35 

35 


35 

35 


/ 

/ 


w 


T 


TABLE  8.39 


179. 


TRANSMISSION  TOWER  > FINAL  DESIGN 


Cross  sections  (In*) 

elements 

SAMCEF 

Gellatly 
& Berkc 
[G13] 

Venkayya 

Schniit 

A Pars hi 

1 33) 

Sclimlt 
f.  Mltir.i 

I 8131 

SOC 

HOC 

IV31 

1 

0.0100 

0.0100 

0.0100 

0.028 

0.010 

0,010 

2-5 

1.9761 

1.9363 

2.0069 

1.942 

1.064 

1.085 

6-9 

3.0108 

2.9946 

2.9631 

3.031 

3.033 

2.006 

10,11 

0.0100 

0.0100 

0.0100 

0.010 

0.010 

0.i)10 

12,13 

0.0100 

0.0100 

0.0100 

0.010 

0.010 

0.010 

lA-17 

0.6848 

0.6840 

0.6376 

0.603 

0.670 

0.6  ’.4 

13-21 

1.6786 

1.6770 

1.6784 

1.67.; 

1.630 

1,677 

21-25 

2.6565 

2.6618 

2.6638 

2.627 

2.670 

2.662 

Weight  (lbs) 

545.21 

545.16 

545.36 

545.49 

545.23 

545.17 

number  of 
analyses 

10 

7 

7 

b 

17 

10 

180 


5 9 3 


57  bars 
42  D.O.F. 


FIGURE  8.40 


I 


TABLE  8.46 

SIMPLIFIED  AIR  BRAKE 


clement 

clement 

SOC 

HOC 

SOC 

1 

367.13 

1010.43 

20 

1156.18 

2 

377.25 

1016.70 

21 

1.00 

3 

1.00 

1.00 

22 

1.00 

4 

364.16 

1006.72 

23 

1.00 

5 

1.00 

2.59 

24 

1.00 

6 

1.23 

2.59 

25 

39,30 

7 

1.00 

1.00 

26 

197.75 

8 

1.00 

1.00 

27 

4.22 

9 

1.00 

3.16 

28 

62.30 

10 

390.03 

1.00 

29 

262.04 

II 

1.00 

1.00 

30 

20.07 

12 

391.94 

4.23 

31 

l.'X) 

13 

1169.00 

1349.89 

32 

M.03 

14 

1.00 

1.00 

33 

479.73 

15 

1.00 

1.00 

34 

1.00 

16 

9.65 

16.66 

35 

l.f)0 

17 

1.00 

1.00 

36 

922.38 

18 

1.00 

1.00 

37 

llb2.33 

19 

1.00 

50.30 

38 

921.93 

method 


deflections  (mm) 


node  I 


node  2 


nod( 


SOC 

HOC 


8.387 

3.393 


7.837 

7.393 


8 

8 
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- FINAL  DESIGNS 


element 

section 

(mm^) 

HOC 

SOC 

1341.36 

39 

l.'» 

1.00 

l.fV) 

40 

1.19 

1.00 

l.oo 

41 

1.00 

7.64 

1.00 

42 

68.20 

1 .00 

1.00 

43 

1.00 

1.00 

28.26 

44 

79.61 

1 .09 

227.. 61 

45 

31.69 

52.51 

1.00 

46 

1169.11 

1354.52 

1.00 

47 

l.oo 

1.00 

232.51 

48 

75.53 

48.54 

3.44 

49 

1.00 

1.00 

1.00 

5') 

914.76 

1)1(}.74 

62.91 

51 

926.11 

1025.11 

473.38 

52 

1.1)0 

1.0) 

1.00 

53 

346.30 

1012.06 

1.00 

54 

44.04 

20.70 

59.26 

55 

13.73 

1.00 

1339.56 

56 

1.00 

1.00 

59.33 

57 

6.42 

5.12 

we ight 

nbr  of 

(kg) 

analyses 

3 


20.51 

15.563 


20 

19 


COMPUTING  TIME  - I.B.M.  370/158 


1 


TABLE 
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structural  optimization  by  finite  element. (U) 

JAN  70  C FLEURYr  0 SANDER  AF0SR-77-31ie 

UNCLASSIFIED  L.T.A.S.  NBR-SA-SO  AFOSR*TR-7S-O670  NL 


An 

Ifi 


10 


15 


Nbr 

ANALYSES 

— ■ ► 


WING  BOX 


FIGURE  8.50 


1 

11 

193. 

TABLE  8.52 

1 

\ 

WING  BOX  - FINAL  DESIGN 

! 

i 

FSD  GOC 


design 

variable 


2 

section  (In  ) 
thickness  (In) 

max.  stress. 

(psl) 

-> 

section  (In'") 
thickness  (in) 

■ax.  stress. 

(psl) 

1.735 

10000 

0.  lOO  )•) 

10000 

0.213 

9139 

0.23409 

7700 

0.195 

- 9467 

O.IOJOO 

- 9503 

0.119 

9175 

0.13936 

3207 

0.022 

8426 

0.02000 

3341 

0.024 

9224 

0.02000 

9 2;35 

0.109 

- 3619 

o.loooo 

- 4105 

0.130 

9193 

0.16670 

10000 

0.109 

5450 

0. 10000 

5975 

0.073 

9165 

0.09005 

3631 

0.041 

9214 

0.03273 

1000) 

0.022 

4937 

0.02f)00 

6106 

0.109 

- 4802 

0.10000 

- 5254 

0.065 

9174 

0.06235 

I'V))) 

0.022 

9030 

0.0200  ) 

9662 

0.045 

9200 

0.04222 

lOOOr) 

90 


I = 

' V ' 

4 deslfii 
viritblts 


e«mp«sit  tltMtnt 


U orthotropic  mofflhronos 


FIGURE  8.53  (continued) 


f 

TABLE  8.54 

COMPOSITE  BOX  BEAM 

General  data 


Material 

e 

s 

alunlnlum 

l>oroii  - epoxy 

Young  tnodulus 

e 

• 

Hi 

Hi- 

1.03  10 

1.05  10^ 

3.0  10  psi 
2.7  10^  psi 

Poisson's  ratio 

• 

« 

Hi 

0.1 

* * 6 

0.21 

Shear  modulus 

• 

s 

G » 

4.033  10 

0.7  10  psi 

Mass  density 

• 

• 

P ■ 

0.1 

0.0725  Ih/in 

Min.  cross  section 

• 

• 

0.0 J5  In 
^ 0.005  in 

/ 

Min.  thickness 

• 

• 

£ =• 

0.005  in 

(per  layer) 

Initial  cross  section 

• 

1.0 

/ 

Initial  thickness 

• 

• 

o 

a - 

^ 0.1 

0.025 

(per  layer) 

Load  tng 


prob 1cm 

node 

lead  comoonent  (lbs) 

X 

y 

Z 

A 

1,2. 3, 4 

0 

0 

- 1000 

B 

1,2 

0 

0 

- 1000 

3,4 

0 

0 

- 975 

(continue) 


tabu;  3.56 


COMPOSITE  BOX  BEAM  - CASE  A - WEIGHT  CONVERGENCE 


SAMCEF  (primal  algorithm) 

weight 

nodal  deflection 

(In) 

(I'^s) 

HQH 

0-12 

20-32 

135.6 

4.313 

2.902 

0.300 

38.28 

6.970 

3.973 

0.300 

34.73 

6.995 

3.930 

0.300 

32.00 

7.000 

3.007 

0.300 

31.76 

6.990 

3.093 

0.300 

30.73 

6.996 

3.994 

0.300 

29.07 

6.997 

3.995 

0.300 

20.33 

6.993 

3.996 

0.300 

28.83 

7.000 

3.906 

0.300 

23.45 

7.000 

3.907 

0.300 

Khot,  Venkaj'ya  & lierkc  | K15I 


weisht  nodal  deflection  (In) 
(lbs) 


186.  )0 
50.50 

40.01 
34.65 
32.43 
31.04 
31.71 
31.55 
31.30 
31.23 
20.48 
28.78 
2.', .02 

28.02 


4.440 

6.262 

6.130 

b.085 


2.030 

3.571 

3.514 

3.402 


20-32 


0.30) 


6.210 

6 . 25b 
f).()84 
6.027 
6.073 
6.974 


3.561 

3.58  2 
3.124 
).0(,7 
4.0  )>) 
4.  ' )0 


i).30O  I 


f).  1 )i) 


O.2O0 


200 


TABLE  8.57 


COMPOSITE  BOX  BEAM  - CASE  A - FINAL  DESIGNS 


201 . 


TABLE  3.58 


COMPOSITE  BOX  BEAM  - CASE  B - WEIGHT  CONVERGENCE 


’I 


iturat ion 

SAMCEF  (dual  algorithm) 

Kliot 

, Venkfiyyn  «i 

Berko  1 41 5| 

weight 

(Ihs) 

nodal  deflection  (in) 

weight 

(lbs) 

nodal  deflection  (in) 

1.2 

3,4 

1,2 

3.4 

1 

42.74' 

14.00 

13.03 

50.57 

14.00 

13.33 

2 

18.31 

13.85 

15.00 

22 .05 

14  .'X) 

13.33 

3 

16.62 

13.03 

15.00 

4 

15.00 

13.07 

15.00 

5 

15.34 

13.07 

15.00 

14.70 

14.00 

13.54 

0 

14.33 

13.08 

15.00 

7 

14.51 

13.00 

15.00 

3 

14.24 

13.90 

15.00 

0 

14.05 

13.09 

15.00 

10 

13.04 

14.00 

15. ')0 

14.76 

14.00 

11.11 

20 

14 .66 

14.00 

14.  lo 

30 

14.62 

14.00 

14.6  3 

40 

14.60 

14.00 

14.82 

50 

14.50 

14.00 

14.01 

60 

14.58 

14.00 

14.96 

70 

14.58 

14.00 

14.93 

30 

14.53 

14.00 

14.99 

00 

14.58 

14.00 

14.99 

l<vo 



14.53 

14.00 

14.96 

CO>f?OSITE  BOX  BEAM  - CASE  B - FINAL  DESIGN 


analyses 


AIRBRAKE  STRUCTURE 


First  simplified  model 
(27  elements) 


1 bar 

26  membranes 
43  D.O.F. 


FIGURE  8.62 


AIRBRAKE  STRUCTURE 


SECOND 


SIMPLIFIED  MODEL 
(125  elements) 


Figui ^ 8.65 


AIRBRAKE 


212 


REFERENCES 


( A1  ] ABADIE,  J. 

"On  the  Kuhn-Tucker  Theorem" 

Nonlinear  Programming,  ABAUIE,  J.  (ed.),  North-Holland,  Amsterdam,  p.  10-36 
(1967) 

( A2  J ABADIE,  J.  & CARPENTIER,  J. 

"Generalization  of  the  Wolfe  reduced  gradient  method  to  the  case  of  nonlinear 
constraints" 

Optimization,  FLETCHER,  R.  (ed.).  Academic  Press,  London,  p.  37-49  (1969) 

[ A3  1 ABLOW,  C.M.  & BRIGHAM,  G. 

"An  analog  solution  of  progranuning  problems" 

Operat.  Res,,  3:383-394  (1955) 

( A4  ) ARMAND,  J.L. 

"Minimum-mass  design  of  plate-llke  structure  for  specified  fundamental  fre- 
quency" 

AIAA  J.,  9(9);1739-1745  (1971) 

( A5  ] ARMAND,  J.L.  & LODIER,  B. 

"Optimal  design  of  bending  elements" 

ASME  Energy  Technology  Conference,  Structural  Optimization  Methods  Session, 
VENKAYYA,  V.B.  (ed.),  Houston,  p.  68-78  (1977) 


( B1  J BECKMAN,  F.S. 

"Tlie  solution  of  linear  equations  by  the  conjugate  gradient  method" 
Mathematical  Methods  for  Digital  Computers,  Vol.  1,  A.  RALSTON  & 11. S.  WILF 
(ed.),  John  Wiley,  New  York  (I960) 

( B2  J BROYDEN,  C.G. 

"Quasi-Newton  methods  and  their  application  to  function  minimization" 

Math.  Comp.,  21:368-381  (1967) 

( B3  I BEALE,  E.M.L. 

"Numerical  Methods" 

Nonlinear  Programming,  ABADIE,  J.  (ed.),  North-Holland,  Amsterdam,  p.  133- 
205  (1967) 

[ B4  ) BROWN,  D.M.  & ANG,  A.H.S. 

"Structural  optimization  by  nonlinear  programming" 

J,  of  the  Structural  Division,  ASCE,  92,  ST6: 319-340  (1966) 

( B5  I BALDUR,  R. 

"Structural  optimization  by  inscribed  hyperspheres" 

J.  of  the  Engr.  Mech.  Division,  ASCE,  93,  EM3:503-508  (1972) 


( B6  ) BARTA,  J. 

"On  the  minimum  v/elght  of  certain  redundant  structures" 

Acta,  Tech.  Acad.  Scl.  Hungar.,  13:67-76  (1957) 

( B7  1 BARNETT,  R.L.  & HERMAN,  P.C, 

"High  performance  structures" 

NASA  Contractor  Report,  CR-1033,  Washington,  D.C,  (1068) 

I B3  1 BERKE,  L.  & KHOT,  N.S. 

"Use  of  optimality  criteria  methods  for  large  scale  systems" 

ACARD  Lecture  Series,  AGARD-LS-70  (1974) 

{ B9  ] BARNETT,  R.L. 

"Minimum  weight  design  of  beams  for  deflection" 

J.  of  the  Engr.  Mech.  Division,  ASCE,  87,  EM1:75-109  (1961) 

( BIOJ  BERKE,  L. 

"An  efficient  approach  to  the  minimum  weight  design  of  deflection  limited 
structures" 

USAF  Technical  Memorandum,  AFFDL-'ni-70-4-FUTR  (1070) 

(Bill  BERKE,  L. 

"Convergence  behaviour  of  iterative  resizing  procedures  based  on  optimality 
criteria" 

USAF  Technical  Memorandum,  AFFDL-TM-72-1-FBR  (1972) 

(B12)  BERKE,  L.  & VENKAYYA,  V.B. 

"Review  of  optimality  criteria  approaches  to  structural  optimization" 

ASME  Winter  Annual  Meeting,  New  York  (1974) 

(BU)  BEREZIN,  I.S.  & ZHIDKOV,  N.P. 

"Computing  methods",  Vol.  2 
Pergamon,  Oxford  (1965) 

[B14)  BARTHOLOMEW,  P.  & ^10RRIS,  A.J. 

"A  unified  approach  to  fully-stressed  des,lgn" 

Engineering  Optimization,  2:3-15  (1976) 

1B15)  BERKE,  L. 

Private  communication  (1977) 


(Cl  J CAUCHY,  A. 

"Methode  gcncrale  pour  la  resolution  des  systemes  d'equatlons  slmultaiices" 
C.R.  Acad,  Sci.  Paris,  25:536-538  (1847) 

(C2  1 CURRY,  H.B. 

"The  method  of  steepest  descent  for  nonlinear  minimization  problems" 

Quart.  Appl.  Math,,  2(3):250-261  (1944) 

(C3  ) CURTIS,  A.R.  4 REID,  J.K. 

"On  the  automatic  scaling  of  matrices  for  gausslan  elimination" 

J.  Inst.  Maths.  Applies.,  10:118-124  (1972) 


214. 


( C4  J CROCKETT,  J.B.  * CHERNOFF,  II. 

"Gradient  methods  of  maxim izat ion"  | 

Pacific  J.  Math.,  5:33-50  (1055)  ] 

1 C5  J CAROLL,  C.W.  | 

"The  created  response  surface  technique  for  optimizing  nonlinear  restrained  | 

systems"  I 

Operat.  Res,,  0:169-184  (1961)  | 

(C6  1 COX,  H.L. 

"The  design  of  structures  of  least  weight"  | 

Pergaroon,  Oxford  (1965)  j 

I C7  J CILLEY,  F.H. 

"The  exact  design  of  statically  indeterminate  frameworks,  an  exposition  of 
its  possibility,  but  futility" 

Trans.  ASCE,  42:353-407  (1900) 

( C3  1 CHERN,  J.M.  & PRAGER,  W. 

"Optimal  design  of  beams  for  prescribed  compliance  under  alternative  loads" 

J.  Optlm.  Theory  Appl.,  5(6):424-431  (1970) 

1C9  ] CHERN,  J.I?.  & PRAGER,  W. 

"Minimum  weight  design  of  statically  determinate  trusses  subject  to  multiple 
constraints" 

Int.  J.  Solids  & Structures,  7:931-940  (1971) 

[CIOJ  CLAUSEN,  T. 

"Ober  die  Form  architektonlscher  Saulen" 

Bulletin  physlco-mathematique  de  I'Acaderaie,  9:279-294  (1R51) 


[D1  ] DAVIDON,  W.C. 

"Variable  metric  method  for  minimization" 

Research  and  Development  Report  ANL-5990  (Revised),  Argonne  National  Labora- 
tory, U.S,  Atomic  Energy  Commission  (1959) 

( D2  J DANTZIG,  G.B. 

"Linear  programming  and  extensions" 

Princeton  University  Press,  Princeton,  N.J,  (1963) 

[ D3  ] DANIEL,  J.W. 

"The  conjugate  gradient  method  for  linear  and  nonlinear  operator  equations" 
SIAM  J.  Numer.  Anal.  4(l):10-26  (1967) 

[D4  ] DIVYER,  W.,  ROSENBAUM,  J.,  SHULMAN,  M.  4 PARDO,  H. 

"Fully  stressed  design  of  airframe  redundant  structures" 

Proc.  Second  Conf.  on  Matrix  Methods  in  Structural  Mechanics,  Dayton,  Ohio 
AFFDL-TR-63-150: 155-181  (1968) 


215. 


(D5  1 DWYER,  W.J.,  EMERTON,  R.K.  & OJALVO,  I.V. 

"An  automated  procedure  for  the  optimization  of  practical  aerospace  structu-  | 
res"  I 

USAF  Technical  Report,  AFFDL-TR-70-118  (1971)  { 

[06  1 DWYER,  W.J.  i 

"An  improved  automated  structural  optimization  program" 

USAF  Technical  Report,  AFFDL-TR-74-96  (1974) 


[FI  J FIACCO,  A.V.  & Me  CORMICK,  G.P. 

"Nonlinear  programming  : sequential  unconstrained  minimization  techniques" 
John  Wiley,  New  York  (196b) 

[ F2  1 FARKAS,  J. 

"Ober  die  Theorie  der  Elnfachen  Ungleichungen" 

J.  Reine  Angew.  Math.,  124:1-27  (1905) 

[ F3  J FLETCHER,  R. 

"Function  minimization  without  evaluating  derivatives  - a review" 

Computer  J.,  8:163-168  (1963) 

[ F4  J FLETCHER,  R.  & REEVES,  C.M, 

"Function  minimization  by  conjugate  gradients" 

Computer  J.,  7:149-154  (1964) 

[F5  J FLETCHER,  R.  & POWELL,  J.D. 

"A  rapidly  convergent  descent  method  for  minimization" 

Computer  J.,  6:163-163  (1963) 

[ F6  ] FOX,  R.L.  & STANTON,  E.L. 

"Developments  in  structural  analysis  by  direct  energy  minimization" 

\ AIAA  J.,  6 (6): 1036-1042  (1968) 

I [F7  ] FADDEEV,  D.K.  & FADDEEVA,  V.N. 

I "Computational  methods  of  linear  algebra" 

^ Freeman  (1963) 

[ [ F8  1 FOX,  R.L. 

I "Unconstrained  minimization  approaches  to  constrained  problems" 

i Structural  Design  Applications  of  Mathematical  Programming  Techniques,  POPE, 

I G.G.  i SCHMIT,  L.A.  (ed.),  AGARDOGRAPH-149,  p.  55-78  (1971) 

L I F9  ] FLETCHER,  R. 

> "Generalized  inverses  for  non-linear  equations  and  optimization" 

Numerical  Methods  for  Non-linear  Algebraic  Equations,  RABINOWITZ,  P.  (ed.), 
Gordon  and  Breach,  London  (1970) 

[ FIOJ  FLETCHER,  R. 

"Minimizing  general  functions  subject  to  linear  constraints" 
i.  Numerical  Methods  for  Non-linear  Optimization,  L00TSM/\,  F.A.  (ed.).  Academic 

[ Press,  p.  279-296  (1971) 


I 


216. 


( FllJ  FRICH,  R. 

"The  logarithmic  potential  method  for  solving  linear  programming  problems" 
Memorandum  of  the  University  Institute  of  Economics,  Oslo  (1955) 

( F12)  FALK,  J.E. 

"Lagrange  multipliers  and  nonlinear  programming" 

J.  Math.  Anal.  Appls.,  19:141-159  (1967) 

1 F13)  FOX,  R.L. 

"Optimization  methods  for  engineering  design" 

Add Ison-Wesley , Reading,  Mass.  (1971) 

[ F14)  FOX,  R.L. 

"Constraint  surface  normals  for  structural  synthesis  techniques" 

AlAA  J.,  3(8):1517-1518  (1965) 

( F15)  FLEURY,  C. 

"Methodes  numeriques  d' optimisation  des  structures" 

LTAS  (Unlverslte  de  Liege),  Rapport  SF-19  (1973) 

[ F16]  FLEURY,  C.  & FRAEIJS  de  VEUBEKE,  B. 

"Structural  optimization" 

Proc.  Sixth  IFIP  Conf.  on  Optlm,  Tech.,  Novosibirsk,  1974  in  "Lecture  Motes 
In  Computer  Science",  Springer  Verlag,  27  : 314-326  (1975) 

( F17J  FLEURY,  C. 

"Applications  pratiques  des  methodes  d 'optimisation  structurale" 

LTAS  (Universlte  de  Liege),  Rapport  SF-30  (1974) 

1 F18]  FLEURY,  C. 

"Optimisation  d'un  spoiler  d'avlon  avec  verification  de  restrictions  de  fle- 
xlbillte" 

LTAS  (Universlte  de  Liege),  Rapport  SF-33  (1975) 

( F19]  FLEURY,  C. 

"Optimisation  des  structures  par  la  raethode  des  elements  finis" 

Coll.  Publ.  Fac.  Sc.  Appl.  (Universlte  de  Liege),  52:63-102  (1976) 

[ F20]  FLEURY,  C.  4 GERADIN,  Mi 

"Application  des  elements  finis  a 1' optimisation  des  structures" 

Methodes  Numeriques  en  Calcul  Sclentlflque  et  Technique,  AFCET-GAt-lNI-GAMI 
(1976) 

[F21]  FLEURY,  C. 

"Algorlthmes  de  prograraraat ion  mathematique  non  llnealre  en  optimisation  strui 
turale" 

LTAS  (Universlte  de  Liege),  Rapport  SF-52  (1976) 

[ F22)  FLEURY,  C. 

"Relations  entre  I'approche  par  crlteres  d'optlmalitc  et  la  programmatlon 
convexe  en  optimisation  structurale" 

Coll.  Publ.  Fac.  Sc.  Appl.  (Unlverslte  de  Liege),  66:3-64  (1977) 

( F23J  FLEURY,  C.  & GERADIN,  M. 

"Optimality  criteria  and  mathematical  programming  in  structural  weight  optl- 
mlzat Ion" 


Computers  & Structures,  3:7-17  (1973) 


217. 


[ F24)  FLEURY,  C.  & SANDER,  G. 

"Relations  between  optimality  criteria  and  mathematical  programming  in  struc- 
tural optimization" 

Proc.  Symp.  "Applications  of  Computer  Methods  in  Engineering",  Univ.  Soutlicri' 
California,  Los  Angeles,  1:507-520  (1977) 

( F25J  FOX,  R.L.  & KAPPOOR,  M.P. 

"Rates  of  change  of  eigenvalues  and  eigenvectors" 

AIAA  J.,  6(12): 2426-2429  (1968) 

I F261  FOX,  R.L.  & KAPPOOR,  M.P. 

"Structural  optimization  in  the  dynamics  response  regime  : a computational 
approach" 

AIAA  J.,  8(10): 1798-1804  (1970) 

I F27J  FOX,  R.L.  & SCHMIT,  L.A. 

"Advances  in  the  integrated  approach  to  structural  synthesis" 

J.  Spacecraft  Rockets,  3:858-366  (1966) 

[ F28]  FRAUENTHAL,  J.C. 

"Constrained  optimal  design  of  circular  plates  against  buckling" 

Harvard  University,  Cambridge,  Report  no.  SM-50  (1971) 

I F29] 

FLEURY, C. 

"Le  dlmensionnement  automatique  des  structures 
elastlques" 

Doctoral  dissertation,  LTAS,Unlversitc  de 
Liege  (1978) 


[Gl  1 GOLDSTEIN,  A. A. 

"Cauchy's  method  of  minimization" 

Numer.  flath.,  4:146-150  (1962) 

(G2  1 GOLDSTEIN,  A.A. 

"Constructive  real  analysis" 

Harper,  New  York  (1967) 

(r,3  J GOLDFARB,  D. 

"Extensions  of  Newton's  method  and  Simplex  methods  for  solving  quadratic  pro- 
grams" 

Numerical  Methods  for  Non-linear  Optimization,  LOOTSMA,  F.A.  (ed.).  Academic 
Press,  p.  239-254  (1971) 

[G4  ] GOLDFARB,  D, 

"Extension  of  Davidon's  variable  metric  method  to  maximization  under  linear 
inequality  and  equality  constraints" 

SIA>!  J.  Appl.  Math.,  17:739-764  (1969) 


218 


[ G5  J GILL,  P.E.  & MURRAY,  W. 

"Newton-type  methods  for  linearly  constrained  optimization" 

Numerical  Methods  for  Constrained  Optimization,  GILL,  P.K.  & MURRAY,  W. 
(ed,).  Academic  Press,  p.  29-66  (1974) 

( G6  J GILL,  P.E.  et  MURRAY,  W. 

"Quasi-Newton  methods  for  linearly  constrained  optimization" 

Numerical  Methods  for  Constrained  Optimization,  GILL,  P.E.  & MURRAY,  W. 
(ed.).  Academic  Press,  p,  67-92  (1974) 

[ G7  1 GRIFFITH,  R.E.  & STEWART,  R.A. 

"A  nonlinear  programming  technique  for  the  optimization  of  continuous  proces- 
sing systems" 

Managmt  Scl.,  7:  379-392  (1961) 

( G8  ] GERARD,  G. 

"Minimum  weight  analysis  of  compression  structures" 

New  York  University  Press,  New  York  (1956) 

I G9  1 GELLATLY,  R.A.,  GALLAGHER,  R.H.  & LUEERACKI,  W.A. 

"Development  of  a procedure  for  automated  synthesis  of  minimum-weight  struc- 
tures" 

USAF  Technical  Report,  AFFDL-TDR-64-141  (1964) 

(GlOj  GELLATLY,  R.A.  & GALLAGHER,  R.H. 

"A  procedure  for  automated  minimum  weight  structural  design  - Part  1,  theoci- 
tlcal  basis" 

Aero.  Quart.,  17(3) :216-230  (1966) 

[GllJ  GELLATLY,  R.A.  & GALLAGHER,  R.H. 

"A  procedure  for  automated  minimum  weight  structural  design  - !’art  II,  appli- 
cations" 

Aero.  Quart.,  17(4) ;332-342  (1966) 

IG12J  GELLATLY,  R.A. 

"Development  of  procedures  for  large  scale  automated  minimum  welglit  structu- 
ral design" 

USAF  Technical  Report,  AFFDL-TR-66-180  (1966) 

[ G13J  GELLATLY,  R.A.  & BERKE,  L. 

"Optimal  structural  design" 

USAF  Technical  Report,  AFFDL-TR-70-165  (1971) 

(G141  GALLAGHER,  R.H. 

"Fully  stressed  design" 

Optimum  Structural  Design  - Theory  and  Applications,  GALLAGHER,  R.H.  & 
ZIENKIEWICZ,  O.C.  (ed.),  John  Wiley,  London,  chapltrc  3 : 19-32  (1973) 

(G15]  GELLATLY,  R.A.,  DUPREE,  D.M.  & BERKE,  L. 

"OPTIM  II  : A magic  compatible  large  scale  automated  minimum  weight  design 
program" 

USAF  Technical  Report,  AFFDL-TR-74-97  (1974) 

[G16)  GELLATLY,  R.A. 

"Tiie  role  of  optimization  in  the  design  of  aircraft  structures" 

Proc.  AGARD  Symposium  on  Structural  Optimization,  Istambul,  AGARD-CP-36-70  : 
paper  9 (1970) 


219. 


( HI  1 IIESTENES,  M.R.  & STIEFEL,  E.L. 

"Methods  of  conjugate  gradients  for  solving  linear  systems" 

J.  Res.  Nat.  Bur.  Standards,  Section  B,  49:409-436  (1952) 

I H2  J HIMHELBLAU,  D.M. 

"Applied  nonlinear  programming" 

Me  Graw-Hlll,  New  York  (1972) 

1 H3  J HESTENES,  M.R. 

"The  conjugate  gradient  method  for  solving  linear  systems" 

Proc.  of  Symposium  In  Applied  Math.,  VI,  Num.  Anal.  : 83-102  (1956) 

( H4  ] HAUZEUR,  M.A. 

"Optimisation  des  structures  - cas  de  restrictions  Imposces  sur  la  flcxlbi- 
lltc" 

These  de  maltrlse  - ENSA,  Paris  (1969) 

IH5  1 HAUG,  E.J.,  PAN,  K.C.  4 STREETER,  T.D. 

"A  computational  method  for  optimal  structural  design 

I.  Piecewise  uniform  structures" 

Int.  J.  Num.  Methods  Engr.,  5:171-184  (1972) 

lH6  J HUANG,  N.C.  & SHEU,  C.Y. 

"Optimal  design  of  an  elastic  column  of  thin-walled  cross  section" 

J.  Appl.  Mech.,  ASME,  25:285-288  (1968) 

[H7  ) HILL,  R. 

"The  mathematical  theory  of  plasticity" 

Oxford  University  Press,  London  (1950) 


I J1  ] JOHNSON,  D.  4 BROTTON,  D.M. 

"Optimum  elastic  design  of  redundant  trusses" 

J.  of  the  Structural  Division,  ASCE,  95,  ST12  : 2589-2610  (1969) 

1 J2  1 JONES,  R.M. 

"Mechanics  of  composite  materials" 

Me  Graw-Hlll,  New  York  (1976) 


( K1  1 KUHN,  U.W.  & TUCKER,  A.W. 

"Nonlinear  programming" 

in  J.  NEYMAN  (ed.).  Proceedings  of  the  Second  Berkeley  Symposium  on  Mathema- 
tical Statistics  and  Probability,  University  of  California  Press,  4J1-493 
(1951) 


220 


[ K2  1 KIEFER,  J. 

"Optimum  sequential  search  and  approximation  methods  under  minimum  ref’ularity 
assumptions" 

SIAM  J.,  5:105-136  (1957) 

( K3  1 KANTOROVITCH,  L.V. 

"Functional  analysis  and  applied  mathematics" 

Transl.  C.D.  BENSTER,  Dept,  of  Math.,  Univ.  of  California,  Los  Anqcles  (1952) 
[ K4  ) KELLEY,  J.E. 

"The  cutting-plane  method  for  solving  convex  programs" 

SIAM  J.,  8(4): 703-7 12  (1960) 

[ K5  1 KLEIN,  B. 

"Direct  use  of  extremal  principles  in  solving  certain  optimizing  problems 
involving  inequalities" 

J.  Operat.  Res.  Soc.  Am.,  3:163-175  (1955) 

( K6  1 KOWALIK,  J.S. 

"Feasible  direction  methods" 

Structural  Design  Applications  of  Mathematical  Programming  Techniques,  POPL, 
G.G.  & SCHMIT.  L.A.  (ed.),  AGARDOGRAPH  149:79-93  (1971) 

(K7  ] KARNES,  R.N.  6 TOCHER,  J.L. 

"Automatic  design  of  optimum  hole  reinforcement" 

BOEING  Report,  D6-23359  (1968) 

(K8  ) KAVLIE,  D.  & MOE,  J. 

"Application  of  nonlinear  programming  to  optimum  grillage  design  with  non 
convex  sets  of  variables" 

Int.  J.  Num.  Methods  Engr.,  l(4):351-373  (1969) 

I K9  J RICHER,  T.P. 

"Optimum  design  - minimum  weight  versus  fully  stressed" 

J.  of  the  Structural  Division,  ASCE,  92,  ST6  : 265-279  (1966) 

(KIOJ  KHOT,  N.S.,  VENKAYYA,  V.B.  & BERKE,  L. 

"Optimum  structural  design  with  stability  constraints" 

Int.  J.  Num.  ‘lethods  Engr.,  10:1097-1114  (1976) 

( Kll]  KAMAT,  M.P.  & SIMITSES,  G.J. 

"Optimal  beam  frequencies  by  the  finite  element  displacement  method" 

Int.  J.  Solids  & Structures,  9:415-429  (1973) 

IK121  KELLER,  J.B. 

"The  shape  of  the  strongest  column" 

Archive  for  Rational  Mechanics  and  Analysis,  5:275-235  (1960) 

(K131  KELLER,  J.B.  & NIORDSON,  F.I. 

"The  tallest  column" 

J.  Nath.  Mech.,  16:433-446  (1966) 

[K14J  KIUSALAAS,  J. 

"Optimal  design  of  structures  with  buckling  constraints" 

Int.  J.  Solids  & Structures,  9:363-878  (1973) 


221 


1 K151  KHOT,  N.S.,  VENKAYYA,  V.B.  & BEP.KE,  L. 

"OpCimum  design  of  composite  structures  with  stress  and  displacement  cons- 
traints" 

AIAA  paper  75-141,  AIAA  13th  Aerospace  Sciences  Meeting,  Pasadena,  Califor- 
nia (1975) 

(K16)  RICHER,  T.P.  & CHAO,  T.L. 

"Minimum  weight  design  of  stiffened  fibre  composite  cylinders" 

J.  Aircraft,  8:562-568  (1971) 

tK17J  KHOT,  N.S.,  VENKAYYA,  V.B.,  JOHNSON,  C.D,  & TISCHLER,  V.A. 

"Optimization  of  fiber  reinforced  composite  structures" 

Int.  J.  Solids  & Structures,  9:1225-1236  (1973) 

( K18]  KIUSALAAS,  J.  & SHAW,  R.C.J. 

"An  algorithm  for  optimal  structural  design  with  frequency  constraints" 

ASME  Energy  Technology  Conference,  Structural  Optimization  Methods  Session 
(ed.  VENKAYYA,  V.B.),  Houston,  p.  5-13  (1977) 


I LI  ] LUENBERGER,  D.G. 

"Introduction  to  linear  and  nonlinear  programming" 

Addlson-Wesley  (1973) 

( L2  I LEVENBERG,  K. 

"A  method  for  the  solution  of  certain  problems  in  least  squares" 

Quart.  Appl,  Math.,  2:164-168  (1942) 

[ L3  ) LANSING,  W.,  DIVYER,  W.,  EMERTON,  R.  & RANALLI,  E. 

"Application  of  fully  stressed  design  procedures  to  wing  and  empennage  struc- 
tures" 

J.  Aircraft,  3(9):683-688  (1971) 


( Ml  1 Me  CORMICK,  G.P. 

"Second  order  conditions  for  constrained  minima" 

SIAM  J.  Appl.  Math.,  15(3) :641-652  (1967) 

[ M2  J MARQUARDT,  D. 

"An  algorithm  for  least  squares  estimation  of  nonlinear  parameters" 

SIAM  J.  Appl.  Math.,  11:431-441  (1963) 

[M3  J MIELE,  A.,  HUANG,  H.Y.  & HEIDEMAN,  J.C. 

"Sequential  gradient  - restoration  algorithm  for  minimization  of  constrained 
functions  - Ordinary  and  conjugate  gradient  methods" 

J.  Optimization  Theory  Appl.,  4:213  (1969) 


( M4  J Me  CORMICK,  G.P. 

"Anti-zigzagging  by  bending" 

Managrat.  Sci.,  15:315-320  (1960) 

[ M5  1 Me  CORMICK,  G.P. 

"The  variable  reduet ion  method  for  nonlinear  programming" 

Managmt.  Sei.,  17:146-160  (1970) 

[ M6  1 !i0SES,  F. 

"Optimum  struetural  design  using  linear  programming" 

J.  of  the  Struetural  Division,  ASCE,  90,  ST6:39-104  (1964) 

( M7  J MOSES,  F.  & ONADA,  S. 

"Minimum  weight  design  of  struetures  with  applieation  to  elastie  grillages" 
Int.  J.  Num.  Methods  Engr.,  1(4): 311-331  (1969) 

[ M8  1 MOE,  J. 

"Penalty  funetion  methods" 

Optimum  Struetural  Design  - Theory  and  Applieations,  GALLAGHER,  R.ll.  A 
ZIENKIEWICZ,  O.C.  (ed.),  John  Wiley,  London,  chapitre  9 ; 141-177  (1973) 

[ M9  ] MARCAL,  P.V.  & GELLATLY,  R.A. 

"Applieation  of  the  treated  response  surfaee  technique  to  structural  optiml- 
zat ion" 

Proc.  of  the  Second  Conf.  on  Matrix  Methods  in  Structural  Mechanics,  AFFDL- 
TR-68-150  (1968) 

(MIOJ  MOE,  J. 

"Nonlinear  programming  in  structural  optimization" 

Course  on  Finite  Element  Methods,  Instituto  dl  Elaborazione  della  Informaziot- 
ne.  Pise  (1971) 

[Mill  maxwell,  C. 

Scientific  Papers,  Vol.  2 (1869)  - rcedite  par  Dover  (1952) 

[M121  MITCHELL,  A.G.M. 

"The  limits  of  economy  of  material  in  framed  structures" 

Phil.  Mag.  (Series  6),  8 (47) ;589-597  (1904) 

[ M13J  MELOSH,  R.J. 

"Structural  analysis,  frailty  evaluation  and  design.  Vol.  1 : Safer  tlteori- 
tical  basis" 

USAF  Technical  Report,  AFFDL-TR-70-15,  1 (1970) 

[M14]  MOSES,  F.  A KINSER,  D.E. 

"Optimum  structural  design  with  failure  probability  constraints" 

AIAA  J.,  5(6):1152-1158  (1967) 

[M15]  MOSES,  F. 

"Optimization  of  structures  with  reliability  constraints" 

Structural  Design  Applications  of  Mathematical  Programming  Techniques,  POPE, 
G.G.  A SCHMIT,  L.A.  (ed.),  AGAUDOGRAPH  149  ; 126-143  (1971) 

[M161  Me  CART,  B.R.,  HAUG,  E.J.  & STREETER,  T.D. 

"Optimal  design  of  structures  with  constraints  on  natural  frequency" 

AIAA  J.,  8(6): 1012-1019  (1970) 


■ 


223. 


I M17]  MIURA,  H.  & SCHMIT,  L.A. 

"Second  order  approximation  of  natural  frequency  constraints  In  structural 
synthesis" 

ASME  Energy  Technology  Conference,  Structural  Optimization  Methods  Session 
VENKAYYA,  V.B.  (ed.),  Houston,  p.  43-54  (1977) 


[ N1  1 NAGTEGAAL,  J.C. 

"A  new  approach  to  optimal  design  of  elastic  structures" 

Comp.  Methods  App.  Mech.  Eng.,  2(3):255-264  (1973) 

[N2  1 NIORDSON,  F.I. 

"On  the  optimal  design  of  a vibrating  beam" 

Quart.  Appl.  Hath.,  23,  1:47-53  (1964) 

[N3  1 NIORDSON,  F.I.  & PEDERSEN,  P. 

"A  review  of  optimal  structural  design" 

DCA^^M  Report  No.  31,  The  Technical  University  of  Denmark  (1972) 


(01  1 ORTEGA,  J.M.  & RHEINBOLDT,  W.C. 

"iterative  solution  of  nonlinear  equations  in  several  variables" 

Academic  Press,  New  York  (1970) 

(02  1 OLHOFF,  N. 

"A  survey  of  the  optimal  design  of  vibrating  structural  elements  - Part  I 
Theory" 

The  Shock  and  Vibration  Digest,  8(3):3-10  (1976) 

[03  ) OLHOFF,  N. 

"A  survey  of  the  optimal  design  of  vibrating  structural  elements  - Part  II 
applications" 

The  Shock  and  Vibration  Digest,  3(9);3-10  (1976) 

[04  ] OLHOFF,  N. 

"Optimal  design  of  vibrating  circular  plates" 

Int.  J.  Solids  & Structures,  6:139-156  (1970) 


224 


I PI  ) POWELL,  M.J.D. 

"Introduction  to  constrained  optimization" 

Numerical  Methods  for  Constrained  Optimization,  GILL,  P.E.  & MURRAY,  W. 
(ed.).  Academic  Press,  p.  1-28  (1974) 

( P2  J POWELL,  M.J.D. 

"An  efficient  method  for  finding  the  minimum  of  a function  of  several  varia- 
bles without  calculating  derivatives" 

Computer  J.,  7(2); 155-162  (1964) 

[ P3  J POLAR,  E. 

"Computational  methods  In  optimization" 

Academic  Press,  New  York  (1971) 

(P4  ] PEARSON,  J.D. 

"Variable  metric  methods  of  minimization" 

Computer  J,,  12:171-178  (I960) 

[P5  ] POWELL,  M.J.D. 

"Rank  one  methods  for  unconstrained  optimization" 

Integer  and  Nonlinear  Prograinnlng,  ABADIE,  J.  (ed.),  North-Holland,  Amster- 
dam, p.  139-156  (1970) 

( P6  ] POWELL,  M.J.D. 

"Some  properties  of  the  variable  metric  algorithm" 

Numerical  Methods  for  Nonlinear  Optimization,  LOOTSMA,  F.A.  (ed.).  Academic 
Press,  p.  1-17  (1972) 

(P7  ] POPE,  G.G. 

"Optimum  design  of  stressed  skin  structures" 

AIAA  J.,  11(11) :1545-1552  (1973) 

[ P8  J POPE,  G.G. 

"The  application  of  linear  programming  techniques  In  the  design  of  optimum 
structures" 

Proc.  AGARD  Symposium  on  Structural  Optimization,  Istambul,  AGARD-CP-36-70  : 
paper  3 (1970) 

(P9  J POPE,  G.G. 

"The  design  of  optimum  structures  of  specified  basis  configuration" 

Int.  J.  Mech.  Scl.,  10:251-265  (1968) 

(PlOl  PRAGER,  W. 

"Optimality  criteria  In  structural  design" 

Proc.  of  the  National  Academy  of  Sciences,  University  of  California,  61(3)  : 
794-796  (1968) 

[Pill  PRAGER,  W.  & TAYLOR,  J.E. 

"Problems  of  optimal  structural  design" 

J.  Appl.  Mech.,  ASME,  35:102-106  (1968) 

IP121  PRAGER,  W. 

"Conditions  for  structural  optimality" 

Computers  A Structures,  2(5);333-840  (1972) 

1P13]  PRAGER,  W.  A MARCAL,  P.V. 

"Optimality  criteria  In  structural  design" 

USAF  Technical  Report,  AFFDL-TR-70-166  (1971) 


225. 


( PRAGER,  W. 

"Optimal  design  of  statically  determinate  beams  for  given  deflection" 
Int.  J.  Mech.  Scl.,  13-893  (1971) 

( P151  PIERSON,  B.L. 

"A  survey  of  optimal  structural  design  under  dynamic  constraints" 

Int.  J.  Num.  Methods  Engr.,  4(4):491-A99  (1972) 


( R1  1 ROSEN,  J.B. 

"The  gradient  projection  method  for  nonlinear  programming  - Part  I : linear 
constraints" 

SIAM  J.,  8(1); 181-217  (1960) 

I R2  ] ROSEN,  J.B. 

"The  gradient  projection  method  for  nonlinear  programming  - Part  II  : non- 
linear constraints" 

SIAM  J.,  9:514-532  (1961) 

( R3  J RYAN,  D.M. 

"Penalty  and  barrier  functions" 

Numerical  Methods  for  Constrained  Optimization,  GILL,  P.E.  & TIURRAV,  U. 
(ed.).  Academic  Press,  p.  175-190  (1974) 

( R4  J REINSCHMlDT,  K.F.,  CORNELL,  A.C.  & BROTCHIE,  J.F. 

"Iterative  design  and  Structural  Optimization" 

J.  of  the  Structural  Division,  ASCE,  92,  ST6  : 281-318  (1966) 

( R5  J RAZANI,  R. 

"The  behaviour  of  the  fully  stressed  design  of  structures  and  Its  relation- 
ship to  minimum  weight  design" 

[ AIAA  J.,  3(12);2262-2268  (1965) 

I [R6  1 ROGERS,  L.C. 

'■  "Derivatives  of  eigenvalues  and  eigenvectors" 

I AIAA  J.,  6:2426-2429  (1968) 

[ [ R7  1 RUBIN,  C.P. 

j "Minimum  weight  design  of  complex  structures  subject  to  a frequency  cons- 

! traint" 

I AIAA  J.,  8(5);923-927  (1970) 

f [ R3  J ROMSTAD,  K.M. 

I "Optimal  design  for  dynamic  response" 

Symposium  on  Structural  Dynamics,  1,  paper  A2  (1970) 

r- 




226. 


I SI  1 STEWART,  III,  G.W. 

"A  modification  of  Davldon's  minimization  method  to  accept  difference  appro- 
ximations of  derivatives" 

J.  Assoc.  Comput*  Mach. ,14:72-33  (1%7) 

( S2  1 STEIN,  M.L. 

"Gradient  methods  In  the  solution  of  systems  of  linear  equations" 

Nat.  Bureau  Standards,  Report  52.7  (1951) 

I S3  1 SARGENT,  R.W.H. 

"Reduced-gradient  and  projection  methods  for  nonlinear  programming" 

1 Numerical  Methods  for  Constrained  Optimization,  GILL,  P.E.  & '^RJRR/W,  W. 

! (ed.).  Academic  Press,  p.  149-174  (1974) 

! [ S4  ] SARGENT,  R.W.H.  & MURTAGH,  B.A.H. 

"Projection  methods  for  nonlinear  programming" 

■ Math.  Prog.,  4:245-268  (1973) 

( S5  ] SHANLEY,  F.R. 

I "Weight-strength  analysis  of  aircraft  structures" 

I Dover,  New  York  (1960) 

i [ S6  J SCHMIT,  L.A. 

"Structural  synthesis  - 1959-1969  : A decade  of  progress" 

Recent  Advances  in  Matrix  Methods  of  Structural  Analysis  and  Design, 
GALLAGHER,  YAMADA  & ODEN  (ed.).  University  of  Alabama  Press  ; 565-634  (1071) 

1 S7  J SCHMIT,  L.A. 

"Structural  design  by  systematic  synthesis" 

Proc.  of  the  Second  National  Conference  on  Electronic  Computation,  Structure' 
Division,  ASCE,  New  York  ; 105-122  (1060) 

( S3  J SCHMIT,  L.A.  a,  FARSHI,  B. 

"Some  approximation  concepts  for  structural  synthesis" 

AIAA  J.,  12(5):692-699  (1974) 

(SO  J STORAASLI,  0.0.  & SOBIESZCZANSKI,  J. 

I "On  the  accuracy  of  the  Taylor  approximation  for  structural  resizing" 

AIAA  J.,  12(2):231-233  (1974) 

(SIO)  SCHMIT,  L.A.  & FOX,  R.L. 

"An  Integrated  approach  to  structural  synthesis  and  analysis" 

AIAA  J.,  3(6):1104-1112  (1965) 

(SllJ  SCHMIDT,  L.C. 

"Fully  stressed  design  of  elastic  redundant  trusses  under  .iltcrnatlve  load  in 
systems" 

Australian  J.  Applied  Science,  9:337-348  (1953) 

(S12J  SVED,  C.  & GINOS,  Z. 

"Structural  optimization  under  multiple  loading" 

Int.  J.  Mech.  Scl.,  10:803-805  (1968) 

(S13J  SHEU,  C.Y.  & SCHMIT,  L.A. 

"Minimum  design  of  elastic  redundant  trusses  under  multiple  static  loading 
conditions" 

AIAA  J.,  10:155-162  (1972) 


r f j 

227.  I 


[ Sl^  SHEU,  C.Y.  & PRAGER,  W. 

"Minimum  weight  design  with  piecewise  constant  specific  stiffness" 

J.  Optlm.  Theory  Appl.,  2(3):179-186  (1968) 

( SIS)  SHIELD,  R.T.  & PRAGER,  M. 

"Optimum  structural  design  for  given  deflection" 

J.  Appl.  Math.  Phys.  (ZAMP),  21:513-523  (1970) 

I S16J  SHEU,  C.Y.  & PRAGER,  W. 

"Recent  developments  In  optimal  structural  design" 

Appl.  Mech.  Reviews,  21 (10) :985-992  (1968) 

(SIT]  SOBIESZCZANSKI,  J.  & LOENDORF,  D. 

"A  mixed  optimization  method  for  automated  design  of  fuselage  structures" 

J.  Aircraft,  9(12) :805-811  (1972) 

( SIS]  SCHMIT,  L.A.  & MIURA,  H. 

"Approximation  concepts  for  efficient  structural  syrthesls" 

NASA  Contractor  Report,  NASA  CR-2552  (1976) 

( S19]  SAMCEF 

"Systeme  d'analyse  des  milieux  contlnus  par  olements  finis" 

LIAS,  Universite  de  Liege,  "Manuel  theorlque"  4 "rtanuel  d 'ut  II  isat  ion" 

( S20]  SCHMIT,  L.A.  & MIURA,  H. 

"A  new  structural  analysis/synthesis  capability  - ACCESS  I" 

AIAA  J.,  14(5):661-671  (1976) 

( S21J  SANDER,  G.  & FLEURY,  C. 

"A  mixed  method  in  structural  optimization" 

ASME  Energy  Technology  Conference,  Structural  Optimization  Mettiods  Session, 
VENKAYYA,  V.B.  (ed.),  Houston,  p.  79-93  (1977) 

[ S22]  SCHMIT,  L.A.  & IIALETT,  R.A. 

"Structural  synthesis  and  design  parameter  hierarchy" 

J.  of  the  Structural  Division,  ASCE,  89,  ST4:269-299  (1963) 

( S23j  SHEU,  C.Y. 

"Elastic  minimum  weight  design  for  specified  fundamental  frequency" 

Int.  J.  Solids  & Structures,  4:953-958  (1968) 

(S241  SIMITSES,  G.J.,  KAMAT,  M.P.  A SMITH,  C.V. 

"Strongest  column  by  the  finite  element  displacement  method" 

AIAA  J.,  11(9):1231-1232  (1973) 

(S251  SCHMIT,  L.A.  & FARSHI,  B. 

"Optimum  laminate  design  for  strength  and  stiffness" 

Int.  J.  Num.  Meth.  Engr.,  7:519-536  (1973) 

[S26]  SCILMIT,  L.A.  & FARSHI,  B. 

"Optimum  design  of  laminated  fibre  composite  plates" 

Int.  J.  Num.  Meth.  Engr.,  11:623-640  (1977) 

IS27J  SANDHU,  R.S. 

"A  survey  of  failure  theories  of  Isotropic  and  anisotropic  materials" 

USAF  Technical  Report,  AFFDL-TR-72-71  (1972) 


228. 


( T1  ] TRAUB,  J.F, 

"Iterative  methods  for  the  solution  of  equations" 

Prentice-Hallp  Englewoods  Cliffs,  New  Jersey  (1964) 

I T2  1 TOCHER,  J.L.  & KARNES,  R.N. 

"The  impact  of  automated  structural  optimization  on  actual  design" 

AIAA  Paper  ^3-361  (1971) 

( T3  ) TAIG,  I.C.  & KERR,  R.I. 

"Optimization  of  aircraft  structures  with  multiple  stiffness  requirements" 

AGARD  Second  Symposium  on  Structural  Optimization,  Milan,  AGARD-CP-123  ; 
paper  16  (1973) 

[ T4  1 TURNER,  M.J.  I 

"Design  of  minimum  mass  structures  with  specified  natural  frequencies"  I 

AIAA  J.,  5(3):406-412  (1967)  I 

( T5  1 TAYLOR,  J.E.  | 

"Minimum  mass  bar  for  axial  vibration  at  specified  natural  frequency"  1 

AIAA  J.,  5(10): 1911-1913  (1967)  I 

( T6  1 TAYLOR,  J.E.  j 

"Optimum  design  of  a vibrating  bar  with  specified  minimum  cross  section"  i 

AIAA  J.,  6(7);1379-1381  (1968) 

[T7  1 TADJliAKUSH,  1.  & KELLER,  J.B. 

"Strongest  columns  and  isoperlmetric  inequalities" 

J.  Appl.  Mech.,  ASME,  84:47-53  (1962) 

IT8  1 TAYLOR,  J.E. 

"The  strongest  column  - An  energy  approach" 

J.  Appl.  Mech.,  ASME,  34:486-487  (1967) 

I T9  ] TAYLOR,  J.E.  & LIU,  C.Y. 

"Optimal  design  of  columns" 

AIAA  J.,  6(8) : 1497-1502  (1968) 

[ TlOl  TSAI,  S.W. 

"Strength  characteristics  of  composite  materials" 

NASA  Contractor  Report,  NASA-CR-224  (1965) 


[VI  1 VANDERPLAATS,  G.N.  & MOSES,  F. 

"Structural  optimization  by  methods  of  feasible  directions" 
Computers  & Structures,  3(4) :739-755  (1973) 

IV2  1 VENKAYYA,  V.B.,  KHOT,  N.S.  & REDDY,  V.S. 

"Energy  distribution  in  an  optimum  structural  design" 

USAF  Technical  Report,  AFFDL-TR-63-156  (1969) 

[V3  J VENKAYYA,  V.B. 

"Design  of  optimum  structures" 

Computers  4 Structures,  1:265-309  (1971) 


P 

i 


229. 


[ 

r 

i 

n 

i 

t 

f 

r 

k 

I 

[ 

f 

t 


( V4  I VENKAYYA,  V.B.,  KHOT,  N.S.  & BERKE,  L. 

"Application  of  optimality  criteria  approaches  to  automated  design  of  large 
practical  structures" 

AGARD  Second  Symposium  on  Structural  Optimization,  Milan,  AOARI)-C?-123  : pa- 
per 3 (1973) 

( V5  J VENKAYYA,  V.B.  & KHOT,  N.S. 

"Design  of  optimum  structures  to  impulse  type  loading" 

AIAA/ASME/SAE  15th  Struct.,  Struct.  Dyn.  and  Matls.  Conf.,  Las  Vegas  (1974) 

( V6  ] VENKAYYA,  V.B.,  KHOT,  N.S.,  TISCHLER,  V.A.  & TAYLOR,  R.K. 

"Design  of  optimum  structures  for  dynamic  loads" 

Third  Air  Force  Conference  on  Matrix  Methods  in  Structural  Mechanics,  WPAFli, 
Ohio  (1971) 

[ V7  ] VENKAYYA,  V.B. 

Private  communication  (1977) 


(W1  J WOLFE,  P. 

"Methods  for  nonlinear  constraints" 

Nonlinear  Programming,  ABADIE,  J.  (ed.),  North-Holland,  Amsterdam,  p.  120- 
131  (1967) 

[W2  I WOLFE,  P. 

"A  duality  theorem  for  nonlinear  programming" 

Quart.  Appl.  Math.,  19:239-244  (1963) 

(W3  1 WILDE,  D.J.  & BEIGHTLER,  C.S. 

"Foundations  of  optimization" 

Prentice-Hall,  Englewoods  Cliffs,  New  Jersey  (1967) 

{W4  J WILSON,  R.B. 

"A  slmplicial  algorithm  for  concave  programming" 

Doctoral  dissertation.  Harvard  University,  Boston  (1963) 

[W5  1 WITTRICK,  W.H. 

"Rates  of  change  of  eigenvalues  with  reference  to  buckling  and  vibration  pro- 
blems" 

J.  of  the  Royal  Aeronautical  Society,  66:590  (1962) 

(W6  J WADDOUPS,  M.E.,  Me  CULLERS,  L.A.,  OLSEN,  F.O.  & ASHTON,  J.E. 

"Structural  synthesis  of  anisotropic  plates" 

AIAA/ASME  11th  Struct.,  Struct.  Dyn.  and  Matls.  Conf.,  Denver  (1970) 


230. 


[Y1  ] YOUNG,  J.W.  & CHRISTIANSEN,  H.N. 

"Synthesis  of  a space  truss  based  on  dynamic  criteria" 

J,  of  the  Structural  Division,  ASCE,  92,  ST6;425-442  (1966) 


[1 


I 


4 

( Z1  ) ZANGWILL,  W.  | 

"Minimizing  a function  without  calculating  derivatives" 

Computer  J.,  10:293-296  (1967) 

[ Z2  J ZOUTENDIJK,  G. 

"Methods  of  feasible  directions" 

Elsevier,  Amsterdam  (1960) 

I Z3  ) ZOUTENDIJK,  G. 

"Nonlinear  programming,  computational  methods" 

Integer  and  Nonlinear  Programming,  ABADIE,  J.  (ed.),  Nortli-Holland,  Amster- 
dam, p.  37-36  (1970) 

( Z4  ) ZARGHAMEE,  M.S. 

"Optimum  frequency  of  structures" 

AIAA  J.,  6(4):749-750  (1968) 

( Z5  I ZARGHAMEE,  M.S. 

"Minimum  \jelght  design  with  stability  constraint" 

J.  of  the  Structural  Division,  ASCE,  96,  ST3:1697-1710  (1970) 


